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ABSTRACT 


This is Volume 1 of three volumes which report the results of a strapdown calibration 
and alignment study performed by the Univac Federal Systems Division for the Guidance 
Laboratory of NASA/ERC. 

This study develops techniques to accomplish laboratory calibration and alignment of 
a strapdown inertial sensing unit (ISU) being configured by NASA/ERC. Calibration 
is accomplished by measuring specific input environments and using the relationship 
of known kinematic input to sensor outputs, to determine the constants of the sensor 
models. The environments used consist of inputs from the earth angular rate, the 
normal reaction force of gravity, and the angular rotation imposed by a test fixture m 
some cases. Techniques are also developed to accomplish alignment by three methods. 
First, Mirror Alignment employs autoeollimators to measure the earth orientation of 
the normals to two mirrors mounted on the ISU. Second, Level Alignment uses an 
autocollimator to measure the azimuth of the normal to one ISU mirror and accelerom- 
eter measurements to determine the orientation of local vertical with respect to the 
body axes. Third, Gyrocompass Alignment determines earth alignment of the ISU by 
gyro and accelerometer measurement of the earth rate and gravity normal force vectors. 

The three volumes of this study are composed as follows: 


• Volume 1 - Development Document. This volume contains the detailed develop- 
ment of the calibration and alignment techniques. The development is presented 
as a rigorous systems engineering task and a step by step development of 
specific solutions is presented. 

• Volume 2 — Procedural and Parametric Trade-off Analyses Document. This 
volume contains the detailed trade-off studies supporting the developments 
given in Volume 1. 

• Volume 3 — Laboratory Procedures Manual. In Volume 3 the implementation 
of the selected procedures is presented. The laboratory procedures are 
presented by use of both detailed step-by-step check sheets and schematic 
representations of the laboratory depicting the entire process at each major 
step m the procedure. The equations to be programmed in the implementation 
of the procedures are contained in this volume. 
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GLOSSARY 


As an aid to understanding the symbolism, we present the following rules of notation 


• Wherever possible symbols will be used which suggest the name of the 
parameter involved. 

• Lower case subscripts are used almost exclusively for indexing over several 
items of the same kind. Examples are the indexes used to identify the three 
gyros, the three accelerometers, the two pulse trains of each accelerometer, 
the two clock scale factors, etc. 

• Lowercase superscripts are used to index over different positions. 

• Uppercase superscripts and subscripts will be used to distinguish between 
parameters of the same kind. For example, T is used to identify a 

BE 

transformation matrix. Lettered superscripts such as BE in T identify 
the particular transformation. 

• An underline will identify a vector. 

• Unit vectors are used to identify lines in space such as instrument axes 
and the axes of all frames of reference. 

• Components of any vector along with any axis is indicated by a dot product 
of that vector with the unit vector along the axis of interest. 

• The Greek sigma (£) will be used for summations. Where the limits of 
summation are clear from the context, they will not be indicated with the 
symbol. 

• The Greek A is always used to indicate a difference. 

• S <ji and C 0 are sometimes used to identify the sine and cosine of the angle 0. 

• A triple line symbol (=) will be used for definitions. 

• A superior 11 " denotes a prior estimate of the quantity. 

• A superior" A "denotes an estimate of the quantity from the estimation routine. 
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SECTION 1 
INTRODUCTION 


This document, m conjunction with two other volumes, describes the achievements of a 
six month study conducted for the. 


Guidance Laboratory 
Electronics Research Center 
National Aeronautics and Space Administration 
Cambridge, Massachusetts 


by the: 


Aerospace Systems Analysis Department 
Umvac Federal Systems Division 
Saint Paul, Minnesota 
A Division of Sperry Rand Corporation 

The purpose of the study is to develop techniques and outline procedures for the labora- 
tory calibration and alignment of a strapdown inertial sensing unit. This document, 
Volume 1. presents a detailed analysis of the calibration and alignment problem and 
develops a specific solution. The nucleus of the study output is the contents of this docu- 
ment. The Procedural and Parametric Trade-off Analyses, Volume 2, is a set of 
addendums which serve to justify decisions made and conclusions reached m the develop- 
ment of specific calibration and alignment techniques. Reference is made to the contents 
of the trade-off document throughout Sections 4 and 5 of this document. The Laboratory 
Procedures Manual. Volume 3, describes specific procedures for an operational im- 
plementation of the solutions obtained in Volume 1. It is an extension of the results of 
Volume 1 into' an operational laboratory situation. The last subsections of Sections 4 
and 5 of this document (Volume 1) form the interface between the study developments and 
the specific procedures found m Volume 3. 

tV 

At the time of this writing, the Guidance Laboratory of NASA/ERC is in the process of 
configuring a strapdown inertial sensing unit which they will use to evaluate many ad- 
vanced concepts. By integrating this ISU with a system computer, they will attain a 
flexible system level research tool for testing analytical concepts, system design con- 
cepts and fabrication concepts. In parallel with the development of the Guidance and 
Navigation System, a laboratory facility is being designed which will contain all of the 
test equipment necessary for conducting the experiments on the strapdown G and N System. 
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Among the more important experiments to be conducted are those which determine the 
feasibility of easily and precisely calibrating and aligning the sensor package m an 
operational laboratory situation. Before such experiments can be conducted, the cali- 
bration and alignment techniques must be developed and definitively documented. The 
Guidance Laboratory contracted Umvac's Aerospace Systems Analysis Department to 
develop and document those techniques. The specific tasks which Univac was contracted 
to accomplish are as follows: 

• To specify mathematical models for the system sensors (gyros and accelerom- 
eters). 

• To define the mathematical description of the sensor package. 

• To develop techniques for the determination of all calibration constants. 

• To develop three techniques for initializing the alignment of the ISU. The three 
techniques involve the use of 

A Optical measurements only 

A Accelerometer measurements for level, and an optical azimuth measurement 
A Accelerometer and gyro measurements only. 

• To accomplish specified trade-off analysis on all calibration and alignment 
techniques. 

• To specify all equations and procedures for the accomplishment of a calibration 
and alignment m the ERC Laboratory. 

• To document, m three volumes, the calibration and alignment developments, 
trade-offs, and procedures. 

The satisfaction of the first four items is accomplished m this Development Document. 

The trade-offs are described in the Procedural and Parametric Trade-off Analysis 
Document and the procedures are outlined m the Laboratory Procedures Man ual . 

# 

The presentation of the calibration and alignment developments in this document, 

Volume 1, is divided into five sections. The purpose of the introduction is to briefly 
state the study problem (accomplished m the above listing of seven items) and to describe 
the developments contained in Volume 1, The purpose of Section 2 is to delineate the 
calibration and alignment requirements. Section 3 presents a system description of 
calibration and alignment with emphasis on the laboratory environment. The specific 
calibration and alignment techniques are then developed in Sections 4 and 5, respectively. 
As an introduction to the scope of this document, the following paragraphs outline the 
developments in these sections. 

The calibration and alignment study tasks have been only generally stated in the preceding 
paragraphs. Before the technique developments can be described, the specific engineering 
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and mathematical requirements of calibration and alignment must be stated. In Section 2 
we accomplish the detailed specification of those requirements. A statement of the cali- 
bration and alignment requirements will be simply presented as a list of parameters to 
be determined in the laboratory. As a lead-in to that listing, Section 2 shows how the 
requirements tie m to the larger system problem of navigation. We accomplish this by 
presenting the general definitions of calibration and alignment as the determination of 
constants required in an operational navigation loop. The mathematics of portions of the 
navigation loop are delineated so that calibration and alignment can be specifically defined 
as the determination of constants contained within the mathematics. 

After specifying the calibration and alignment requirements m Section 2, Section 3 directs 
our attention to the laboratory environment m which the calibration and alignment is to be 
accomplished. As an introduction to the environment, we present m the first subsection 
of Section 3 functional system descriptions of both calibration and alignment. The func- 
tional description of the ERC laboratory calibration is presented in comparison with what 
we call an Ideal Calibration. The comparison of the ideal with the actual ERC laboratory 
calibration serves to illustrate those compromises necessary in the development of a test 
laboratory. The functional description of alignment presents those separate operations 
required m an operational alignment. Three functional diagrams are presented m 
Section 3, one for each of the three alternative alignment techniques. All functional 
descriptions serve to define those measurements, other than inertial instrument measure- 
ments, which are required to accomplish the calibration or alignment. The additional 
measurements correspond to an independent measure of the kinematic environment. The 
determination of those additional measurements is the subject of the second subsection of 
Section 3. Section 3 is concluded with a brief description of the hardware available in the 
laboratory, and the interfaces between those pieces of hardware. 

The calibration technique developments in Section 4 are directed toward specifying the 
details of the calibration functions which are generally defined m Section 3. The basis of 
calibration is presented in Section 3 as the input of environment and inertial instrument 
measurements into computations which are a function of those measurements and the un- 
known calibration constants. The general equations from which the computations are 
evolved are developed in the initial subsection of Section 4. Those general equations are 
developed by introducing the parameters which identify the laboratory kinematic en- 
vironment and the ISU geometry into the inertial instrument mathematical models. Subse- 
quent to the development of the general equations, particular choices of test table orienta- 
tion are used to define the "Positions" to be used for the determination of all calibration 
constants. The chosen positions are shown to produce significant reductions m the com- 
plexity of the general calibration equations. With the aid of these reductions it is 
possible to solve for the calibration constants by a series of relatively simple experiments. 
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In the third subsection of Section 4 the calibration computations are tabulated. The 
quantization and instrument and environment noise considerations are described in con- 
junction with the tabulation. The fourth subsection of Section 4 describes those laboratory 
activities required prior to the actual calibration. All such activities are related to either 
the survey of the location of the ISU relative to the test table, or the compensations for the 
small low frequency motion of the test table base. The last subsection of Section 4 forms 
the tie between this Development Document and the details of calibration implementation 
presented in the Laboratory Procedures Manual, Volume 3» In Volume 1 the implementa- 
tion of calibration is only briefly described, the details being left as the subject of 
Volume 3. 

The alignment techniques developed in Section 5 expand the functional descriptions of 
alignment as presented in Section 3 into a set of alignment techniques. Alignment is 
broken into three separate routines, preprocessing of sensor outputs, the application of 
chosen estimation procedures to the preprocessed outputs, and calculation of alignment 
matrices from the estimated values. Since the preprocessing and alignment matrix 
calculations are developed in Section 2, the major emphasis m Section 5 is centered on 
the estimation problem. 

Before describing the development of an estimation technique, the basic functional re- 
quirements and the preprocessing computations are presented, respectively, m the first 
two subsections. The third subsection describes a detailed development of models for 
the environmental disturbance and sensor noise. The next two subsections are then de- 
voted to the development of two approaches for estimation m Level Alignment and Gyro- 
compass Alignment. The first approach develops a procedure for estimating average 
values of the gravity and earth rate vectors, while the second approach leads to estimates 
of instantaneous values of these vectors. Estimation techniques are developed using three 
basic statistical procedures: simple average, least squares, and posterior mean. From 
these estimates the average and instantaneous values of the alignment matrices are then 
obtained. The last subsection of Section 5 describes explicit equations for the recom- 
mended alignment techniques, and ties the results of Section 5 to the procedural details 
of alignment described m the Laboratory Procedures Manual. 
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SECTION 2 

CALIBRATION AND ALIGNMENT REQUIREMENTS 


The purpose of this study is to determine a procedure for the calibration of the NASA/EEC 
strapdown inertial sensing unit (ISU) and to delineate three operational laboratory tech- 
niques for the initial alignment of the same inertial sensor unit. Clearly, the initial task 
in this, or any study, is to carefully describe the problem as a specifically defined study 
task. This we propose to do m this section of our report. 

The key words in the above general statement of the study purpose are the words "calibra- 
tion" and "alignment". The first activity in this section will be to develop (in Section 2. 1) 
the definitions of those key words. Our approach to the development of those definitions 
is to present a description of an operational navigation loop and, as a conclusion to that 
description, to present calibration and alignment as the determination of constants re- 
quired as inputs to the navigation loop. There are alternative approaches to the defini- 
tions of these terms but we feel our approach is optimum in that it clarifies the necessary 
relationship between the calibration and alignment problem and the larger system problem 
of inertially navigating a propelled vehicle. 

Subsequent to the navigation- system definition of calibration and alignment we will, in 
Sections 2. 2 and 2. 3, describe the calibration and alignment requirements as they relate 
to the ERC strapdown inertial sensing unit. Section 2. 2 describes the calibration re- 
quirements, and Section 2. 3 describes the alignment requirements. 

The development of the calibration requirements m Section 2. 2 will be directed toward 
the tabulation of the instrument constants and instrument- to-body-axes transformation 
matrix constants which are necessary m an operational navigation loop. The first 
activity in that section will be the description of the geometry of the ERC ISU. This will 
be followed by a description of the inertial instruments contained m that ISU. The m- 
strument-to-body-axes transformation matrices will then be described. All of the 
described equipment and geometry will then be used to develop the navigation loop "Pre- 
processing Computations". Finally, the constants in the Preprocessing Computations 
will be defined as the constants to be obtained m calibration. 

In Section 2. 3 the alignment requirements will be described as the real-time measurement 
of the ISU fixed or earth-fixed coordinates of two vectors. Three alternative choices of 
these two vectors will be presented. The geometry of alignment will also be presented. 
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The final presentation will be the specific alignment matrix mathematics corresponding to 
the three alternative alignment techniques. 


2. 1 DEFINITION OF CALIBRATION AND ALIGNMENT 

The necessity for a calibration and alignment of an inertial sensing unit is directly related 
to the use of the ISU in a vehicle guidance system. More specifically, calibration and 
alignment requirements are related to the necessity for a real-time transformation of the 
ISU instrument outputs into a best estimate of a vehicle's velocity and position. In this 
section we will, from a discussion of real-time inertial guidance activities, define the 
general calibration and alignment requirements. 

The functional diagram shown m Figure 2-1 serves as a description of the initial activities 
m a real-time navigation loop. That diagram will be the focal point of our attention for 
the remainder of this subsection. Figure 2-1 shows only that portion of the navigation 
loop which transforms instrument outputs into estimates of velocity and position. (The 
remaining portions of the loop are the guidance logic, automatic control, and dynamic 
response which are not shown. ) 

The input to the ISU is the kinematic environment of the vehicle and ISU as represented 
by the applied acceleration a and angular velocity w. The outputs of the ISU are (usually) 
sequences of pulse counts taken over small intervals of time. These outputs are the in- 
puts to the computer. The computer's immediate task is to convert those measurements 
into a knowledge of velocity and position. The velocity (v) and position (r) must be rep- 
resented as components (v^D^) and (r*D^) m the frame (D^) in which one chooses to 
navigate (D^ represents a triple, k = 1, 2, 3, of unit vectors directed along the orthogonal 
navigation axes). 

The initial activity m the conversion to velocity and position is the transformation of the 
pulse counts into estimates of the integrals of the instrument-axes components of applied 
acceleration and angular velocity. The instrument axes are represented by the triads 
A^ and of (in general) nonorthogonal unit vectors directed along the input axes of the 
accelerometers and gyros, respectively. 

The second activity, in the conversion to velocity and position, is the transformation of 
the integrals of the instrument-axes components of applied acceleration and angular 
velocity into integrals of body-axes components. The body axes (B^) are a triad of 
orthogonal unit vectors which are fixed to the ISU. These body axes can be defined m 
various ways. They can be defined by use of any two of the instrument axes or they can, 


2-2 


















as they are in this study, be defined by use of two mirror normals. (The manner in 
which this definition is accomplished for the EEC system will be found m Section 2. 2. ) 

The final activity m the conversion of instrument outputs into navigation- axes components 
of velocity and position begins with the input of the body-axis integrals of applied accel- 
eration and angular velocity to the translational and rotational differential equations of 
motion. The numerical solution of those differential equations yields the desired velocity 
and position. The solution of the rotational differential equations serves to transform 
the argument of the translational differential equations into navigation- axes components. 
The output of the translational differential equations solution is then the desired com- 
ponents of velocity and position. 

It is noted that various constants are required from computer memory as inputs into all 
routines. The initial routine requires those instrument constants which scale and 
correct the instrument outputs. The second routine requires the nonorthogonal three-by- 

A £1 

three matrices, Q and Q , which transform the integrals of the instrument-axes 
components into the integrals of the body-axes components. The third routine, being 
the solution of differential equations, requires initial conditions. The initial condition 
for the rotational differential equation solution is an initial body- to- navigation- axes, 
transformation matrix. The initial condition of the translational differential equation 
solution is an initial knowledge of navigation- axes components of velocity and position. 

A knowledge of all of these constants is required prior to any operational use of the ISU» 
The development of the numerical values of these constants can be divided into three 
separate problems; and the statement of two of these problems can be used as a defini- 
tion of calibration and alignment. 

The problem of determining the instrument constants used in the first routine, and the 
A G 

Q and Q matrices used in the second routine, will be considered m this report as the 
problem of calibration. The problem of determining the initial body-to-navigation-axes 
transformation matrix will be referred to as the alignment problem. The remaining 
problem of initializing velocity and position is an operational problem, which is not 
within the scope of this work. 

We will extend these definitions to the subject EEC strapdown ISU. Specifically, we will 
delineate more detailed definitions m terms of the geometry and instruments charac- 
terizing the ERC system. Section 2. 2 will treat calibration, and 2. 3 alignment. 
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2.2 CALIBRATION REQUIREMENTS 


It was seen in the preceding subsection that calibration is defined as the determination 
of those instrument constants and constant matrix transformations required m the 
transformation of instrument outputs into integrals of body-axes components of applied 
acceleration and angular velocity. In this section we will describe the equations in the 
navigation routines which utilize the calibration constants. We will specify those 
equations with the assumption that the ERC strapdown ISU is the subject sensor unit. 
From that description, we can then specifically describe the calibration requirements 
as the determination of the constants contained within those navigation routines. 

The desired equations are directly deducible from the geometry of the ISU and the 

mathematical models of the instruments. We therefore begin the presentation in this 

section by describing the geometry of the strapdown ISU, followed by a description of 

the accelerometers and gyros contained within the ISU. Following those descriptions 
A G 

we will define the Q and Q matrices. Next we will employ all of this information to 
develop the desired equations; and finally we will utilize those equations m the tabulation 
of the required calibration constants. 

2.2.1 ISU Geometry 

The ERC ISU is a strapdown sensing unit containing. 

• Three vibrating- string accelerometers 

• Three single degree of freedom gyros 

• One mirror cube 

• Associated structural and electronic devices. 

The strapdown ISU has been specified such that the accelerometer input axes (A^.), the 
gyro input axes (G ^), and three mirror normals (M^) are nominally orthogonal and 
nominally aligned. In implementing the specification, there will naturally be deviations 
of small angles between the supposedly aligned instrument and mirror axes. In 
Figure 2-2 an exaggerated representation of those deviations from nominal is shown. It 
will be assumed in this study that the cosines of the angles between supposedly aligned 
vectors are equal to one and the cosines of the angles between supposedly orthogonal 
vectors are equal to small first order numbers. 

Additional geometry required m subsequent developments is the nominal location of the 
gyro output (O^) and spin (S^) axes relative to the input axes already described. Those 
nominal locations are shown in Figure 2-3. 
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GYRO ONE 


Figure 2-3. Nominal Orientation Between Gyros 
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It is always convenient, m inertial navigation, to define, from the geometry of the ISU, a 
set of orthogonal unit vectors which represent a common "body set of axes” to which all 
accelerometer and gyro outputs can be referred. The definition of those body axes is 
usually arbitrary; that is, any two ISU fixed vectors can be used. For the purpose of this 
study we will utilize two mirror normals. The body axes for the subject ISU are defined 
by 




B 2 = (M 1 x Mg) x Mj/ 

B 3 = (M x xM 2 )/|M t x M 2 


x Mg 


This definition is shown schematically m Figure 2-4. 
2. 2. 2 Accelerometer Model 


A schematic of the accelerometer is shown m Figure 2-5. The accelerometer consists 
of two masses separated by a spring and supported for centering purposes by two strings 
(S^ and Sg) and ligaments normal to and Sg. When the accelerometer is at rest or 
moving with constant velocity, the sum of forces acting on the masses is zero. When the 
instrument is accelerated, the sum of forces will adjust to cause the masses to move with 
the same acceleration. Strings and Sg will change in tension as a function of the com- 
ponent of acceleration along the strings. (This direction is the sensitive axis of the in- 
strument. ) Since the resonant frequency of a vibrating- string is a function of its tension, 
the frequency of strings S^ and Sg may be read and converted to acceleration along the 
sensitive axis. 


The math model of the accelerometer is presented on Chart 2-1. The outputs from the 
accelerometer are the pulse counts, and Ng, representing the number of zero 
crossings from strings Sj and Sg in the time interval t^-t^. Since the counting process 
can start and terminate at a fixed time for any sample, a quantization error (represented 
by Eq) of up to two counts (one per string) may occur. A fixed bias (D Q ) is assumed. 

The entire output is multiplied by the scale factor D^. The second and third order co- 
efficients (Dg and Dg) are extremely small. 

In developing the model (see Appendix A) several assumptions were made. The most 
critical are 


1. The accelerometer has negligible instrument noise. 

2. The effects of terms higher than the third order are negligible. 


2-8 







Model: ARMA D4E Vibrating String Accelerometer 

Axis: A is a unit vector directed along strings S^ and Sg 

(the sensitive axis) 

Figure 2-5. A Schematic Diagram of the Accelerometer 


2-10 










CHART 2-1 


THE FUNDAMENTAL ACCELEROMETER MODEL 


THE ACCELEROMETER MODEL IS: 



a is the acceleration applied to the accelerometer 

t s=t£ t^ is the time interval over which a is measured 

A is a unit vector directed along the input axis of the accelerometer 

and Ng are the number of zero crossings detected m ^ t ^ t^ 
from both strings of the accelerometer 

Eq is the instrument quantization error due to the fact that t a and t^ 
do not correspond to zero crossings 

is the accelerometer scale factor 
Dq is the accelerometer bias 
Dg is the second order coefficient 
Dg is the third order coefficient 
fg and f^ are string frequencies in pulses/second 


2-11 









3. There are no cross coupling effects. 

4. The strings are cohnear. 

2. 2. 3 Gyro Model 

The gyroscopes used m the ERC strapdown ISU are Honeywell GG334A single degree of 
freedom, pulse rebalance gyros. The gyros contain a gimballecl rotor as shown m 
Figure 2-6. The rotor spins at a high angular rate. The gimbal is restrained by the 
gimbal bearings to rotate with respect to the case about the O axis only as shown in the 
figure. Any angular motion of the gyro case about the input axis, G, will generate a 
gyroscopic torque that tends to rotate the gimbal about O. A signal generator measures 
the gimbal deflection. The deflection is compared at a 3. 6 KHz rate with two equal 
thresholds of opposite sign and a positive, negative, or zero pulse is generated, based 
on the results of the comparison. This signal is sampled by the readout electronics and 
fed to a torque generator where a torque pulse is generated to offset the deflection. 

The model of the gyro is given on the accompanying chart. (See Appendix B for a denva- 
N 

tion. ) £ 6 v is the net count of positive and negative rebalance torques. A<i> is the scale 
k-1 * t N 

factor of the instrument. The term (I (w* G) dt is the desired information from the in- 

**0 — 

strument. and is equal to the integral of the angular velocity component along the sensi- 
tive axis. R is a fixed bias term. The three terms with coefficients, Bp Bq and B g 
are due to the fact that the center of force of the gimbal support differs from the gimbal 
center of mass, causing a torque proportional to acceleration (mass unbalance effect). 
Terms with coefficients Cjp Cg g , Cj g , Cq S and Cjq arise because of the deformations 
of the gimbal, caused by acceleration forces that produce mass unbalance effects. The 
term with Qjj coefficient is due to scale factor nonlmearities. The Qjg term is due to 
the differences of moments of inertia about S and O. The term containing J is the 
effect of dynamic coupling because of finite gimbal inertia. 

2. 2. 4 Q Matrices 

A G 

In Section 2. 1 we defined the Q and Q matrices as those constant matrices which 
transform the integral of the instrument- axes components of applied acceleration and 
angular velocity into the integrals of the body-axes components of the same vectors. 

In this subsection we will specifically define those matrices. 

A G 

First, the Q and Q matrices, as suggested by the superscripts, transform, re- 
spectively, the integrals of the accelerometer-axes components and the integrals of the 

A P 

gyro-axes components. Second, Q and Q , being constant matrices, transform all 
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Model: Honeywell GG 334A single -degree -of -freedom, pulse 
rebalance gyroscope. 

Axes: S is a unit vector along the spin axis of the rotor. 

O is a unit vector along the output axis as defined by 
the gimbal. 

G = O x S is the sensitive axis of the gyro. 


Figure 2-6. A Schematic Diagram of the Gyro 
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CHART 2-2 


THE FUNDAMENTAL GYHO MODEL 


THE GYRO MODEL IS: 


E 6 k - / N (u>.G) dt + / N R+Bj(a-G) + B 0 (a*0) + Bg(a-S) + C n (a*G) 2 + C sg (a.S) 2 

- k=1 J Vl I— 


+ Cjg(a G) (a S) + C ns {a-0) (a-S) + C in (a*G) (a-O) 




? d 

Qjj( as* G) + (o)'S) + J ^ (tj>* O) dt + An + Eq 


u. is the angular velocity applied to the gyro 
a is the acceleration applied to the gyro 

tp s t ^ t^ is the time interval over which a and u- are measured 

t N - t Q = Nr, where N is an integer, andr is the gyro sampling 
period 

S is a unit vector along the spin axis of the rotor 

O is a unit vector directed along the output axis as defined by the 
gimbal 

G is a unit vector along OxS (that is, the sensitive axis of the gyro) 

6 k is the kth gyro pulse, equal to +1, -1, or 0 for positive, negative, 
or no pulse 

A<f> is the gyro scale factor 
R is the gyro bias 

Bj Bq and B g are the gyro unbalance coefficients 

C II ^SS C IS ^OS and ^10 are tke gyr0 COm P li ance coefficients 

Qjg and Qjj are dynamic coupling coefficients due to gimbal deflection 
and scale factor nonlinearity, respectively 

J is the angular rate coefficient 

An is the effect of gyro noise over the Ctg, t N ] interval 
Eq is the gyro quantization error 
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triples of vector components between the frames, not just the integrals of those vector 

components. (That is, the Q matrices can be taken in and out of the integral at willo ) 

A G 

Therefore, we can say that Q and Q are defined by: 


-1 




Ail 

§2 

= 

q a 


A 2 

A 


_ _ 


_-3_ 

’5f 




V 

§2 

= 

q g 


—2 

*3 




-3 


It will be seen m Section 4 tnat calibration determines not Q A and Q G but (Q A )~'*‘ and 
C —1 

(Q ) • We therefore need to deduce the matrices from their inverses. From the 

A G 

geometry presented m Section 2. 2„ 1 it is seen that the Q and Q matrices are approxi- 
mately identity matrices. This fact makes the deduction of the matrices from their in- 
verses quite simple. In the accompanying chart that deduction is presented. Note, in 
Chart 2-3, that the inverses appear, at quick glance, to be orthogonal (that is the ele- 
ments are "direction cosines"). This apparent orthogonality results from being 
orthogonal. However and are not, m general, orthogonal, and therefore the in- 
verse matrices are also not orthogonal. 

A G 

In Appendix C alternate forms of the Q and Q matrices are presented. Those forms 
are functions of the~separation-angles between the unit-vectors contained within the ISU. 
Even though we will not specifically present techniques for finding separation-angles, the 
reader may be interested in those forms for the purpose of deducing separation angles 
from the calibration-determined (Q A )"^ and (Q*"*) - '*" elements. 

2.2.5 Preprocessing Computations 

In this section we will show how the ISU geometry and mathematical models lead to the 
specific equations found in the initial computational routines of a navigation loop. Those 
equations, which we call the Preprocessing Computations, include all of the constants 
which must be determined during a laboratory calibration. 

Referring to the flow diagram presented m Figure 2-1 we see that the initial routine in 
the navigation loop is the transformation of the instrument outputs into a knowledge of the 
integrals of the instrument-axes components of applied acceleration and angular velocity. 
Referring to the instrument models, we see that the models represent functional rela- 
tionships between the inputs and outputs of the initial navigation routine. (This statement, 
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CHART 2-3 


Determination of Q A and Q G from (Q A ) -1 and (Q G ) 1 


• The calibration determines not Q A and Q G but (Q A ) -1 and (Q G ) _1 . 

• The inverse matrices have the form. 

(Aj_ . B 1 ) (Aj_ • B 2 ) (Aj_ • Bg) 

(Q A ) -1 - (A 2 - Bj) (A 2 . B 2 ) (A 2 - Bg) 

(Ag • Bj) (Ag • Bg) (Ag • Bg) 


(G x • B x ) (G x ' B 2 ) (G 1 • Bg) 

(Q G ) -i = (g 2 - Bj) (g 2 -b 2 ) (G 2 -B 3 ) 

< g 3-=2> (5s- 

Because of the excellent mechanical specifications on the strapdown ISU, each 
of the above matrices will have ones on the diagonals and first order small 
quantities on the off-diagonal. That is, each matrix can be written as. 


where I is the identity matrix and E is a small off-diagonal matrix. 

The inverse of (I + E) (to first order) is (I - E). 

A G 

The Q and Q matrices can, therefore, be written as: 


1 -(A^Bg) -(Aj-Bg) 

-(A 2 *Bj) 1 -(Ag* Bg) 

-(Ag - §i) -(Ag*B 2 ) 1 


= 


( g 2 * *l> 


“(Gj-Bg) -(Gj’Bg) 
1 * — 3 ) 


-(Qa - Bl) -(Gg-Bjj) 


where all elements within the matrices are found m calibration. 





of course, assumes that quantization and noise are neglected*) The mathematical model, 
functional relationships are not, however, explicitly in the form: routine output = f (rou- 
tine input). To the contrary, the relationships are the inverse: routine input = g (routine 
output). The desired navigation equations must therefore result from an inversion of the 
instrument mathematical models. 

The inversion of the models is quite a simple matter. This is due to the fact that each 
instrument is designed to be a linear instrument, therefore, all "nonlinear" terms are 
the result of design deficiencies and therefore are quite small relative to the proportional 
term plus bias. We conclude then that all "nonlinear" terms can be approximated by 
functions of the instrument outputs. The following discussion shows how this is accom- 
plished. 

The accelerometer model (neglecting quantization and noise) has the form: 

An = D^J (a • A) dt + DjDqJ dt + higher order terms 

where An is the difference in the number of zero crossings detected from each string m 
the tune period over which the integrations in the equation are made. In the following 
discussions, this time period will always equal At. 

As a first order approximation, 

. 1 
j (a • A) dt = (An - D. D At) 

D 1 

Let us define 

(a* A) At s J(a.A)dt 
J (a • A) dt 1 

or (a. A) = — =-= = (AN - D<D At) 

J dt D x At 1 ° 


We can see that (a • A) is, from the mean value theorem of calculus, a value of (a • A) 
somewhere m the time period of integration. 

Referring to the second order term 


J(a- A) 2 dt 
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which can be written 


(a - A) 2 At 

We see that, when the period of integration is chosen short enough that (a • A) is essentially 
constant over each period, then it can be assumed that. 

(a -A) 2 = (aT* A) 2 


(This is a very good approximation, considering the fact that the coefficient cf the square 
term is quite small.) From all of the above statements we can therefore infer: 


J (a- A) 2 dt 


D 1 At 


(An - Dj_D At) 


At 


Similarly, the cubic term can be written as 


| (a - A) 3 dt 


' 1 

( AN - D jD At) 

Dj At 1 0 


3 


At 


The expression for (a* A) can also be used to determine the gyro unbalance integrals from 
the approximations: 


where 


J(a.G k )dt 


J(a- A k )dt 



0 0 1 


J (a. A x )dt 

I(a-P k ) dt = 

10 0 


J(a- A 2 )dt 


10 0 
L. — 


J(a- A 3 )dt_ 


0-10 


J (a - A|) dt 

J(a-S k )dt = 

0 0-1 


J (a* Ag) dt 


0 10 


J (a • Ag) dt 


(See Figure 2-3.) 


J(a. A k )dt 



(We have refrained from using the instrument index k until it was absolutely essential. 
This served to keep the notation as simple as possible. ) 
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The compliance integrals are found to be similar to the higher-order accelerometer terms. 
We note that all of the gyro approximations utilize the outputs of the accelerometers to 
compensate for the gyro acceleration-sensitive terms. 

After all approximations, we will have equations of the form: 

f (a • A^) dt = f (accelerometer outputs) 

J (u;* G^) dt = g (accelerometer and gyro outputs) 

A G 

The Q and Q matrices can then be used to find. 

J(a.B k )dt = 2 qA J(a-A^)dt 

Kj 

jVB k )dt = SqG J(co.GJdt 

The above statements lead to the complete set of computations, which are found on the 
following chart. We will henceforth, in this document, refer to those equations as the 
Preprocessing Computations. The following nomenclature is required for the understand- 
ing of the Preprocessmg Computations: 

■Hi 

• (Sy^i and (2y)]j2 are the counts from the one and two strings of the k Lil 
accelerometer. 

in 

• ( s 6) k is the count from the k gyro. 

T 

» S 2 is the clock scale factor (the subscript 2 serves to distinguish this scale 
factor from another used in calibration). 

T 

® (Sn 2 ) is the count from the system clock. All other terms have been 
previously defined. 

2.2.6 Calibration Requirements 

The Preprocessmg Computations developed in the preceding subsection are seen to be a 
function of a great number of constants. Those constants were, in Section 2. 1, defined 
as the constants to be calibrated. An explicit statement of the calibration requirements 
is therefore the determination of the quantitative value of the constants contained within, 
the Preprocessing routme . In Chart 2-5 those constants are lusted. As a matter of 
convenience, the nominal values, ranges, and precision requirements, where available, 
are given. 
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CHART 2-4 


PREPROCESSING COMPUTATIONS 
Inputs (Sy) k2 ,(Sy) kl , (U6) k ,and (unj )for k = 1, 2, 3 

The outputs B k )clt and J^+ At (a-BjJdt (k = 1, 2, 3) are given by the following 

computations 


* p k £ [ Csr) k2 - ter> kl ] 

* *> k E(£6) k 

. At -s|(En|) 

* L(«.G k )At: = P k (A$) k - (R) k At 

« i%A k )At; - ^/(D l)k - (D 0 ) r At 

* (^5 k ) 5 C(u;-G k ) At j/ At 

* (a • Aj.) £ [ (a • A k ) At]/ At 


* (a • G r ) s (a- A fc ) 



~o o r 

Ha- A x ) 

<*5? - 

10 0 

<r^ ) 


1 0 0 

L€^). 


~0 -1 o' 

’(a* Aj) 

<*I? = 

0 0-1 



0 10 

(a- A 3 ) 


• „ t t " At (ii-Gk)d t = l (<£-G k )At j-*. (B I ) k (a.G k ) + (B 0 ) k (a.O k ) + (B s ) k (a»S k ) .At 
- u( c n ) k ( a. G k ) 2 t- (C ss ) k (a-S k ) 2 : At 

“ '-(Gig) k (a i Gj.)(a. S^.)- (pQc;) k (a- O k )(a S k ) + (CjQ) k (a-G k )(a-O k ) .At 
" *- (Qjj) k (uJ‘G k ) -r (Qxs)rC — - k )( -‘-k^ -^t 
V' t t_At (a-A k )dt = I(a.A k )At: - (D 2 ) k (a 7 A^) 2 At - O g ) k (a- A k ) 3 At 

.j^ t+A V ± -B k )clt = EQ^ Lt t ’ At (w.G J )dt 

•J^V^dt = EQ* .^t-^Ca.A^Jdt 
where 


1 

-(G X .B 2 ) 

-(Gj-Bj) 

-(g 2 -Bi) 

1 



-(G 3 .B 2 ) 

1 

1 

”(—].* Sa) 


-(A 2 -§!) 

i 

-(a 2 * § 3 ) 


-(A 3 -B 2 ) 

1 
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CHART 2-5 


LISTING OF CALIBRATION PARAMETERS 
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Note that a great number of the constants have a nominal value of zero. Note also that 
the accelerometer bias term is the only term which has a nominal value which appreciably 
affects an instrument output (that is, affects the output over and above the effect of the 
linear term). The relatively large value of the bias term accounts for the bias being 
utilized in the preceding subsection as a part of the higher-order term approximations. 

i 

Only the precisions of the bias and scale factor are given. The scale factor precision is 
presented as a relative requirement (ratio of uncertainty to magnitude). The bias 
precision is presented as an absolute error. It will be shown in the trade-off document 
that the errors in each of the other terms act as either a scale -factor -like or bias-like 
error. Therefore, all precisions are inferred from either the scale factor or bias 
precision. 

2. 3 ALIGNMENT REQUIREMENTS 

In Section 2. 1 alignment was defined as the initialization of the matrix which transforms 
from an ISU-fixed frame of reference (body axes) to a navigation frame. In this section 
our purpose is to explicitly state the requirements for determining that matrix. Three 
alternative techniques will be presented. The definition of the three alignment techniques 
will be presented in Section 2.3. 1. The alignment requirement associated with each 
technique will be found to be the measurement of either the body -axes or earth-axes 
components of two system vectors. In Section 2, 3.2 the geometry associated with the 
alignment techniques will be presented. The explicit functional form of the alignment 
matrix for all three techniques will be delineated in Section 2.3.3. 

2.3.1 Definition of Three Alignment Techniques 

Alignment has been defmed as the initialization of the body-to-navigation-axes transfor- 
mation matrix. For the subject ISU, the body axes (B fc ) are defmed by the normals to 
two ISU-fixed mirrors (see Section 2,2. 1). For the purpose of this study, we will assume 
that the navigation axes are aligned with a set of local-level earth axes (E k ), where 

E i is directed up (along the line of local gravity) 

Eg is directed east (normal to the local meridian) 

Eg is directed north (normal to E^ and Eg) . 
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(Throughout the text we will, at times, also refer to the earth axes as U, E, and N, where 

U = E 1 
E = E 2 
N = E 3 

We can now more explicitly define the alignment problem as the determination (at some 
time, say t) of the 3x3 orthogonal matrix T, where T is defined by: 


1 

ltd 

jifc 

1 


"u* 




V 

a 

= 

E 

= 

T 


—2 

A 


N 




1 

CO 

ffll 

r 


_ _ 


- 



There are numerous techniques for determining the elements of the matrix T„ Each 
technique considered in this report is based upon expressing the matrix functionally in 
terras of the components of two vectors (which are known, in an operational situation) 
in both the body- and earth-fixed frames. The typical operational situation would be 
an earth-fixed orientation of the ISU. Assuming that the operationally available meas- 
uring devices to be used during alignment are any combination of three gyros, three 
accelerometers, or two two-degree -of -freedom autocollimators, then the vectors which 
can be used to functionally define the T matrix are: 


• The unit mirror normals and M 3 , which can be measured in the earth frame 
by the autocollimators and are known in the body frame because they define the 
body frame 

• The local environment vectors g and w E ("gravity"* and earth rate), which can 
be determined in the body frame by the accelerometers and gyros and which are 
known in the earth frame because they explicitly define the earth frame, that is, 



where U = 


and W = 


U 

(WxU)/| W X U| 
Ux(WxU)/| Wx U 


g/ }S 


.E 




The T matrix can be expressed in terms of the components of any two of the four above- 

O' 

mentioned vectors (M 1} M 2 , g, and w ); but, as a contract requirement, only the 
following three combinations are of interest in this study: 

*See Section 3. 1. 1 for a definition of "gravity". 
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• and Mg 

• M-^ and g 

• g and w . 

The names of the techniques which implement the use of the components of these three 
vector combinations are, respectively: 

• Mirror Alignment 

• Level Alignment 

• Gyrocompass* 

In the next subsection we will present the geometry of the four vectors (M ^ Mg, g, 
and w® ), and in the following subsections we will present the explicit mathematical 
relationships between the components ctf those vectors and the elements of the T matrix. 

2.3.2 Alignment Geometry 

In the preceding subsection we described alignment in terms of the determination of the 
components of two vectors in both the body- and earth-fixed frames. We chose, as 
alternatives, the vector combinations 

• M i and Mg 

• and g 

• g and w . 

In this subsection we present the geometric relationship between the four vectors which 
are considered in our three techniques. 

The required geometry is shown in Figure 2-7. The following comments explain the 
notation: 

• U is a unit vector directed up; that is U = g /j g| • 

• W is a unit vector directed along earth-rate; that is W = w E /j « E j. 

• a i andQ!2 are, respectively, "azimuths" of the one and two mirror 
normal, as determined by an autocollimator. 

• 01 and 02 are, respectively, the "zeniths" of the one and two mirror 
normal, as determined by an autocollimator. 
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Figure 2-7. Earth and Mirror Axes 
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• X is the angle between U and W, which we will refer to as "local latitude". 

• Pis a unit vector directed along the projection of the one mirror normal 
into the local horizontal plane. 

2.3.3 The Three Functional Forms of the Alignment Matrix 

In this subsection we present the explicit mathematical relationships between the elements 
of the T matrix and the measurable body or earth referenced components of 

• M-j and M 2 for Mirror Alignment 

• Mi and g for Level Alignment 

• g and t*) E for Gyrocompass 

These functional relationships are presented on the three accompanying charts. We will 
not, in this subsection, meticulously derive the relationships; but will instead present 
sufficient information such that the derivations are obvious. Let us consider each chart 
in turn. 


Mirror Alignment 

This derivation is quite easily explained. In Section 2. 2. 1 we defined the body-axes as: 

§1 = Mi 

(M^ x Mglxh^ 




®3 




xM, 
(Mj x M 2 ) 


MjxMf, 

Therefore any vector, say U, can be written m body coordinates as 


U = (U- + 


U * (M 1 xM 2 ) xMj U • (M 1 X M 2 ) 


x Mg 


*2 + 


M^xMg 


h 


and, after substituting the identity 


U = ExN into the last component 
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CHART ^2-6 


MIRROR ALIGNMENT MATRIX 


| Inputs 8^, Oj, Sg and 

From these quantities the alignment matrix is given by: 


(U-Mi) 

(MjXUHMjXMg) 

(ExN)-(M^x Mg) 

iMj xm 2 I 

|M : x Mg 1 

(E-M^ 

(M 1 xE)-(M 1 xM 2 ) 

(NxO).(MjXM 2 ) 

iMi xm 2 I 

| M i x M 2 ! 

(N.M^ 

(Mj xN)-(MjXM^) 

(U x E)*(Mj x Mg) 

iMi xm 2 I 

|MiXM 2 I 


where 

IMjX M 2 | = [1 - (M x * M 2 ) 2 ] 1//2 

(M x • M 2 ) = (M^ U)(M 2 . U) + (M^ E){M 2 - E) + (fc^. N)(M 2 - N) 




cos 9 ^ 


(u-m 2 ) 


cos 8 2 

(E-Mj) 

= 

cos sm^ 


(e-m 2 ) 

= 

cos o^2 sin ® 2 

(N-Mj) 


sin j sin 9 ^ 

y 

(n-m 2 ) 


sin « 2 sin 6 2 


An optional technique might utilize the value of lIV^ x Mgl from a previous alignment 
thus eliminating the aforementioned dot product and square root operations. 





and interchanging the dot and crass in the second component, we have 


U = [(U-M^Bj 


(M^UWMixMj,) 


ExN)-(M lX M 2 ) 

| MjxMgj 

§2 + 

l^xMgl 


*3 


The three bracketed quantities are the elements of the first row of the T matrix. In the 
Mirror Alignment chart we see those elements in the first row, and similar elements for 
E and N in the second and third row. The relationships between those elements and the 
azimuth and zeniths as determined by optical equipment are listed below the matrix. The 
azimuth and zenith relationships are obvious from Figure 2-7. 

Level Alignment 


This derivation is quite simple if one separates the problem into three parts by defining 
matrices (Tp T 2 , and Tg say), where the three matrices are defined by 



U 

P xu 
P 




U 

T 1 


P x U 



P 




"Ml 

T 2 


U 

M-. x U 

_ _ 


- i _ 


Mi 




*1 

u 

= 

T 3 


— 2 

M-£ xU 




L ^3 


Obviously T = Tj T 2 Tg. 

The three matrices shown in the Level Alignment chart are Tp Tg, and Tg respectively. 
The derivation of the T ^ matrix is obvious from the definition of P in Figure 2-7. The 
derivation of T 2 is based upon the fact that 

PxU=(MiX U)/J MjXU 
and P= Ux(M 1 xU)/|m 1 x U | 
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CHART 2-7 


LEVEL ALIGNMENT MATRIX 


Inputs (g- B t ). (g . B 2 ), (g * B s ) and a. 


From these quantities the alignment matrix is given by: 




10 0 


0 1 0 


1 0 0 

T 

= 

0 sina^ cos a . 


1 

0 0 iM.xUl 


(U-B^ (y*B2) (U-Bg) 





1 (Mi-U) - 1 - 





0 -coso^ smttj 


|m : xuI iMjxul 0 


0 -(U-E^) (U-Bg) 


where 

• (Mi* u) = (u- b x ) 

• IMjxUN Cl - (M^U 2 ] 1 / 2 

0 (U* Bjj) = (£' 

<• g = [ (g- B 1 ) 2 + (g- B 2 ) 2 + (g • B 3 ) 2 ] 1/2 

An optional technique might utilize any of the following additional inputs: 

• The zemtlTahgle (S^) of mirror one might be utilized to find (M^ • U) from 

(Mj- U) = cos 0 1 

« The magmtude of gravity (g) might be supplied from a local survey. This piece 
of information can be utilized to reduce the number of required accelerometers 
to two. 
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The derivation of Tg is based upon the fact that 


Ml = b x 


U = (U* Bj) B x - (U- B 2 ) B 2 4 (u- Bg) Bg 
MiXU = ( M-^xU* Bj) + (MjXU- § 2 ) Bg + (Mj x U- Bg) Bg 


The third identity becomes, after substituting 


M^U = - (U* Bg) Bg + (U- Bg) Bg 

Belov/ the matrix expression we see the obvious relationships between the elements of the 
Tg and Tg matrices and the body components of g(as determined by accelerometers). 
Below those relationships, we see alternate methods that utilize the zenith cosines and 
sines from an optical measurement, and g from a survey. 

Gyrocompass 


In this derivation we express the T matrix as a product of two matrices (T^ and Tg say) 
where T^ and Tg are defined by: 


5.1 




W 

5 2 


T 4 


u 

_ ^3. 


- - 


Wx U 

w" 


- 


s_r 

u 

= 

T 5 


*2 

WxU 




*3 


That separation is shown on the gyrocompass chart. The matrix T^ is obtained from the 
identities 



CHART 2-8 


GYROCOMPASS MATRIX 
Inputs (go (g. B 2 ) ; (g . B 3 ) } (w E . Bj) (w E - B^ and (o. E . Bg) 

Prom these quantities the alignment matrix is given by: 

0 1 0 (W-B^) (W* Bg) (W-Bg) 

1 

0 0 IWxUl (U'B 2 ) (U-Bg) 

1 (w-u) 

I W x U I " Iwxu] 0 (W X U)> (§ 2 x Bg) (Wx U) • (BgX B x ) (Wx U) • (B ^Bg) 

where 

• (WoU) = (W.B^CU.B^^ (W-B 2 )(y«B 2 )+ (W»Bg)(U.Bg) 

• I Wx uj = Cl - (W-U ) 2 ] 1 / 2 

• (W.BjP - (w E . B k )/w E 

• (U’B^ -(S'^/g 

• W E =C (« E * B t ) 2 + (w E - B 2 ) 2 + (co E . Bg) 2 I 1/2 

• g = [ (g • B 1 ) 2 + (g . B 2 ) 2 + (g • Bg) 2 ] l/2 

An optional technique might utilize any of the following additional inputs:- 

• The local latitude (X) might be utilized to find (W • U) from 

(W • U) = cos X 

0 The magnitude of gravity (g) might be supplied from a local survey. 

8 The magnitude of earth rate (w®) might be supplied from a local survey. 

A use of all additional inputs could reduce the number of necessary instruments to three 
(either two accelerometers and one gyro, or one accelerometer and two gyros). 
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The rows of the matrix are obviously made up of the components of W, U, and WxU. 

Below the matrix expression is found the relationships between the elements of the and 
T,. matrices, and the accelerometer and gyro determined body-axes components of g 

0 t? 

and oj . Below those relationships is found a discussion of alternate techniques utilizing 
g, W E and x. 
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SECTION 3 

SYSTEM DESCRIPTION 


The calibration and alignment requirements for the NASA/ERC strapdown inertial 
sensing unit were presented in Section 2. The calibration requirements were defined 
as the determination of the inertial instrument model constants and the elements of 
the matrices which transform between the instrument and body frames of reference. 
Alignment was defined as the real-time initialization of the body-to-earth transformation 
matrix. The alignment requirements were defined as the measurement of the body 
and/or earth-frame components of two system vectors. Three alternative choices of 
sets of vectors were introduced. These alternatives characterized the three alignment 
techniques: Mirror Alignment, Level Alignment, and Gyrocompass. As a necessary 
aid to the satisfaction of the calibration and alignment requirements, various pieces of 
laboratory equipment are needed. It is also necessary to understand fully the nature of 
the kinematic environment in which the equipment and ISU are located. In this section 
our purpose is to describe that equipment and environment, beginning with their relation- 
ships with the problems of calibration and alignment. 

3. 1 FUNCTIONAL DESCRIPTIONS OF CALIBRATION AND ALIGNMENT 

The following paragraphs describe the functional activities of calibration and alignment. 
(Section 3. 1. 1 discusses calibration, and 3. 1.2 discusses alignment.) These functional 
descriptions serve as a definition of the required inputs to the calibration and alignment 
evaluations which come from sources other than the ISU. 

3. 1. 1 Calibration 

In the following discussion we indicate ihe functional requirements for determining the 
calibration numbers. Our discussion will be quite general, the major purpose being to 
introduce the reason why the equipment described m subsequent subsections is required. 

As an aid to our presentation, we find it useful to compare the ERC laboratory calibration 
with an "Ideal Calibration". This comparison serves to indicate the compromises which 
are necessary in defining an operational calibration laboratory. 

On the following chart (entitled Calibration) we present two calibration functional diagrams. 
We refer to the two techniques represented by those diagrams as Ideal Calibration and ERC 
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CHART 3-1 


CALIBRATION 


IDEAL CALIBRATION (OR STANDARDIZATION) 


Real-time 
environ- 
ment inputs 
into both 
the ISU 
and the 
standard 


instrument Mapping 


ISU outputs as a 
function of standard 
outputs over the 


entire expected 
operational range 
of the ISU 


Standard 


The integral of 
acceleration and 
angular rate 


ERC LABORATORY CALIBI 


Assumed 
Math Models 


Real-Time 

Inputs 



INSTRUMENT 

Computation 


of 

Constants 


Instrument Constants 
and Q Matrices 



Measurement 
Devices Used 
In Real Time 


il of the 
er Rate 


Non Real-Time 
Inputs 



Store 


Average g, «j and 
deviations from those 
averages as measured 
sometime prior to 
the real-time 
calibration. 


3-2 







Laboratory Calibration. The Ideal Calibration diagram represents the manner in which 
calibration would be accomplished if an unlimited amount of time and money were avail- 
able. The ERC Laboratory Calibration diagram represents the maimer in which the 
calibration will be accomplished under more realistic constraints. 

If one wished to calibrate a "black-box" (ISU) ideally, he would operate in the manner 
indicated in the Ideal Calibration functional diagram. Subsequent to the development 
of an ISU, whose outputs are designed to vary over a range of kinematic environment 
inputs (a and to), one would wish to determine the quantitative relationship between the 
ISU outputs and that environment. However, that environment exists only conceptually, 
and not quantitatively, until a device is available which is defined as the measurer of these 
kinematic quantities. That device is referred to as the "Standard" measuring device. 

With the Standard available, it is then possible to calibrate the ISU by placing both the 
ISU and the Standard in the same environment and mapping the output of both over the 
range of the kinematic quantities which are considered to be significant. This mapping 
would take the form of a table of ISU and Standard outputs over the required operational 
range of the ISU. The mapping would necessarily be accomplished in a frame which 
characterizes the ISU and is known relative to the Standard, In our case, this frame 
will be the body axes as defined by the ISU mirror normals. 

We see from the aforementioned statements that calibration is nothing more than the 
implementing of the requirement that the ISU behave as the Standard would under the 
same kinematic conditions. Thus, after calibration the ISU will have been "standardized". 
Subsequent to the standardization, it is assumed that the ISU can measure the kinematic 
environment, as the Standard would under the same conditions. This is accomplished by 
a transformation of the ISU outputs into a measure of the environment by use of the mapping 
information. 

This is all rather interesting but not, operationally, very feasible. First, such a Standard 
is not available in the laboratory, and even if it were, time would not allow for a mapping 
over the entire operational range of the kinematic inputs. Secondly, m the case of applied 
acceleration, the typical operational range of the kinematic inputs cannot be easily 
generated m the laboratory. (A centrifuge would .be required for accelerations higher 
than one g.) Thirdly, it is not always feasible to have even a substitute for the Standard 
operating at the same time as the ISU. All of these problems explain the deviations of 
the ERC Laboratory Calibration from the Ideal Calibration. Before elaborating those 
differences, a description of the ERC environment is necessary. 
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The kinematic inputs found in the ERC laboratory include: 


1) The applied accelerations and angular velocities characteristic of any point 
on the earth's surface 

2) The local deviations from those accelerations and angular velocities due to such 
things as earthquakes and cultural noise 

3) The generated environments, caused by the ERC test table. 

The first category includes earth-rate and the applied acceleration (normal- specific 
force) which negates the acceleration due to gravity in a "static" orientation relative 
to the earth. (This acceleration is often confused with gravity. It is, on the average, 
equal in magnitude and opposite in sign to gravity. The very common convention is to 
refer to this applied acceleration as g. We will, in the remainder of this document, also 
refer to it as g. Note, however, that we always direct g away from the surface of the 
earth.) The second category will be referred to as "noise". The laboratory test table 
(see Figure 3-1), mentioned m the third category, has a motor-driven capability of 
rotating at speeds up to several thousand earth rates. Such rotations will develop angular 
velocity and angular -velocity-related acceleration inputs to the ISU (the ISU being always 
attached, during calibration, to the table). 

As suggested in the discussion of the Ideal Calibration, it is necessary that an independent 
measure (Standard Output) of the laboratory environment be available in order that 
calibration can be accomplished. This independent measure, even when a substitute for 
the Standard is used, should be accomplished at the same time that the ISU is yielding 
outputs. In the case of the angular velocity of the test table, an independent, real-time 
measure will be accomplished through the use of the output of the test table resolver. 

Tp 

The measurements of the g and u vectors are, however, accomplished at some time 
prior to calibration and "stored" for use during calibration. The storage of the direction 
of those two vectors is evidenced in the location of such things as optical lines, resolver 
zeros, etc. , and the magnitudes by storage of numbers in a computer memory. Informa- 
tion about the noise is stored in the form of graphs showing characteristics such as 
power spectral densities. Because g, to , and noise are not measured in real time, it 
is assumed that their behavior is the same at the time of calibration as it was at the 
time of measurement, therefore, they can be considered a good approximation of a real- 
time measurement. 

To this point we have described the manner in which the Standard output is evidenced in 
the ERC laboratory. We require only one more statement, in this presentation of the 
calibration functional activity, about the substitution of ERC Laboratory Calibration for 
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Ideal Calibration. The Ideal Calibration was described as a standardization over the 
entire range of inputs. Fortunately, m the case of Laboratory Calibration, some of the 
standardization has been accomplished by instrument designers prior to the placement 
of the accelerometers and gyros in the ISU. A great deal of time and effort has already 
been devoted to the development of a functional relationship between the output of the 
instruments and their environment inputs. Those functional relationships are referred 
to as mathematical models. As seen earlier, the instrument models contain many con- 
stants. Because of the availability of the models, it is only necessary in calibration to 
map a number of environments equal to the number of calibration constants for the de- 
termination of those constants. It is assumed that a knowledge of the models, and the 
model constants, serves to interpolate the mapping between the chosen environments. 

In this subsection we discussed the independent measurements required as an aid to 
calibration. In Section 3.2 we will discuss how those measurements are specifically 
developed as inputs into the determination of the calibration constants. 

3. 1.2 Alignment 

In Section 2. 1 alignment was defined as the initialization of the matrix which transforms 
from an ISU-fixed set of axes to a navigation set of axes. In Section 2.2 the ISU-fixed 
axes were defined by two ISU-fixed mirror normals, and in Section 2. 3 the navigation 
axes were defined as an earth-fixed, local-level frame of reference. Further, in 
Section 2.3, three alternate mathematical forms of the alignment matrix (T) were derived. 
Each form showed a requirement for a different set of optical or inertial-mstrument 
measurements as an input into the quantitative determination of the alignment matrix. 

In this section we will discuss the techniques for determining each set of inputs. 

As an aid to this disuession we present in the accompanying chart, entitled "Alignment 
Functional Diagrams", a schematic of each of the three alignment techniques. In the 
remaining paragraphs of this section we discuss, m turn, the contents of each functional 
diagram . 

Mirror Alignment 

The routine labeled Alignment Matrix Computations represents the computations described 
in the Mirror Alignment Chart found m Section 2. 3. 3. As shown in Section 2.3. 3 those 
computations require, as inputs, the optically determined azimuth and zenith of both the 
one and two mirror normals. (In practice, the actual optical measurements might be 
angles other than the azimuth and/or zenith angle. It is always an easy matter, however, 
to convert the actual measurements into the required azimuth and zenith. ) The optical 
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measurements will be manually transferred to the matrix computation routine, which will be 
part of a digital computer program m the ERC facility. Because optical measurements are 
extremely accurate by design, we do not include in the functional diagram any data-filtering 
function. 

Level Alignment 

The second functional diagram represents the Level Alignment procedure. The Alignment 
Matrix Computations indicated m the last block in the diagram represent the computations 
found m the Level Alignment chart described in Section 2. 3. 3. That routine requires, as 
inputs, the body-axes components of g and the azimuth of the one mirror. In Section 2. 3 it 
was mentioned that the body-axes components of g will be available as the result of acceler- 
ometer measurements and the azimuth of the one mirror as an optical measurement. At the 
left side of the diagram we see the input of these measurements. The optically obtained 
azimuth goes directly to the matrix routine {as it did m Mirror Alignment). The accelerom- 
eter inputs, however, will require further processing, since they will be in the form of 
three digital pulse counts. 

We saw, in Section 2, 2.5. that the Preprocessing Computations convert such counts into 
integrals of body-axes components of the applied acceleration inputs to the accelerometers. 
However, those computations assumed no quantization and instrument noise. Therefore, 
the transformation of the outputs of the Preprocessing Computations into the desired 
body- axes components of g would require four additional operations in order to accomplish 
a good estimation of g- . These are: 

• A differentiation of the integral outputs of the Preprocessing Computations 

• A compensation for instrument quantization 

• A compensation for instrument noise 

» Separation of g from random environmental accelerations. 

If the ISU were to be in a stationary orientation relative to g during alignment {that is, if the 
accelerometer input were a constant g acceleration), the first operation would be simply a 
division of the Preprocessing output by the total time of integration (say At). Additionally, 
the compensation for instrument quantization could then be accomplished by simply waiting 
sufficiently long such that the quantization residual would be arbitrarily small. However, 
the ISU at the ERC facility will not be m a constant g ennronment, Due to such things as 
local vehicle motion, personnel movement, etc. . the ISU will be. in fact, in the presence of 
the nominal local g plus "noisy” vibrations. If some a priori knowledge of that noisy envi- 
ronment is available, it is possible to accomplish some of the aforementioned compensations 
by the development of mathematical filtering operations on the "Preprocessing” outputs. 
Those mathematical operations are presented in Chart 3-2 as the Estimation Matrix Com- 
putations and Estimation Routine. The former involves the computation of constants prior 
to the actual alignment. The input to those computations is the a prion noise information. 
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The Estimation Routine represents the on-line operations on the Preprocessing outputs. The 
outputs of the Estimation Routine are the required estimates of the body-axes components of g. 

The Preprocessing routine has been completely defined m Section 2.2.5, ana the Alignment 
Matrix Computations have been defined in Section 2. 3.3; hence, the development of the 
alignment techniques presented in Section 5 of this report will be preoccupied with the 
Estnnation Matrix Computations and the Estimation Routine. In Section 3.2.2 we will 
present a discussion of the a priori noise information which is required as inputs to the 
estnnation routines described m Section 5. 

Gyrocompass 


The third functional diagram in Chart 3-2 represents the operational gyrocompass pro- 
cedures. The Alignment Matrix Computations shown m the diagram were presented in 
Section 2.3.3. Required inputs are the body-axes components of g and At the 

left-hand side of the diagram we see the inputs of accelerometer and gyro readouts 
required for the determination of the body -axes components of g and w®. in the preceding 
section there is a discussion relating to the transformation of the accelerometer outputs 
into an estimate of the body-axes components of g. This discussion also applies to Gyro- 
compass with the following modifications 

* In Gyrocompass the Preprocessing computations as presented m Section 2.2.5 
will be used entirely, whereas the Level Alignment uses only the accelerometer- 
related computations. 

° The estimation routines will operate on both gyro and accelerometer dala. 

3.2 ENVIRONMENT MODEL 

In the preceding discussions we showed that independent environment measurements are 
required for calibration, and a priori noise information is required for alignment. 

We indicated the manner in which that information is available at the ERC facility. In this 
section we will show specifically how the required measurement information is made 
quantitatively available to the calibration and alignment computational routines. 

The reference environment information takes on different forms, and therefore can be 
discussed independently. First, there are the stored g and u-' E vectors which must be 
expressed in terms of body-axes components for calibration purposes, There is the -c 
vector (angular velocity of the test table), which must also be expressed in body -axes 
components. We will refer to these three vectors as the deterministic environment. We 
will show in Section 3.2, 1 how the body-axes components of the deterministic environment 
are obtained as a function of test table gimbal angles. The remaining environment inputs 
have been referred to as random noise. They will be described in Section 3.2.2. 


3-9 



3.2.1 Deterministic Environment 


E T 

The deterministic environment is made up of the three vectors g, ^ , and <£ , which 
are assumed known in frames well surveyed in the laboratory. Our purpose in this 
section is to show how those vectors are transformed from the laboratory frames into 
ISU oody-axes components. The transformation will be accomplished through the use of 
quantitative measurements taken from both the laboratory test table and the system 
autocollnnators. We begin our discussion (in 3.2. 1. 1) by describing the geometry of those 
pieces of equipment. We then define the transformations between the many rigid bodies 
making up the equipment. Finally, (m Section 3.2. 1. 3) we will develop the operational 
transformation of the deterministic environment into body -axes components. 

3.2.1. 1 Laboratory Geometry 

The geometry of the test table and autocollnnators is the geometry which enables us to 
E T 

transform g, to , and oj into body-axes components in the ERC laboratory. In 
Figure 3-2 we present a schematic of this geometry. This figure is a repeat of Figure 3-1, 
with the addition of the defined laboratory frames. Chart 3-3 presents the definitions of 
the frames indicated in Figure 3-2. A few comments are necessary as an aid to the 
understanding of Figure 3-2 and the chart containing the definitions of the frames. 

• All frames are defined by orthogonal unit vectors directed along the frame axes. 

• The Sk frame is not explicitly defined in Figure 3-2. The explicit definition 
will depend upon the (at this time) unknown geometry of the autocollnnators. 

(The lack of an explicit definition has, however, proved to be no burden in the 
work that follows. ) 

• The body axes are not shown in Figure 3-2, because their relative orientation 
depends upon the manner in which the ISU is attached to the inner-axis rigid 
body. 

° The Fk frame will be required to line up with the !k frame. This alignment 
will naturally be with the laboratory frame, which, in turn, is thought to be 
coincident with the earth axes defined by g and wE, We will see in Section 4.4.4 
that this alignment will have to be corrected periodically by the use of bubble 
levels. 

° Each adjacent pair of test table frames is assumed to have a common axis. 

• The four test table rotation angles are defined as cu, and Cb, as 

shown in Figure 3-2. 

• The test table orientation shown in Figure 3-2 is the zero orientation — that is, 
the orientation when all resolvers yield a zero output. 

° The ISU will be attached to the section labeled "test subject" throughout the 
entire calibration. 
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Figure 3-2. Laboratory Frames 
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A triad of unit vectors fixed to the body con- 
taining the rotary axis 
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taining the inner axis 

• ?i ?2 -3 


A triad of unit vectors defining the body axes 
as defined by the mirror normals (see 
Section 2. 2. 1) 
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3.2. 1.2 Definition of Transformations 

Let us choose the and frames for our example. The relationship between the unit 
vectors of the two frames is given by: 



th RT 

where a multiplication of the k row of the 3x3 matrix T with the T column represents 

— 1 L 

expressed in the T,__ frame. If x is any vector known m the frame and one wishes 
to express that vector in the R^ frame, we dot the above definition with x, yielding: 



te-Tjf 

(5-T 2 ) 

T 3 ) 


If, further, we wish to transform to the frame, we have: 


(x-Fj) 




and so forth. 


3.2. 1.3 Operational Transformations 

With the geometry information now completely described, it is possible to show how 
E T 

g, « , and w are transformed into body-axes components. The accompanying chart 

shows how that transformation is accomplished. Note that the chart specifies that the 

test table gimbal angles will always be used in determining the transformation, instead 

of autocollimator surveys. This is purely a matter of convenience. It certainly would 

be cumbersome to survey via the autocollimators over the 4F steradians in which the 

mirror normals can be located: Besides, the test table was designed to accomplish 

the necessary transformations. Note that the autocollimators are absolutely essential 

RT 

for one very important operation, namely the determination of the matrix T . 
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During calibration there is a requirement for the transformation of w , g and a> into 
body axes components. Gravity and earth rate are vectors explicitly known m the 
frame and u) F is a vector explicitly measured in the R^. frame. The body axes com- 
ponents of these vectors can, therefore, be written: 


(w E . Bfc) = V T?? a- E . 


E. 


(s T %)-? , "w “ t -5, 


BE RD 

The* matrices T and T are, therefore, required. These matrices, as a matter of 
convenience, will always be found as a function of test table gimbal angles. Therefore 

rpBR _ ,j,BI rpIO rpOR 

and t EE — x R ^ x^^ q'OR rpUT x^^ 

BI 

The matrix T is a constant which must be found from an initial survey (see Section 
4. 4. 1). The remaining matrices can be seen from the previous definitions to be: 
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o“ 


“~0 

1 

0“ 

lT RT ] = 

C0 2 

0 

-S0 2 

ii 

■ — i 

H 

E-< 

C0J 

0 

-S0 i 


-S0 2 

0 

-c 


-S0 4 

0 

-C0 L 


LT FE ] = 


10 0 
0 10 
0 0 1 


(Note, the frames are defined so that each matrix has the same form. ) 
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3.2.2 Random Environment 


The laboratory test environment introduces a random translational and rotational "noise" 
input. In describing the noise inputs, we will utilize the worst-case model. Namely, the 
translational input and angular input are assumed to be independent, and the components 
of the random input vectors are assumed to be independent. 

The random translational inputs about the up-down, east-west, and north-south axes 
have been assumed to be statistically independent and identically distributed. Each input 
has a power spectrum illustrated in Figure 3-3.* It follows that the translational motions 
along any three perpendicular axes are uncorrelated and have the spectrum given in 
Figure 3-3. The nominal input (local gravity) is assumed to be the long-term average 
input. Assuming ergodicity of the expectation, the translational inputs have zero expec- 
tations. The random translational inputs do not produce a significant output from the 
gyros. 

The random rotational input produces a rotation about an axis in the horizontal plane. 

The rotation about the vertical axis can be neglected. The random angular inputs about 
the east -west axis and the north-south axis are statistically independent, and each has 
the power spectrum illustrated in Figure 3-4. These inputs are assumed to have zero 
expectations. 

The spectra given m Figures 3-3 and 3-4 are the basis of numerical computations involving 
environment noise. The selection of the "Recommended Alignment Techniques" in 
Section 5. 7 assumed an environment as indicated in Figures 3-3 and 3-4. The alignment 
processing techniques derived in Sections 5.4 and 5. 5 use the power spectra of the 
translational and rotational noise inputs but do not depend on the specific numbers given 
m Figures 3-3 and 3-4. 

3. 3 HARDWARE DESCRIPTION AND INTERFACE 

The material presented to this point has been introductory in nature. That is, all dis- 
cussions were either related to the detailed statement of the calibration and alignment 
problem, the definition of terms, the description of necessary equipment, or the descrip- 
tion of the laboratory environment. In this section we will complete the presentation of 
introductory material by covering two descriptive tasks which aid in the understanding 
of the calibration and alignment development in Sections 4 and 5 and the operational 


* Spectra data is given by H. Weinstock in "Limitations on Inertial Sensor Testing 
Produced by Test Platform Vibration", NASA Electronics Research Center, 
Cambridge, NASA TN D -3683, 1966. 
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procedures presented m the Laboratory Procedures Manual. First, there is laboratory 
equipment which is necessary for the implementation of calibration and alignment pro- 
cedures and which has not yet been described. Second, the laboratory equipment has not 
been presented as a system as a whole, with all of its interconnections. In this section 
the description of that equipment and interface will be accomplished through the use of 
a system diagram. In Section 3. 3. 1 that diagram is introduced but without the lines 
denoting the equipment interface. In conjunction with the introduction of that diagram, 
we will present brief descriptions of the equipment which each block in the diagram 
represents. In Section 3. 3.2 the useful interfaces between all equipment will be presented 
in tabular form. 

3.3.1 System Diagram and Equipment Description 

The master system flow diagram is shown in Figure 3-5. This diagram will be used as 
an aid in the Laboratory Procedures Manual to describe the system activities during 
various phases of calibration and alignment. In those applications of the diagram, interface 
lines will be added to indicate specific modes of operation. Data flows are indicated by 
narrow lines; dynamic or monitor interfaces are indicated by wide lines. A brief descrip- 
tion of each of the boxes represented on die master system diagram, m Figure 3-5, follows: 

Input/Output Console - The input/output console consists of the equipment that provides 
a manual computer interface. Included m the input/output console are* computer control 
panel, keyboard and typewriter, paper tape reader and punch, and the display panel. 

Operator — The operator in this system must perform many of the tasks of control and 
data transfer. The box "operator" includes not only the person(s) directing the laboratory, 
but also his worksheets, instructions, and notes. 

Systems Control and Monitor — This box represents the equipment, capability, and 
activity used to monitor and control the system during calibration and alignment. 

Frequency Counters — Six frequency counters are available for use in calibration to 
measure instrument output. These counters measure the number of counts on one pulse 
train for a fixed number of counts on another. One of the two trams may be a difference 
tram formed from two inputs. The frequency counters are used in calibration, because 
they can read the leading edge of one pulse tram and thus substantially reduce the 
quantization error relative to the use of the computer registers. 

Auxliary Data Sources - These include data sources available to the operator but not 
sufficiently well defined as equipment or measuring devices to be represented individually 
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Figure 3-5. Laboratory Flow Diagram 
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Oil the system flow diagram. Examples of these sources are bubble levels and survey 

E 

information on the magnitude of g and . 

Autocollimators — The two-degree -of -freedom autocolliraators are available to measure 
the earth-fixed coordinates of the ISU mirror normals. 

Resolvers — These resolvers measure the orientation of the test table. The angles 
and 0^ are static resolver readouts on the trunnion, outer and inner axes of the 
test table. The angle * s a rotary axis readout which can be used m either a static 
or dynamic mode . 

Master Oscillator — This is the central timing source of the system. The master oscillator 
includes countdown circuitry. 

Gyros, Accelerometers — These are the instruments contained within the strapdown ISU 
(see Section 2.2). 

Interface Electronics Unit (IEU)* — The IEU allows the computer to sample outputs of the 
inertial instruments and tuner. The IEU contains accumulating registers for each of the 
inputs shown m the diagram and the capability to periodically interrupt the computer to 
allow for sampling and resetting (without loss of data) of each of the registers. 

Computer* — The computer schematically indicated in the system diagram is the laboratory 
computer Honeywell DDP- 124. Other portions of the data processing shown may, however, 
be performed on other computers at the discretion of the programmers and operators. 
Blocks shown within the computer represent functions used m both calibration and align- 
ment, Shown are programs to input and output data from and to the console, a program 
to input data from the IEU, and memory buffers for input data and output data (the results 
of computations). Space has been left within the computer block to allow representation 
of the various data processmg tasks. 

3.3.2 Equipment Interface 


Figure 3-6 illustrates the principle data paths that might be of interest daring calibration 
and alignment. Each of the paths is numbered and described by number in Chart 3-5. 
These paths represent the calibration or alignment data flow. 


*More detailed descriptions are presented in Appendix C, 
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Figure 3-6. Laboratory Flow Diagram 
System Data Interfaces 
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CHART 3-5 


Description of Equipment Data Paths 


Data 

Interface # 

Type 

Data Description 

1 

Coded 

Input — Includes magnetic and paper tape, keyboard, 
and display and computer control panel inputs. Data 
represents information and control from or through 
the computer plus data or program filed on magnetic 
or paper tape. 

Output — Display, typewriter, or paper tape panel 
output to the operator. 

2 

Visual or tape 

Various types of data 

3 

Visual 

Two angles/autocollimator 

4 

Visual 

Four angles test table position 

5 

Pulsed 

Rotary axis motion probably one pulse 

6 

Timing 

Frequency less than 1 MHz 

7 

Timing 

2. 034 MHz 

8 

Timing 

3. 6 KHz 

9 

0, +1 or -1 pulses 
at 3. 6 KHz 

Gyro output to frequency counters and IEU 

10 

Pulsed 

Zero crossing pulses from each of two vibrating 
strings per accelerometer to frequency counters 
and IEU 

11 

Count 

For inputs shown — count would be number of time 
pulses per n turns 

12 

Count 

For input shown - count would be number of gyro 
pulses (signed) for n turns of table 

13 

Count 

For inputs shown - count would be number of time 
pulses per ini gyro pulses 

14 

Count 

For inputs shown — count would be numbers of time 
pulses per N accelerometer pulses 

15 

Visual 

Status and monitor information plus output from 
the counters 

Low visual display and printed 

16 

Manual 

Manual input of data to computer (includes key 
punching, mounting of tapes and punching of 
buttons) 
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CHART 3-5 


Description of Equipment Data Paths (Continued) 


Data 

Interface # 

Type 

Data Description 

17 

Binary Data 

Counts from IEU registers. Input m succession 
with data valid for same period of time. 

18 

Binary Data 

Counts from IEU registers. May be summations 
of data from several successive transfer across 
interface #17 

19 

Binary Data 

Input data shown in buffer 

20 

Binary Data 

Output data shown in buffer 
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Manual and monitor interfaces are shown in Figure 3-7. The manual interfaces correspond 
to operator activities during various portions of the calibration and alignment procedure. 
Monitor is performed during the many procedures to verify the operation of equipment being 
used at that time. Explanations for each interface on Figure 3-7 are presented on the 
figure. 
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SECTION 4 

DEVELOPMENT OF CALIBRATION TECHNIQUES 


The functional description of calibration m Chart 3-1 in Section 3. 1. 1 indicated inputs of 
instrument data and environment measurements to a computation routine, which in turn 
outputs the calibration constants. The obvious sources of the relationships contained 
within that routine are the instrument mathematical models. In Section 4. 1 we will use 
the instrument models to develop general equations from which the relations found in the 
computation routine are evolved The general equations will be seen to contain certain 
controllable parameters which describe the input environment in terms of the test table 
controllable orientation and angular speed. In Section 4. 2 a number of sets of particular 
values for the control parameters will be chosen such that the general equations reduce 
considerably inform. Each set of particular values corresponds to a different calibra- 
tion "Position*'. It will be shown that the determination of any constant can be accom- 
plished by the simultaneous solution of at most two of the reduced equations. In 
Section 4. 3 the complete set of calibration computations will be delineated The relations 
employed in the computation routine correspond to the solution of the calibration con- 
stants from the data gathered in Positions 1 through 15. 

In Section 4. 4 we will describe the operations and computations required prior to the 
collection of calibration data. An example of a precalibration operation is the determina- 
tion of the orientation of the ISU body axes relative to the test table inner-gimbal frame 
(i. e. , T ). In Section 4.5 a brief discussion of the implementation of the proposed 
techniques will be presented. The discussion of the implementation of the calibration 
techniques is directed towards clarifying the relation between the developments in this 
document and the operational procedures described in the Laboratory Procedures Manual. 

Before proceeding to the development of the calibration techniques, it is appropriate to 
describe those incentives which motivated our specific choices of calibration techniques: 

• The determination, of any calibration constant should be made as insensitive as 
possible to the imprecision of any other constant(s). 

• It is advisable to use as few different test table orientations as possible; and 
where the orientations are different, to try to make the orientations differ from 
one another by as little adjustment of the table as possible. The satisfaction of 
this requirement serves two purposes. First, it will allow for the simultaneous 
calibration of many instruments. Secondly, by limiting the number of table 
orientations, the amount of manual activity will be limited, thus minimizing 
calibration time, and also the chances of human error. 
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• The imprecisions of the test table orientations should have a limited influence on 
the values of the calibration constants. 

• The calibration should be made as flexible as possible. We wish to present the 
calibration in such a way that additional experiments can be accomplished with a 
minimum number of changes to existing calibration procedures 

• We wish to accomplish our calibration with little or no data filtering. We would 
like to minimize the effects of noise by judicious choices of approaches other 
than involved software processing. 

• The computation program should be as simple as possible. 

• Data collection time should be limited to about 10 minutes, and the total calibra- 
tion time to less than eight hours. 

• The precision of calibration, as a function of time, should be apparent from 
error analyses accomplished on the resultant techniques. 

In the discussions which follow, it will be found that it is possible to satisfy a majority of 
the above requirements. 

4. 1 DEVELOPMENT OF GENERAL CALIBRATION EQUATIONS 


The routine that accomplishes the evaluation of the calibration constants we indicated in 
Chart 3-1 as a routine entitled ‘'Computation of Constants". In this section we will develop 
the general equations from which the computational routine is developed Those general 
equations will be seen to contain the control parameters which describe the environment 
inputs. The chosen control parameters are the angular speed of the test table, the first 
two gimbal angles (e j and 9 ^) of the test table, and the T matrix (which is a function 
of the Og and $4 gimbal angles), hi Section 4. 2 we will show how particular choices of 
these control parameters result in relations from which the calibration constants can be 
extracted. 

The presentation in this section is divided into two parts. In Section 4. 1. 1 we develop 
the general equations for the three system gyros, and in Section 4 1. 2 we develop the 
general equations for the three accelerometers. 

4. 1. 1 Gyro Equations 


The development of the general equations begins with a presentation of the Fundamental 
Gyro Model. After introducing the ERC environment and geometry into that model we 
will have developed equations which are a function of, among other things, the angular 
speed of the test table and test table orientation parameters. In subsequent subsections 
we will show how the control of those test table parameters is employed in the determina- 
tion of the required instrument calibration constants 
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The presentation m this subsection contains a great deal of mathematics So as not to 
interfere with the prose, we will present the mathematical development as a series of 
eight charts. In each chart, after the first, the equations found on the preceding chart 
will be modified to indicate certain assumptions about the environment In the dis- 
cussions which follow we describe, in turn, the assumptions and the related mathematics 
presented on each chart. 

The Fundamental Gyro Model (Chart 4-1) 


The development of the generalized gyro calibration equations begins with the gyro 
mathematical model. That model was described m Section 2. 2. 3 and presented as 
Chart 2-2. Chart 2-2 is repeated here as Chart 4-1. The gyro mathematical model 
describes the relationship between the output (So) of the gyro and the input kinematic 
environment (a and u) over a time period Iq-4:^. 

Introduction of Laboratory Environment (Chart 4-2) 


We first introduce into the gyro mathematical model the vector representation of the ERC 
laboratory kinematic environment At the top of the chart the kinematic inputs are listed 
Note that every possible input has been listed. This is done so that, at one point in the 
development, there exists an expression which assumes nothing about the negative effect 
of any possible input. Note also that the environment description assumes that the gyro 
is subjected to an input angular velocity in * as generated by the test table rotary axis 
motor 

Approximations (Chart 4-3) 


The next step is to neglect those kinematic inputs to the gyro which can reasonably be 

expected to have a negligible effect on the gyro output. A gyro is designed to be nominally 

a linear angular velocity measuring device; therefore all acceleration-sensitive terms 

are small. The noise acceleration and the test table-induced accelerations are also 

small relative to the nominal g input. Therefore the effects of these small accelerations 

are second order in all unbalance and compliance terms and are assumed negligible. 

Similarly, it is assumed that the small angular velocity noise terms can be neglected in 

all angular velocity- sensitive terms other than the linear term Note that in Chart 4-3 

E T 

we have arranged the equations such that only the deterministic g, w , and w inputs 
exist to the right of the equality. Note also that the equation has been divided by the 
gyro scale factor. 
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Introduction of Bocfy Axes and Instrument Indexing (Chart 4-4) 


The equations presented thus far have referred to a single gyro. At this point we intro- 
duce the index i defining three gyros (i = 1 , 2, 3). 

We mentioned in Section 3.2 that for calibration purposes, the vectors g, u> , and w 
must be known in body-axis components. In Chart 4-4 we introduce the body axes, and 
express the integrals of these vectors in terms of body-axis components. Because the 
nominally known instrument- to-body- axis transformations differ from the actual trans- 
formations by small numbers, it is assumed that the nominal values can be used in other 
than the proportional angular velocity term. Note that the first three elements in the 

fh pi 

equation represent the i 111 row of the (Q ) matrix, scaled by the gyro scale factor 
Note also that these elements are assumed constants whereas, on a microscopic scale, 
they are time-varying within the limit cycle amplitude of the instruments. 

The function found to the left of the equality has, at this point in the development, been 

p 

defined as the triple Ppi = 1, 2, 3). This vector, which we will refer to as the gyro 

processing vector, contains the instrument readout term plus the quantization and noise 
G 

terms. The vector, and the approximations made in its evaluation during calibration, 
will be discussed in Section 4. 3. 1. 

Choice of Body Axes (Chart 4-5) 

The next step in the development is to introduce the ERC ISU nominal geometry. Those 
transformations which describe that geometry were defined in Section 2. 2. 1. Because 
the orientations of the output and spin axes are not cyclic, a general index equation 
cannot be developed. Therefore a separate equation for each gyro is presented in 
Chart 4- 5. 

At this point the general equations relate the gyro processing vector (which includes 
measurable gyro readout, noise, and quantization error) to the measurable body-axis 

TP T 1 

components of the environment (described by g, uo , and w ). In the remaining charts 
we will relate the body-axis components of the environment to the controllable test 
table parameters. 

Integral Evaluations (Charts 4-6, 4-7, and 4-8) 

The preceding chart listed the required equations as three expressions which are linear 
in the unknown calibration constants. The coefficients of those unknown constants are 
presented as integrals of body- axis components of _g, u , and u . In the discussions of 
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calibration techniques which follow this subsection we will show how calibration is 
accomplished by a control of the values of these integrals. Before we can show this we 
must functionally relate the integrals to the controllable test table geometry. The 
following three charts develop those functions. 

TP T 1 

In the first chart (4-6) the integrals of the body-axis components of g, _to, and cu are 
expressed as transformations from the frames in which they are well known. The 
definitions of the transformation geometry and notation were explained in Section 3. 2. 1. 

In Chart 4-7 the earth- axes components of g and to and the rotary axes components of 

n-i 'T'Tc 1 _ m 

w are introduced. Additionally, the T and T RT matrices are expressed as functions 
of the and ^ gimbal angles. (See Section 3. 2. 1 for definitions of this geometry ) 

With these equalities introduced, we can now extract all but the time-varying parameter 
(0g) from the integrands of the equations 

In the final chart (4-8) the integrals are combined as the calibration constant coefficients. 
The equalities listed at the top of the chart allow the integrals to be separated into sums 
of monotonically increasing terms, harmonic terms, and terms which are functions of 
terminal conditions only. The only harmonic terms are those which contain integrals of 
sines and cosines. The monotonic increasing terms are those containing At, and the 
terminal condition terms are those which contain A’s other than At. 


Charts 4-5 and 4-8 constitute the required general gyro calibration equations . Note that 

our result is a set of three functional relationships among: the ’processing vector” 

G E 

P^, the unknown calibration constants; the magnitudes of gravity (g), earth rate (w ), 

and latitude (\); the total time of integration (At), and the controllable test table 

parameters, which are 


and 


0^ — The trunnion axis angle 

^2 — The total angle of revolution about the rotary axis 


T d0g 

w = The speed of the test table 

dt 

T - The matrix which transforms from the rotary axis frame to the 
body axes for the m™ calibration position. 

(See Section 3 for definitions of all test table geometry, ) 


In Section 4. 2 we will show how gyro calibration is accomplished by a control of the test 
table parameters contained within the functional relationships. 
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CHART 4-1 


THE FUNDAMENTAL GYRO MODEL 


THE GYRO MODEL IS 


E5 k - / N (us.G) dt + / N R+Bj(a-G) + B Q (a-0) + R g (a-S) + C n (a-G) 2 + C gs (a-S) 2 

- k= l J tfi t n L 


c Ig(a-G) (a-S) + C os (a- O) (a-S) + C IQ (a-G) (a-O) 


+ Q n ( to* G) + Qjg(co-G) (to*S) + J 


dt <ar 2) 


dt + An + Eq 


u- is the angular velocity applied to the gyro 
a is the acceleration applied to the gyro 

tg s t £ t N is the time interval over which a and a. are measured 
= Nr, where N is an integer, andr is the gyro sampling 

period 

S is a unit vector along the spin axis of the rotor 

O is a unit vector directed along the output axis as defined by the 
gimbal 

G is a unit vector along O x S (that is, the sensitive axis of the gyro) 

6^ is the kth gyro pulse, equal to +1, -1, or 0 for positive, negative, 
or no pulse 

is the gyro scale factor 
R is the gyro bias 

Bj Bq and B g are the gyro unbalance coefficients 

^II ^SS ^IS C OS and ^ro are the ® yro COrn Pl iarlCe coefficients 

Qjg and Qjj are dynamic coupling coefficients due to gimbal deflection 
and scale factor nonlinearity, respectively 

J is the angular rate coefficient 

An is the effect of gyro noise over the [tg, t N l interval 
Eq is the gyro quantization error 
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CHART 4-2 


INTRODUCTION OF LABORATORY ENVIRONMENT 



t TT EE T T E E 

• ^ = w + w = (« + A ^ ) + (a- + Aw ) is the total applied angular 

velocity 

TT 

where: w is the total test table angular velocity 

T 

te is the measured test table angular velocity 
EE 

is the total laboratory angular velocity 
E 

a is the assumed (surveyed) laboratory angular velocity 
t TT TT • TT 

« a = g + u- x (w x r) + -o x r + Aa is the total applied specific force 

where* g is the assumed (surveyed) laboratory specific force 
TT TT 

a x (**. x r) is the centripetal specific force due to the table 
motion ! 

■ TT 

w x r is the angular rate specific force due to the table motion 

Aa is the deviation of the assumed laboratory specific force from 
the true 

TT TT * TT 

• a: x (to x r) and to x r are formal expressions, as r is not explicitly 

defined. 


INTRODUCING THESE ASSUMPTIONS INTO THE EQUATION IN CHART 4-1, WE HAVE 


A « jy^kj = I t ' N |t E + i T + A u E + A u. T ) gJ dt 

+ j jj^R + Bj J^(g + jjj TT x (o> TT x r) + t x r + A a) oj + jjs * x (fc TT x r) + C"*' xr+ Aa) pj 

+ Bg |fe+ « IT ii(i TT \rlu Tr 5i* Aa) S + C u jfe + S; TT x(fc TT x r)+ x rT x r+ Aa) pj 2 
+ C ss ^t« TT x( t TT xr)tJ T 'ii*lj| S 2 

+ Cjj (g + a> TT x (ij* T x r) + r + Aa) G j j(g+ x TT x {t^ TT x r) + (u^ T xr + Aa)§| 

+ C QS (g + x TT x (a TT x r) + x TT x r + Aa) pj j^(g + (a TT x r) + u.' TT x r + Aa) sj 

+ C| 0 (g * u TT x (o. TT x r) + C fT x r + Aa) pj j(g+ u TT x (x TT x r) + o? TT xr + Aa) p| 

+ Qjj (u.® i- o. T + A w £ ->■ A i T ) G 2 

+ Q b (a® + a T t fi u. E t A J) G |(u. E + + A a® + A <;. T ) 3 j + J d/dt jju.® + a T + Aa E + A J) pj | dt 


+ An + Eq 
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II M2 


APPROXIMATIONS 

ASSUME 


CHART 4-3 


• Terms containing the integral of components of A a always have negligible 
effect on the gyro readout. 

• Terms containing the integral of components of u- x (u> x r) and 

♦ *p r r 

x r have a negligible effect on the gyro readout. 

1 « TT x (w TT x r) | < (0. 22) 2 x 0. 5 = 0. 024 ft. /sec. 2 = 0. 008 g. 

E 

• Terms containing Qjg or Qjj and the integrals of components of A a? or 

A have a negligible effect on the gyro readout. 

E 

• Terms containing J and the integrals of the rate of change of A cu and 

T 

A u have a negligible effect on the gyro readout. 


INTRODUCING THESE ASSUMPTIONS INTO THE EQUATION IN CHART 4-2 WE HAVE 
(AFTER SOME ARRANGEMENTS): 
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CHART 4-4 


INTRODUCTION OF BODY AXES AND INSTRUMENT INDEXING 


ASSUME I • B^ Bg Bg are a triad of orthogonal unit vectors which have a fixed 

orientation relative to G • and S x where 1 is the instrument index equal 
to 1, 2, or 3. 

Therefore, ^ = S (G^ B^ , ^ = 2(0-^)^, ^ = S^-B^ 

k k k 

* (O^Bj.) and (S B^) differ from their nominal values (G 

(OpB^) 11 and (SpB^) 11 by small numbers. Those differences are sufficiently 
small as to only affect the gyro readout via the proportional angular velocity 



(G [ Bjl/AC 


(Gj 


(G, B 3 ),A* 


B 0/A4>[, 


C U' A *|| 


c icy i4 [i 


c os'«[i 


V^|i 


WHERE 


; N («. E B,;ai » r‘ N i« T B.jdt \ 

*o <0 ‘ 1 

; ,f V E § 2 idt - ;' N ( W T B 2 >dt j 

‘o ’*0 ” * 


J jj <2i 5k>" j*” ( £ 5k* dt | 

|rE(G 1& ) ! '(GiB r )''- ,N (gB k )(gB r )d t l 
1 k r t Q J 

'EB{S 1 B 1> ) n (S 1 B r ) 11 ‘ N (ft BJdt | 

,Kr l o ► 

jEE(GiB k )''(q i B r )"-' N ( g B k )( E B r )dt| 

( K r t Q 1 

|S£(5 | B k ) n (S 1 B r ) n ^( g B k )( g B r ) dt J 

l£E(0 | | k ) n (S 1 B/- t «( s B k )( g Br)dt} 

|EE(G| Bp” ( Sl gp" (y Bp^ B r )dt + r‘ N (a E B k )(t T B r )dt + >< t T Bpfco 6 Bp, it. ;‘% T B k )( 2 T B r )dt } 
] f T ;<- r A 2*'" <§. 5 r >“ (V 5 §*><i E &><“■ ;V B k )(- T B r )dt. r tN (u. T gp^E gp* . 9<><± T 5P* j 

)E(2,B k ) n ;' N d - t ( ffi E B k )Mt T B k ) d.[ 


{ U \ i - (1/ impart, - (IMdjpEqli - (1/M), f*** (a^E . Au T ) Cj dt 
l L “ J t 0 
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CHART 4-5 


[Tool |C 0 T1 p -1 < 
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CHART 4-6 


INTEGRAL EVALUATIONS 
ASSUME 


• g and ^. T are located as shown m Figure 3-2 

• AH transformations between frames follow the convention given by the 
following example 1 

Bj - £ T BE R 

. m im 

• | N dt = t N - t Q +€£<= At + € c where g c is the clock quantization error 
t 0 

• The effect of € c is negligible 



THE NINE INTEGRALS CAN NOW BE WRITTEN AS: 


( yE . dt 
*0 " " 

= J^E EE t? R T RT T TE (w E .E ) dt 
J + mnp im mn np v _ *:rp' 

t 0 

A (^ T .b ) dt 
to ~ ~ 

= t RR (to T .R ) dt 

J , m im — m ; 

*0 

dt 

to 

= f* 1 * £££T BR T RT T TE (g.E ) dt 
J , mnp im x mn nps -p' 

to 

to 

= r EEEt? R T RT T TE {g.E ")j js£E T BR T RT T TE (g. E )"] dt 

J t [mnp 1 km mn np 'l_pj |qus rq qu us _s'J 


- t* 1 * Is ££ t BR T RT T TE (w E .E ) fs EE t BR t rt t te (oj E .E )1 dt 
[mnp i km 1 mn A np _pj [qus A rq A qu us v _ _s'j 


= [sEE I® R T RT T TE (w E .E ) S t BR (w T .R )1 dt 

[mnp 1 km i mn J -npT -pj [g rq v _ _q'J 

r 0 

^“ T -?k^" T *5r) dt 

- (“ T -V] dt 

0 L 

,’tN 

% dt 

to 

- tN-to" 44 

^dt[<“ E -5kH“ T -5k>] dt 

to 

- E EE T BR T R, f Tf E (co E E 1 + E T BR (cof R R 

* mnp T ktn T mn T npt“ -jy + ^kq^ .K q ;j 
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CHART 4-8 


INTEGRAL EVALUATIONS (Continued) 


= 1/2 - (cos 2}> 2 ) / 2 
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= & sm (dj-t- X) 
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^ T rl i k3 * T kl X r3 } SU1 2 u rl A k2 kl X r2 ; 


T W HiaU1 *2 + T k 3 R Acos0 2 


whore 0. = 1 if k = r 



4. 1. 2 Accelerometer Equations 


As with the preceding gyro equation development, we develop the general accelerometer 
calibration equations by introducing the laboratory geometry and environment into the 
instrument model. A great number of the comments relative to the development of the 
general gyro equations apply equally well to the development of the general accelerometer 
equations. 

As in Section 4. 1. 1 , we discuss each of the charts in turn. 

Fundamental Accelerometer Model (Chart 4-9) 

The accelerometer mathematical model was introduced in Section 2. 2. 2. Chart 4-9 is a 
repeat of Chart 2-1, showing the input/output relationships for a vibrating- string accel- 
erometer. The notation presented is self-explanatory. 

Introduction of Laboratory Environment (Chart 4- 10) 

In Chart 4-9 the accelerometer output is seen to be influenced by only applied accelera- 
tion inputs. Note that the accelerometers are assumed to be in a stationary attitude 
relative to the earth. The stationary attitude assumption dictates that all accelerometer 
calibrations will be accomplished without a use of the dynamic rotational ability of the 
test table. The main reason for this constraint is the fact that a motion of the test table 
introduces undesirable angular velocity-related accelerations. (See Section 2. 1 of the 
trade-off document.) 

Approximations (Chart 4-11) 

The environment approximations are self-explanatory. All neglected terms are assumed 
to have a second order effect on the accelerometer readout 

Introduction of Body Axes and Instrument Indexing (Chart 4-12) 

The comments presented in the gyro equation development apply equally well here. 

A Choice of Body Axes (Chart 4-13) 

These equations are the desired general form. In Section 4. 2 the determination of the 
calibration constants will be shown to be dictated by a control of the parameters found in 
the Environment Evaluation part of Chart 4-13. 
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CHART 4-9 


THE FUNDAMENTAL ACCELEROMETER MODEL 


THE ACCELEROMETER MODEL IS: 



a is the acceleration applied to the accelerometer 

t ^t^ t K is the time interval over which a is measured 

A is a unit vector directed along the input axis of the accelerometer 

N^ and Ng are the number of zero crossings detected in t ^tst^ 
from both strings of the accelerometer 

Eq is the instrument quantization error due to the fact that t and t^ 
do not correspond to zero crossings 

D^ is the accelerometer scale factor 

Dq is the accelerometer bias 

Dg is the second order coefficient 

Dg is the third order coefficient 

fg and f^ are string frequencies in pulses/ second 
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CHART 4-10 


INTRODUCTION OF LABORATORY ENVIRONMENT 
ASSUME 

• a L = g + Aa is the total applied specific force 
Where : 

g is the assumed (surveyed) laboratory specific force 

Aa is the deviation of the assumed laboratory specific force from 
the true 

INTRODUCING THESE ASSUMPTIONS INTO THE EQUATION IN CHART 4-9 WE HAVE: 



(THAT IS to TT x (w TT x r) and i) TT xr WILL NEVER BE SENSED) 
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CHART 4-11 


APPROXIMATIONS 



• Terms containing the product of the integrals of components of Aa with 
Dg or Dg have a negligible effect on the accelerometer readout. 

• fe* A) is a constant over the time interval t £ t 5 t b 

• J t& dt = At + € 

fc 0 

where : 

At = Nr, N is an integer and r is the clock period 
is the clock quantization error 

• Terms containing the product of €q with Dq, Dg, or Dg have a negligible 
effect on the accelerometer readout. 


WITH THESE ASSUMPTIONS, THE ACCELEROMETER EQUATION BECOMES: 
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CHART 4-12 


INTRODUCTION OF BODY AXES AND INSTRUMENT INDEXING 



* B-p B 2 , and Bg are a triad of orthogonal unit vectors which have a constant 

orientation relative to A., where 1 is the instrument index equal to 1, 2, 
or 3. 

Therefore, 3 

A, “*l(A 1 'B k )B k 

* (A^B^.) differ from their nominal values (A.'B^) n by small numbers. Those 

differences are sufficiently small as to only affect the accelerometer readout 
via the proportional acceleration term. 

* The effect of ^ is negligible. 


pf = ) (g - b 1 )&t [ 

+ } i'£-5 2 )Al! 

+ {Blass' }. j(l-B 3 )At| 

+ { D i D o}, l at ! 

+ { DlD2 h \[l r ( ^'5k )n( ^-5rt n (£-B k )(g-B r )]At J 

+ KM, i[kr m ( - l '^ )n( ^‘^ > ^'^ )n(£ '^ K£ ‘- r )(S 'M At i 


if = (N 2 - N x ) -LD^J. J* (Aa-A x ) dt + Eq 


• The first three terms in the right hand side of the equation include the 
effect of the misalignment. 

• The second order cross couplings due to the misalignment have been 
neglected by the second assumption. 
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CHART 4-13 


A CHOICE OF NOMINAL INSTRUMENT AXES 


ASSUME 



pf = {D 1 (A 1 .B 1 )} i {( g .B 1 )4t}pA. {^-^{(g-B^At} P A = {D^.B^^tg.B^At} 

- {DitAj-Bj)} + {Dj^-B^J^tg-B^at} + (d^Aj-^) } 3 {(g-^jat} 

+ 'DjfAj.Bj)} {(E-B 3 )flt} + {(g.BjJAt} + {c/Aj-y} {(g-BjlAt} 

i Z V 


+ ! At l 

+ {<sV At } 

- { D l D 3} 1 {W 3 "} 


+ { d i d o}, l 4t i 


i D l D o}, 


{ D ft } 2 {(I-B 2 ) 2a ‘} + { D 1 D 2 } 3 {(I'l/M 

{ d i° 3 } 2 {<£-52> 3a *} + K d 3 }, {(i^M 


ENVIRONMENT EVALUATION 


ASSUME 


• g is located as shown in Figure 3-2„ 

• All transformations between frames follow the convention given by the 
following example: 

B = S T BR R 
—l m lm — m 


THE 1 th BODY AXIS COMPONENT OF g IS. 


(g-B.) = S£S T BR T RT T TE (g«E ) 
m n p im mn np v a tip' 

(Test Table readout) 

= SS T? S T SE (g°E ) 
m n im mn v a _n' 

(Autocolimator alignment) 
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4. 2 CHOICES OF CALIBRATION ENVIRONMENTS 


The general gyro equations are presented in Charts 4-5 and 4-8 and the general accel- 
erometer equations are presented in Chart 4-13. These equations represent the functional 
relationships among the instrument outputs, the input environment measurements, and the 
calibration constants. 

Each instrument equation is linear in n unknown calibration constants. Therefore, it is 
possible to determine the numerical value of all calibration constants contained within any 
equation from the simultaneous solution of n equations corresponding to n different 
measurements of instrument outputs and input environments. Such a technique of con- 
stant determination would involve the inversion of an n x n matrix When n is large, as 
it is in these instrument equations, matrix inversion is very cumbersome. 

There is, fortunately, an easier technique for determining the calibration constants 
That technique involves the control of the environment inputs (by a control of the test 
table parameters) such that the instrument outputs would be insensitive to a large number 
of terms. This corresponds to the adjustment of the environment- sensitive coefficients 
of a large number of constants in the general calibration equations to zero. If it were 
possible to null all but one, the determination of the remaining constant would naturally 
be trivial. In this system, however, it is not possible to null all but one but we can m 
many cases null all but a few coefficients. In the subsections that follow we will apply 
this "nulling technique" to the calibration of the ERC ISU. The result will be a set of 
equations from which any calibration constant can be determined by the simultaneous 
solution of at most two equations. Each equation corresponds to the input/output rela- 
tionship for an instrument subjected to a particular environment, by control of the test 
table parameters. 

We begin our presentation, m Section 4. 2. 1, by dictating the environments and developing 

p _ i 

the equations from which the gyro scale factor and (Q ) matrix can be determined. In 
Section 4. 2. 2 which follows we will show how to calibrate the gyro unbalance, bias, and 
square compliance terms. In Section 4. 2. 3 we will show how to determine the com- 
pliance-product coefficients, and in Section 4. 2. 4 we will complete the discussion of 
gyro calibration by describing the experiments for investigating the gyro scale factor non- 
linearity and J term. The discussion of the calibration of gyro constants m any sub- 
section will assume that the constants discussed in previous subsections are well known 
from previous calibrations. 

The description of accelerometer calibration begins in Section 4. 2. 5 with a description 
of the calibration of all but the cubic term, hi Section 4 2. 6 we complete the calibration 
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developments by describing the determination of the remaining accelerometer cubic 
terms. 

4.2.1 Determination of Gyro Scale Factor and Misalignment 

In this subsection we show how particular choices of input environments reduce the general 

gyro equations to forms which enable a relatively simple calibration of the gyro scale 

G -1 ^ 

factors and elements of the (Q ) matrix. (The matrix elements are sometimes referred 

to as "misalignments" from the nominal ISU design ) Our attention is directed to the 

general gyro equations found on Charts 4-5 and 4-8. We will dictate choices of the test 

table parameters found in the integrals shown in Chart 4-8 such that the desired angular 

velocity- sensitive terms predominate. 

We see that many of the integrals found in Chart 4-8 are functions of harmonic terms as 
well as terms which increase monotomcally with time The harmonic terms are terms- 
involving integrals of trigonometric functions of dg. Such integrals are bounded m value; 

m d0^ 

as a matter of fact, if = can be made constant, the harmonic terms would equate 

dt 

to zero for any multiple of whole turns (0g = 2mr) of the table. Under such conditions a 
large number of the terms in Chart 4-8 would disappear. In Chart 4-14 we see the sub- 
stitution of the integrals into Chart 4-5 under the condition of whole turns of the table, 
while rotating at a constant speed. {See Section 2. 1 of the trade-off document for further 
comment about whole-turn equations. ) The assumptions made in the equations in 
Chart 4-14 are shown at the top of the chart. The condition on the transient terms re- 
quires additional comment. 

The ERC table will have a precision limitation on its ability to rotate at a constant speed. 

rp n-« w A r T l 

That limitation is two parts in ten thousand, that is Aw /w ^ 2 x 10 , where Aw is 

the error in the speed of the table, and to is the speed of the table. Assuming that a 
maximum error of plus £ is evidenced m a first half turn, and a maximum error of minus 
C is evidenced in a second half turn, then 

27 T 

!' sin dt = 

0 " 


(* sin(w^ + A w T)tdt t f sin(w^ - Aw^Jtdt 

0 7T 



4-21 



Daring all scale factor and misalignment determination experiments the maximum possible 
speed will be used. That speed will be just below the saturation level of the gyros, which 
is 15 degrees per second. (See Section 2. 1 of the Trade-Off Document for further com- 
ments. ) Under such conditions: 

• The proportional transient terms go as 

co 3 A = 120 x 10 deg = 0. 04 sec (per revolution) 

• The unbalance transient terms go as 

BgA = 8 x 10“ 7 deg 

• The compliance transient terms go as 

Cg 2 A = 0. 32 x 10“ 7 deg 

which are all obviously very small and can be neglected (under the assumption, of 
course, that the above analysis typifies the worst-case deviation from a constant speed). 

Referring again to Chart 4-14, we see that a horizontal position of the test table rotary 
axes (i. e. , ^ = 90°) would null all imbalance terms. Chart 4-15 introduces that condi- 
tion The remaining test table control parameters in Chart 4-15 are the first column of 
the T BR matrix. (The first column of the T BR matrix dictates the orientation of the 
table rotary axes (R^) relative to the ISU body axes. ) The orientation of is a function 
of the inner and outer gimbal angles (0, and 5h). Having two gimbal angle degrees of 
freedom dictates that any values of the first column of T can be requested. (Equating 
the T BR choices to and settings is the subject of Section 4. 4 2. ) In Charts 4-16, 

4-17, and 4-18 we show the calibration equations for six choices of the first column of 
T . All choices are shown at the top of the charts. We see that 

• Chart 4-16 corresponds to the alignment of the first body axis with the rotary 
axis in both the plus and minus sense. 

• Chart 4-17 corresponds to the alignment of the second body axis with the rotary 
axis in both the plus and minus sense. 

• Chart 4-18 corresponds to the alignment of the third body axis with the rotary 
axis in both the plus and minus sense 

In Charts 4-16, 4-17, and 4-18 the test table parameters have been completely specified 
The first table gimbal angle is equated to 90°. The second gimbal is rotating over 
whole turns at a constant speed. And the third and fourth gimbals are implicitly specified 
by choices of the first column of t BR . We will refer to these six orientations as Posi- 
tions 1 through 6, respectively. 
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We note in Chart 4-16 that only the first terms in the three gyro equations have opposite 
signs in Positions 1 and 2. Therefore, those terms can be isolated as a function of Pjjf 
(k = 1, 2, 3) by simply subtracting the equation for Position 1 from the equation for Posi- 
tion 2. From known environment inputs and known Pr 1 values the first gyro constants can 

r* K. 

therefore be determined. (The manner in which Pjj? is known is the subject of Section 
4. 3. 1. ) The second and third constants can be found by similar uses of Charts 4-17 and 
4-18, respectively. The final equations for the gyro scale factors and (Q u ) elements 
are found in Section 4. 3. 2. As a matter of convenience the computational equations from 
this and subsequent subsections are tabulated at a single point in this document, which is 
Section 4. 3. 2. 

4. 2. 2 Determination of R, Bp Bq, B g , Cjj and C gg 

Subsequent to the calibration of the principal angular velocity sensitive terms (scale 
P -1 

factors and (Q ) elements), the gyro equations predominantly contain, as unknowns, 
acceleration- sensitive terms (l. e. , unbalance and compliance coefficients). This pre- 
dominance is even more evidenced when the angular velocity input is controlled to a small 
constant value. Under that condition the remaining angular velocity terms (Qjj and Qjg) 
become relatively unimportant as influences on the gyro outputs. These points suggest 
that the calibration of the unknown unbalance and compliance coefficients should be ac- 
complished under the conditions of extremely small angular velocity inputs. Not only 
will the Qjj and Qj g be negligible, but also the imprecision in the already calibrated scale 
factor and (Q^) ^ elements will have a minimum influence on the precision of the unbalance 
and compliance coefficients to be determined 

In Chart 4-19 we present the general gyro calibration equations under the influence of the 
minimum practical angular velocity environment. That angular velocity input is earth rate, 
that is, the table is stationary relative to the laboratory. We say minimum 'practical" 
environment because it would be possible to rotate the table at near minus earth rate, thus 
reducing the total angular velocity input below earth rate, but earth rate alone is so small 
that there appears to be no reason to try to regulate the speed of the table to a small number 

At the top of Chart 4-19 we present our gimbal angle choices of Oj = 0 and ^ = 90°. There 
are several reasons for these choices. First, it must be pointed out that we are interested, 
for the purpose of calibrating acceleration-sensitive coefficients, in controlling the 
orientation of only the input £ vector relative to the body axes. To completely control one 
vector relative to the body axes requires only two orientation degrees of freedom. Two of 
the four test table degrees of freedom can therefore be chosen for matters of convenience. 
We choose the particular values of 6j and as shown in Chart 4-19, for the following 
reasons of convenience: 
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CHART 4- 14 


1) All transient terms are negligible 

2) Terms of 0(ctF)^ are negligible 

3 ) Terms of CKu^Atfg) are negligible 
4 } <t>2 15 constant 

For whole turns of ^ the gyro equations are 
given on this page. 
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- . CHART 4-16 

GYRO CALIBRATION WHOLE -TORN EQUATIONS POSITIONS 1 AND 2 j 
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GYRO CALIBRATION WHOLE- TURN EQUATIONS POSITIONS S AND 4 
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CHART 4-18 


GYRO CALIBRATION WHOLE -TURN EQUATIONS POSITIONS 5 AND 6 
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• The value <t>\ = 0° is chosen because it results in the same 03 and 04 settings 
for the six positions required for the subject calibration as required for the first 
six positions. 

• The value - 90° was chosen for two reasons, the first reason is that it places 
the 03 gimbal in the north-south direction, and that gimbal, m conjunction with 
tiie east-west <t>i gimbal, can be used for small angle corrections of the table 
base motions as measured by bubble levels (see Section 4. 4. 4), The second 
reason for the choice of 09 = 90° is that it results in only the second column of 
TBR being required xn calibration computations* This results in a minimum 
amount of data handling during precalibration survey activities. 


The equations m Chart 4-19 contain, as control parameters, only the first and second 
column of the T^® matrix. The acceleration- sensitive terms are, however, a function 
of the first column only. The table orientation control will therefore be preoccupied with 
that column. In Charts 4-20, 4-21 and 4-22 we introduce the choices for the first column 
of T corresponding to Positions 7 through 12. These choices for the first column of 
T are the same, respectively, as they were for Positions 1 through 6 . With the 
assumption that the first three gyro coefficients are known from the calibration described 
in the preceding subsection, we see that Charts 4-20, 4-21 and 4-22 present six equations 
m the six unknowns: 


R > B p b q, Bg, Cjp Cg g 

We note that at most two equations are required for the solution of any required unknown. 
In Section 4. 3 the solution of the equations for the six unknowns is presented. 

4. 2. 3 Determination of C IO , C Ig and C og 

The three product- compliance coefficients (Cjq, Cj g and Cqq) were not evidenced in any 
equation for Positions 1 through 12. None of those positions senses the minimum of two 
body-axes components of acceleration required for the detection of product-compliance 
coefficients. In this subsection we choose three additional laboratory-fixed orientations 
(Positions 13, 14 and 15), each position detecting two (and only two) body-axes components 
of g. On Charts 4-23. 4-24 and 4-25 we present the instrument equations for those three 
positions. In each equation on those three charts there exists only one unknown product- 
compliance coefficient. The solution for that coefficient, m terms of the known input 

r* 

environment vector P , and the previously determined calibration constants, is found 

— RR 

in Section 4. 3. The determination of the second column of the T matrix is considered 
a precalibration activity. In Section 4. 4. 3 that activity is described. 

4. 2. 4 Gyro Nonlinearity and J Term Experiments 

There are three remaining constants to be described, namely Qjj, Qj g and J. The 
Qjj constant is intended to represent the scale factor nonlinearity. In the following 
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CHART 4-19 
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GYRO CALIBRATION FIXED ORIENTATION POSITIONS 7 AND 8 
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GYRO CALIBRATION FIXED ORIENTATION POSITIONS 9 AMD 10 
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discussions we will describe not only how that constant is found, but also any other higher 
order w sensitive coefficients that are evidenced in the gyro readout. The Qjg calibration 
is not presented. It is assumed that this term is too small to be detected. The last 
paragraphs in this subsection discuss the calibration of the J term. It will be seen that 
this calibration requires test table speed controls not found in any other calibration. 

Gyro Nonlinearity 

The Qjj- term in the gyro equation was described in Section 2 as the nonlinearity term. 

This term is intended to describe, in conjunction with the scale factor term, the output 
rate (say P) as a function of the input (say w ) as: 

P = A + BW + CM 2 

rather than the more familiar 

P = A + Bin 

The interpretation by many is that this term introduces a nonconstant scale factor as- 
sumption, that is 

P = A + (B + C w) <o 

Regardless of the interpretation, it seems appropriate to assume nothing about the highest 
power of .os and in fact to try to conduct experiments to find all coefficients (say A J{ ) where 

• 2 n 

P = Aq + A ^ co + w = * * • A n io 

Such experiments are very simply described but would probably be somewhat time con- 
suming to implement. 

Let us direct our attention to the equations for the one gyro m Position 1, the two gyro m 
Position 3, and the three gyro m Position 5. We see that each equation can be written: 

P 2 
— An + A. to+ A<jto 
At u i ^ 

where Aq = R/ A# + f (acceleration) 

A 1 = * 1 / A $ 

A 2 = Qjj/ 
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If experiments are conducted where only tn is changed (and not acceleration), Aq would 
always be a constant. We see that a variation of the table speed only, in the positions 
mentioned, accomplishes this need. Also the equations can be generalized to contain higher 
order w terms. That is: 


P/ At = A 0 + A^to + A 2 T • * • A n w n 
The experiments can now be delineated: 


• Collect gyro data £6 from Position 1, 3, or 5 for n different speeds of the 
test table. 

• Collect the table speed data by measuring Av^ and At. 


u> =- 


A02 

At 


for constant to 1 and whole turns of 


• Let the total speed imposed on the gyro be described by 

^ ^2 F 

to = + to sin X 

At 


• Let P be given by 


P = 


(Sfi) 

At 


that is, data is collected sufficiently long such that quantization and noise are 
negligible (see Section 2. 2 of the Trade-Off Document). 

• Plot P against to 



• Analyze the plot to find f, where P = f (to). 


J Term 

The environments chosen for the determination of the J terms are shown in Chart 4-26. 
Note that the positions chosen correspond to Positions 7 and 11 (which were used in 
Section 4. 2. 2 ) with the modification of a rotating table. The gyro data will be collected 
over a period durmg which the angular speed has changed. This calibration is the only 
one which requires the determination of the time-varying integrals: 
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chart 4-26 


GYRO CALIBRATION ROTATING TABLE POSITIONS 7 AND 11 J 
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where y k = ul E sm X J'^ sm + T^j[ N cos tyit k - 1, 2 or 3 
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• J sin^ 2 dt 

• J cos ^ 


The manner in which these integrals are determined will depend completely upon the man- 
ner in which the test table angular acceleration is commanded during the experiment. The 
calibration procedures will therefore be dictated by the commanded angular acceleration 
profile. There appear to be only three interesting alternatives: 


• The first alternative, which appears to be the best, is when the angular 
acceleration can be controlled to a desired function of time In that event 
the aforementioned integrals could be evaluated prior to the data collection 

A good example of a commanded profile might be a constant angular acceleration 
over a short data collection time. 

• Another alternative, almost as good as the above, would be an angular ac- 
celeration profile which is an analytic function, but not known until the time 
of the experiment. An example would be the ability to command a constant 
angular acceleration, but not any given constant. In this event the integrals 
would be evaluated after data collection. 

• The least attractive alternative would be when the profile cannot be commanded 
as a clean analytic function. The J term experiment could be conducted under 
such circumstances, but there would be a requirement for the 02 resolver 

to be collected in real time for the purpose of evaluating the integrals. 


There is no reason to specify which of the above alternatives is to be used until the test 
table is evaluated to discern its ability to control angular accelerations. As a consequence, 
the J term equations in Section 4. 3. 2 are not specified as the equations to be programmed, 
as are the other calibration constant equations. Instead they are presented as functions 
of terms which will be described as functions of the angular acceleration profile at that 
time when the control characteristics of the table are better known. As a matter of con- 
venience to the reader the form of the equation is presented for the case when a constant 
angular acceleration profile is commanded. 


4. 2. 5 Determination of Accelerometer Coefficients 


The general accelerometer calibration equations were developed in Section 4 1. 2, and 
the results presented on Chart 4-13. We recall that it was assumed in the development 
of those equations that the test table would always be stationary (relative to the laboratory) 
during the entire accelerometer calibration. We recall also that nine of the gyro positions 
were also stationary. Analysis shows that the nine stationary gyro positions are very 
good choices for the entire accelerometer calibration. 
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We mentioned in Section 4. 2. 2 that, subsequent to the calibration of the principal angular 
velocity-sensitive coefficients, the gyro acceleration-sensitive terms predominate as un- 
knowns. Thus, when concerned with calibration only, we can treat the gyro as an acceler- 
ometer. We note that the three unbalance terms in the gyro equation play the same role, 
functionally, as the three scale factor and (Q^) terms in the accelerometer equation. 

We also note that the square term in the accelerometer equation appears functionally 
the same as the square compliance terms in the gyro equation. These facts, and the fact 
that the gyro and accelerometer input axes are nominally aligned, results in the use of the 
same positions for the determination of acceleration-sensitive coefficients in both the gyro 
and accelerometer equations. 

In Chart 4-27 we see the accelerometer equations when dj = 0° and ©2 = 90°. As with 
the gyro calibration discussed in Section 4. 2. 2, the accelerometer calibration requires 
only two table angle degrees of freedom. Therefore, two of the four table degrees of 
freedom can be arbitrarily chosen. The particular values of ® ^ and dg shown in Chart 4-27 
are chosen for the same reasons mentioned in Section 4. 2. 2 (where the gyro bias, unbal- 
ance. and square compliance calibration is described). In Charts 4-28 and 4-29 we present 
the accelerometer equations for Positions 7 through 12. We note that all but the cubic 
term can be explicitly extracted from these equations. (The cubic term always appears m 
any equation with the scale factor term and therefore cannot be separated from the scale 
factor term. ) In Section 4. 3. 2 the explicit solution for the accelerometer bias, square 
coefficient and off-diagonal (Q^)"^ matrix elements are presented. Three additional sets 
of equations are presented which relate the scale factor and cubic term combination to the 
instrument and environment measurements. 

4. 2. 6 Determination of Accelerometer Cubic Term 

It was noted in Section 4. 2. 5 that the six positions (7-12) did not allow for the explicit 
evaluation of the cubic or scale factor terms. We therefore require additional positions 
for the extraction of the cubic terms Positions 13, 14, and 15 (described in the calibration 
of the gyro product compliance coefficients) are appropriate as the additional positions. 

In Charts 4-30 and 4-31 we present the accelerometer equations for those positions. The 
equations for the solution of the cubic terms are presented in Section 4-3. 

For each accelerometer, either of the two positions in which the corresponding body axis 
is nominally oriented 45° off the vertical may be used. Therefore any two of the three 
positions may be chosen to complete the calibration. The equations presented in Section 
4. 3 utilize Positions 13 and 14. 
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CHART 4-27 
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CHART 4-28 



ACCELEROMETER CALIBRATION POSITIONS 7 AND 8 
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CHART 4-29 


ACCELEROMETER CALIBRATION POSITIONS 11 AND 12 
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CHART 4-30 



ACCELEROMETER CALIBRATION POSITION 13 
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CHART 4-31 


ACCELEROME TER CALIBRATION POSITION 15 
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4.3 CALIBRATION EQUATIONS 


We noted in Section 4. 2 that the simultaneous solution of at most two equations would yield 
the value of a required calibration constant as a function of the input environment and the 
defined or vectors. The and P-j^ vectors were seen in Section 4. 1 to be a func- 
tion of the instrument outputs, noise, and quantization error. In the first subsection which 

p A 

follows we will approximate the P^f and P k vectors by the instrument outputs only. Those 
approximations will be incorporated into hie determination of the calibration computations, 
which are tabulated in Section 4.3.2, 


4. 3. 1 Processing 

We showed in Section 4. 2 that each calibration constant can be solved for as a function of 

A G 

at most two input environments, two P-^ or two P^T components, and previously determined 
constants. That is, the solution for any constant (say y) can be written as: 


where 


Ax + B 


x = — 


A = 



environment coefficient of y 


P = P^ or P^ for k = 1, 2, or 3 


B = function of other calibration constants and environment 
inputs. 


In Section 4. 1 we defined the P vectors as a function of instrument readout, quantization 
error, and noise. We suggested in the introduction to Section 4 that we wish to approxi- 
mate the P vectors as functions of instrument readouts only. We would like therefore to 
collect the instrument data in such a way that the effects of quantization and noise fall be- 
low some required threshold. Fortunately the noise can be represented by random pro- 
cesses with bounded means and variances. The quantization error is by its very nature 
also bounded. On the other hand, for a nominally constant input environment, the instru- 
ment output is a monotone increasing function of the observation time. Thus by choosing 
sufficiently long observation intervals, the percent error in the assumption that the instru- 
ment output equals P can be made arbitrarily small. 
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In Section 2 of the trade-off document the analysis which leads to the above conclusions is 
presented. The results of that analysis are presented in form of graphs in which the 
precisions of the calibration constants (with the assumption of neglected noise and quanti- 
zation in or ) are plotted against time. Those graphs will be used in the Laboratory 
Procedures Manual to determine the calibration time required to obtain a desired precision 
in any constant. 

4. 3. 2 Computation of Constants 

With the approximations described in the previous section, it is now possible to solve 
explicitly for the calibration constants m terms of well-known instrument and environment 
measurements. The equations are presented in tabular form on the Calibration Equations 
Charts. These equations are those which are to be programmed. An exception is the 
J term equations which will not be in program form until the time when the test table is 
evaluated (see Section 4.2.4). 

It has been noted in Section 4.2.5 that the accelerometer third order term cannot be 
separated from the scale factor by a choice of positions. In the following set of equations 
there are two equations given for each accelerometer that relate the scale factor term 
[D^(A.- B )] to the third order term (DjDg)j. If a simultaneous solution of the two equa- 
tions is used to determine the scale factor and the third order term, then the scale factor 
will be sensitive to errors in the bias and the second order term. These terms appear on 
the equation listed second in each of the three sets of two equations. This may be avoided 
by determining (D^Dg^by simultaneous solution or other methods) and using this value to 
solve the first equation for [D j . This value is subject to the accuracy of other 

terms only through the extremely small term containing Dg. (D^) is then given by the square 
root of the sum of the squares of [DdA «B )j for j = 1, 2, 3. 

-LI j 1 


In developing these equations several equalities are used which introduce previously un- 
mentioned parameters. The following comments describe those parameters and their 
nomenclature: 

• Because we are dealing with the measurement of quantities m many different 
positions, a position index must be introduced. The numerical superscripts in 
all equations refer to the positions 

• The vector (k = 1, 2, 3) is approximated by the gyro readouts, therefore 

p k - < E6 >k 
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GYRO CALIBRATION EQUATIONS 
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GYRO CALIBRATION EQUATIONS (Continued) 


Gyro Two (Continued) 
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GYRO CALIBRATION EQUATIONS (Continued) 
Gyro Three (Continued) 
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GYRO CALIBRATION EQUATIONS (Continued) 
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GYRO CALIBRATION EQUATIONS (Continued) 
Gyro Three 
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GYRO CALIBRATION EQUATIONS (Continued) 
Gyro Two 
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GYRO CALIBRATION EQUATIONS (Continued) 


Product Compliance 


EC Ig /A®] t — 


(S6){ 3 


^ Kf 13 J 


Gyro One 
2 V2“ 


VZ r 1 

[R/A^ +— lBj/ASL - CB g /A$] 

£ L 


2(J P r 

uGj^. B 1 /A<a>] 1 [VI72 cos A + T RR13 sin A] + [(G^ B^/AS^ 

g 2 

[ViTz cos A + T® R13 sin a] + [Gj. Bg/A*^ [T RR13 sin a] 


+ lCjj/ A®3 1 + [Cgg/A#^ 


ICjo/A^Ii = + 


P 6 >i 4 1 2 


g 2 S* 1(^)0141 g2 


— [R/A#^ [<• Bj/A$j x + CBq/A^Ij] 

or^ S 


lC 0S^1 - " 


• G x - Bj/AC]^ [yiTFcos A + T RR14 sin a] + CG^ Bg/A#^ 

g 

[ T BRl 4 sm x j + -Q^. Bg/A*^ [VTTTcos A + T® 2 R14 sm a] - [C n /A$J 1 

2 (^i 5 2 _ V2~ r 

+ _ lr/mA + — -LBq/A + [B n /A?>] 1 

«A,T ^„TxG15 J2 1 „ L S 1 O' lj 


2 

(S6)J 5 


^ TxG15 
[F n i } i . 

0 E 

2 W r 



(\/I71T cos A + T 22 R15 sin a) + LGy Bg/A^ jji/2 cos A + T RR15 smA]J 


+ lC ss /A?>^ 


4-55 



GYRO CALIBRATION EQUATIONS (Continued) 
Gyro Two 
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GYRO CALIBRATION EQUATIONS (Continued) 
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GYRO CALIBRATION EQUATIONS (Continued) 
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GYRO CALIBRATION EQUATIONS (Continued) 
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ACCELEROMETER CALIBRATION EQUATIONS 
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ACCELEROMETER CALIBRATION EQUATIONS (Continued) 


Bias and Second Order Term 
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ACCELEROMETER CALIBRATION EQUATIONS (Continued) 
Accelerometer Three 
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ACCELEROMETER CALIBRATION EQUATIONS (Continued) 
Accelerometer Three 
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where ( 7 ' 6 ) K corresponds to the pulse count from the kth gyro over the total 
time of data collection. 

• The vector p 4 (k = 1, 2, 3 ) is approximated by the accelerometer readouts, 
therefore 


pA => . Py)kl 

where (Sy)^ (£y)kl correspond to the pulse counts from the second and 
first string of the kth accelerometer, respectively. 

• The total time of data collecting will be recorded as a count ( s nt) from the 
system clock, therefore 

H =*> ( s nj) 

T 

where S , is the scale factor of the system clock. The subscript 1 serves to 
distinguish the clock pulse train used in calibration from the pulse tram used in 
the Preprocessing computations. 

• The total test table angle (A^) will be recorded as a number of whole turns 
(or a number of fractions of whole turns), therefore 

Adg =0 S 0 (^n^) 

where (^n 0 ) is the number of increments of angular displacement, and is 
the scale factor which converts the number of increments to a finite angle. 


In Section 3. 3 we presented a description of the laboratory facility with all of its measure- 
ment and computational devices. In Figure 3-6 we showed the possible equipment interfaces. 
For the purpose of calibration all instrument data collections will be accomplished with 
the frequency counters shown in Figure 3-6. (See Section 2. 1 of the trade-off document 
for the reason why the counters are used. ) The specific employment of the counters for 
all positions is described in the Laboratory Procedures Manual in the sections entitled 
Fundamental Modes. Also found m the Fundamental Modes sections are all events in the 
collections, transfers, and computations during calibration in the form of flow diagrams 
accompanied by descriptions of the activities 

4. 4 PRECALIBRATION REQUIREMENTS 


The required constants contained within the equations tabulated in Section 4 3 were pre- 
sented as functions of instrument outputs and parameters describing the environment inputs. 
Before the data can be collected which is necessary as inputs into the equations, several 
initial survey tasks must be accomplished. 
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• The environment selections were presented in Section 4. 2 as choices of the 
and 0 2 gimbal angles and the first column of the T BB matrix The test table 
orientation is, however, controlled bv choices of four gimbal angles; therefore, 
the choices of the first column of T BB - must be expressed in terms of the gimbal 
angles, which we know from Section 3 to be 0^ and 0^. 

• The T BR matrix is a function of the T B * matrix as well as the 03 and 04 gimbal 
angles. We must, therefore, determine T B * before we can equate the first 
column of T bR to 0 g and 04 . 

• The gyro bias and compliance constants are seen (see Section 4. 3. 2) to be a 
function of (among other things) the second column of the T B & matrix. That 
column can be determined (once T B I. 0 3 , and 0 4 are known) by a use of equalities 
presented in Section 3. 

TP 

• In all previous developments it was assumed that T (the transformation between 
the test table base frame and the earth frame) was equal to the identity matrix. 

In the operational laboratory this matrix will deviate (by small numbers) from the 
identity matrix. It is, however, possible to correct for the deviation by the use 
of bubble levels. 


In the following subsections we present the manner m which all of the above tasks are ac- 
complished The order of presentation is the chronological order in which these tasks 
should be accomplished in the laboratory. 

4. 4. 1 T BI Survey 

The initial activity subsequent to the attachment of the ISU to the test table is the determi- 
nation of the orientation of the ISU body axes relative to the test table frames. This cor- 
responds to the evaluation of the T matrix (see Section 3.2. 1). In Chart 4-32 we see how 
this is accomplished. 

In Chart 4-32 we refer to the test table orientation used in the determination of T BI as 
Position Zero. Position Zero can be any orientation, but the zero orientation shown in 
Figure 3-2 might be the most convenient for it results in 

( T TO T OR T RT T TF T FE ) _1 = I (the identity matrix) 

BS SE BE 

(The matrix product T T = T is functionally equal to the Mirror Alignment Matrix 
shown in Chart 2-6. As mentioned in the discussion of the Mirror Alignment Matrix is 
Section 3. 1.2 the particular evaluation of this matrix depends upon the particular geometric 
angles which are outputted from the autocollimators. At the time that the form of those 
outputs are known, the exact form of T B ® T^ F can easily be determined, ) 

HI 

The T Survey activity is formalized as a "Pre calibration" activity in the Laboratory 
Procedures Manual. It will most probably be accomplished very near the time that the 
ISU is placed on the test table. It will probably not be necessary to repeat this survey 
except when the ISU is removed from the table and then "rebolted" in another orientation. 
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CHART 4-32 


T Determination 

Given, a set of resolver readings 0 ° 0 ° 0 ®, and <b\ from position zero and the 
ta nc S n for ? ation from body to earth axes via autocollimator readings 

/pSEOj 

Find: the matrix T BI which transforms from the body axes to the axes fixed 

to the inner gimbal of the test table. 

1. From the laboratory geometry definitions described in Section 3 we have: 

T BI = ( t BS0 t SE 0) (^ T rOO T ORO T RTO T TFO T FEj-l 
BSO SEO 

where T T is given oy the autocollimators and 
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4. 4. 2 Test Table Resolver Settings 


The calibration selections in Section 4. 2 were accomplished by dictating values for 0 -., 

0n and the first column of the T matrix. In implementing those selections it is neces- 
" BR 

sary to equate the T choices to the controllable 0o and 0 4 gimbal angles. In Chart 4-33 

^ ^ BR BI 

we see the functional relationship between the first column of T and 0 3 and cp^ (and T 
which is known from the procedures developed m Section 4. 4. 1). 

In Chart 4-34 we present the 0 4 , 0 g, 0 3 , and 0 4 settings for all fifteen calibration positions. 

Included in that chart are the equations for the determination of 0 Q and 0, from the choices 

BR ^ 

of the first column of T . Those equations are special cases of the arithmetic contained 

in Chart 4-33 for Positions 1, 3, 5, 13, 14, and 15. A duplicate of Chart 4-34 is presented 
in the Laboratory Procedures Manual. The numerical solutions for 0 3 and 0 4 must be ac- 
complished and placed in the chart before calibration can be accomplished. 

4. 4. 3 T BBm Determination 

The gyro bias and compliance computations presented in Section 4. 3 are functions of 
(among other things) the second column of the T matrix for Positions 1, 3, 5, 13, 14, 

■DTD 

and 15. The T matrix is a function of $ 3 and 0 4 (see Section 3. 2). The gimbal angles 

0o and 0. are known by the use of the computations presented in Chart 4-33. In Chart 

^ ^ BR 

4-35 we present the computations which develop the required second column of T from 

the known 0 3 and 0^ andgles (and T B ^ as given by the computations presented m Chart 4-32). 
Although not required for calibration purposes, the computations of the third column of 

DD 

the T matrix are presented for information purposes in Chart 4-35 The first column 
BR 

of t is known because it was utilized as the environment selection control parameter. 
4.4.4 Bubble Level Corrections 

Throughout the development of the calibration techniques it was assumed that the test table 
base frame is aligned with the earth axes, that is: 

T = I (the identity matrix). 

In practice this matrix will deviate from identity, due to such things as solar heating of 
the building and settling of the building. The resultant low frequency motion of the base 
relative to the earth can be corrected (immediately before calibration data collection) by 
the use of buble levels. There are three ways in which the corrections can be implemented. 
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CHART 4-33 


0™ 0™ Equations 


rivpn t BI ^BRm rpBRm T BRm 
» i ll » i 21 » 31 

Find: 0 ™ and 0 ™ 


1. We know that 


,pBRm __ ^,BI rpIOm r^ORm 


or (rpBI)TrpBRm_ T IOm rpORm 


■RT 

where T is given from a prior survey and 
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sin 0 “cos 0™ 
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COS 0 ^COS 0“ 





_ 


2. Solving for the first column of T IOm T ORm we have 



which gives the desired functional relationships. 
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CHART 4-34 


Position 1 


Position 7 


TEST TABLE RESOLVER SETTINGS 
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CHART 4-35 


rpBRm Det ermina t 10n 


Given: 0 “, 0 “, and T BI 

Find: T BRm 

1. We know that 

rpBRm _ r pBI t pIOm r j,ORm 

DJ 

where T is given from a prior survey and 


rpIOm^ORm 
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2. The first column of T BRm is the calibration control parameter. The values of 
this column have already been included m the calibration equations. 

O LJ yyn 

3. The second column of T which is the only column required for inclusion into 
the calibration equations, is: 

rpBRm _ rpBI ,m rpBI . ,m 

T 12 " T 12 cos 0 4 " T 13 sin 0 4 


rpBRm _ rpBI _,m ^BI A m 
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T 32 “ T 32 cos 0 4 " T 33 sm 0 4 


These equations need only be solved 
for m = 1, 3, 5, 13, 14, and 15. 


DO TV! 

4. The third column of T , which is not required in the calibration equations, is 
given for information: 
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• The first and most obvious way is to incorporate 'leveling screws” on the base 
of the table. 

BI 

• Another technique is to store the T matrix in the computer as other than an 
identity matrix, and to input the bubble level readings taken before data collection. 

• The third technique is to place two of the table gimbal angles in the direction of 
the two bubble level degrees of freedom, and correct for the lean of the table by 
small corrections to the gimbals. 

We choose die third technique because it is the easiest to implement. 

The bubble level corrections for the rotating ISU in Positions 1 through 6 need only be 
imposed on the trunnion axis (i. e. 0j). As a matter of fact, if the motion of the table 
base is less than, say, 20 s'econds no correction is required. This is because the prin- 
cipal input to the gyros in Positions 1 through 6 is the table speed, which is unaffected 
by the motion of the base. In Positions 7 through 15, however, the principal input to the 
instruments is g; therefore the level corrections are essential. This is why, in every 
position after the sixth, we placed the outer axis in the north- south direction. Assuming 
that the two bubble levels are in the east-west and north-south directions, level corrections 
can be imposed on the 0 -^and 0 3 gimbals. 

4. 5 IMPLEMENTATION OF CALIBRATION TECHNIQUES 

In this subsection we begin the formulation of all of the aforementioned topics into a for- 
malized laboratory calibration procedure. In formulating these procedures we will intro- 
duce several housekeeping operations which have not been previously mentioned. 

The discussions in this subsection are directed toward a presentation of the chronology 
of events which occur during the calibration. We will not go into great detail, the detailed 
presentation being left as the task of the Laboratory Procedures Manual. This subsection 
should be considered as the interface between this Development Document and the Labora- 
tory Procedures Manual. 

The calibration procedures are divided into four separate activities. 

• Turn-on 

• Precalibration 

• Calibration 

• Computation. 
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Whether accomplishing a complete or partial calibration, the above four activities occur 
in the order presented. In the following paragraphs we will briefly describe each activity: 

Tum-On 


These activities include all of the various housekeeping tasks which include such things as: 

• Power-on to equipment 

• Monitor equipment operation % 

The detailed specifications of these activities cannot be tabulated until the laboratory facility 
is completely defined, hi the Laboratory Procedures Manual a space has been allocated 
in Part I for inclusion of the details of turn-on to be specified at the time when the labora- 
tory is configured. 

Precalibration 


At some time, between the placement of the ISU on the test table and the initiation of cali- 
bration , the following system survey activities must be accomplished: 

• Determine T®^ 

• Find 0g and 0^ for all calibration positions 

• Find the second column of the TBR matrix for all calibration positions 

• Store g, uP, X, S T , and S.^ 

The first activity locates the ISU relative to the test table. The second activity determines 
the inner and outer test table gimbal angles settings for all calibration positions. These 
settings were shown in Section 4. 4. 2 to be a function of the T BI matrix. The third pre- 

T5TD 

calibration activity computes the second column of the T matrix, which was shown in 
Section 4. 3. 2 to be necessary for computing the gyro bias and compliance terms. The 
fourth activity records system numbers required in the calibration equations. All of these 
activities are described m detail in Part I of the Laboratory Procedures Manual. 

Calibration 

At any time subsequent to the completion of the Turn-on and Precalibration activities, the 
ISU can be calibrated. We formally define the calibration activities as the completion of 
the following list of activities for any or all positions: 
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• Connect instruments to frequency counters 

• Set frequency counters 

• Set test table resolvers with bubble level corrections 

• Set table speed 

• Collect data 

• Transfer data to computer. 

A complete calibration would accomplish the above activities for all positions The order 
of positions is completely arbitrary; but the packages of six (1 to 6 and 7 to 12) will proba- 
bly be accomplished m numerical order. A partial calibration need only accomplish these 
tasks for positions required for the determination of the required constants. The details 
of the above activities are found in Part II of the Laboratory Procedures Manual. 

Computation 


Computation is very simply the solution of any or all of the equations in Section 4. 3. 2. A 
complete calibration requires all computations and a partial calibration would require the 
solution of only a few of the equations in Section 4. 3. 2. The details of the computation 
procedures are described in Section II-4 of the Laboratory Procedures Manual. 
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SECTION 5' 



SECTION 5 

DEVELOPMENT OF ALIGNMENT TECHNIQUES 


In Section 2 alignment was defined as the initialization of the matrix which transforms from 
an ISU-fixed set of axes to a navigational set of axes. The ISU axes were defined by use of 
two mirror normals and the navigation axes were defined as an earth-fixed, local-level 
frame. Transformation of this alignment problem to any other alignment problem using 
different ISU and navigational frames is then a simple problem of coordinate transformation. 
The discussion in Section 2 revealed that alignment could be accomplished by measurement, 
in body and/or earth-frame, of the components of two system vectors. Three different 
choices of these vectors lead to the three alignment techniques: Mirror Alignment, Level 
Alignment, and Gyrocompass. 

Further analyses in Sections 2 and 3 lead to the functional description of the three alignment 
techniques shown in Chart 5-1. In this description each technique is further broken down 
into four basic types of computational routines. These are: 

o Preprocessing Computations 

® Estimation Routine 

• Estimation Matrix Computations 

• Alignment Matrix Computations. 

These routines have as inputs certain a priori information, calibration constants, instru- 
ment outputs, and/or outputs from other routines as indicated in Chart 5-1. 

Before beginning the detailed development of alignment it is important to note several points 
which dictate the viewpoint adopted in the remainder of this section. First note that there 
are basically three types of routines indicated in Chart 5-1. They are Preprocessing, 
Estimation (including both Estimation Routine and Estimation Matrix Computations), and 
Alignment Matrix routines. The mathematics of the Preprocessing Routine was developed 
in subsection 2.2.5 and will be considered only briefly here (m Section 5. 1). The Alignment 
Matrix routine uses estimated values of g‘B, and mirror azimuth (Level Alignment) or 
g • R and to (Gyrocompass) to initialize the alignment' matrix T. This relatively 

straightforward mathematical problem has been discussed in Section 2 and is considered 
again rather briefly in Section 5. 5. The remaining routine, Estimation, is the major sub- 
ject of discussion in this section. Before developing various estimation techniques, the 
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CHART 5-1 
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environment and sensor noises are modelled m Section 5. 2. Then in Section 5. 3 and 5. 4 
specific estimation techniques are developed under different assumptions. The explicit 
equations for the recommended alignment techniques are summarized m Section 5. 6. 

Second, it should be noted that there is no detailed discussion of Mirror Alignment since 
this section emphasizes estimation which is not relevant to the mirror alignment problem. 
The Alignment Matrix Calculation discussion of Section 5. 5 is, of course, applicable to the 
mirror approach when mirror azimuth and zenith angles are given. 

Third, the reader should be forewarned of the emphasis on Level Alignment over Gyro- 
compass in this section. It was found, not unexpectedly, that the alignment errors m 
Gyrocompass may easily be two orders of magnitude larger than those expected in Level 
Alignment. This is, of course, mainly due to the low signal-to -noise ratio of the earth- 
rate signal in gyro quantization noise. Further details of this comparison of Level Align- 
ment versus Gyrocompass are given in Section 5 of the trade-off document (Volume 2) 
where a Monte Carlo simulation of an alignment problem is used to obtain quantitative 
results. 

Finally note that Section 5 of the trade-off document justifies many of the comments included 
below. Section 5 of Volume 2 includes further discussion of the assumptions required and 
the results of a simulation of the proposed estimation techniques. It is important to note 
that in several places important assumptions have been made with little justification when 
the data was not available to include completely realistic values. The collection of ac- 
curate data about the noise environment is a very difficult and expensive problem. How- 
ever, all results of this study have been presented m such a manner that when more ac- 
curate data is available, modifications can easily be made. 

5. 1 PREPROCESSING COMPUTATIONS 

The Preprocessing Computations yield integrals of angular velocity and acceleration in the 
body frame • 

t+At t+At 

J w-B dt | a-B^dt , j - 1, 2, 3. 

t t 
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The inputs to the computations are the counter outputs (Sy)^, (Sy) kl , (£ 6 ) k , and (S n J), 
k = 1,2, 3. The calibration constants for each sensor are required, along with the trans- 

A P 

formations from the sensor axes to body axes Q , Q . The notation is defined in Section 
2. The preprocessing equations are developed in subsection 2.2.5. They are reproduced 
in Chart 5-2. In this chart, dots at the left hand margin indicate the alternative computa- 
tions used for Level Alignment or for Gyrocompass. Note that At was assumed to be 
small in the development of the equations in Section 2.2.5. This restriction is not required 
in the Preprocessing Computations used in alignment since the ISU is relatively stationary. 

5.2 ENVIRONMENT AND SENSOR NOISE MODELS 

Before developing processing techniques, we must describe the effect on sensor outputs of 
various random inputs: environment translational acceleration and rotation, accelerometer 
noise, and gyro noise. This section is a continuation of Section 3.2,2 which describes the 
general characteristics of the environment noise. In the development of estimation tech- 
niques, quantization errors are not included. Several of the resulting techniques are tested 
with a Monte Carlo simulation to determine the effect of quantization, computer word length, 
and anomalous noise inputs. The results are presented in Section 5 of the trade-off docu- 
ment. In the following paragraphs, we first describe the environmental components of the 
sensor inputs (5.2. 1) and subsequently describe their effects on the observed sensor out- 
puts (5.2.2). 

5.2. 1 Sensor Input Acceleration and Angular Velocity 

It is convenient to define a "level frame" that moves with the ISU and whose average orien- 
tation is collinear with the earth axes, as indicated in Figure 5-1. The body axes are fixed 
relative to the level axes. If there were no environment disturbance, the level frame would 
coincide with the earth frame. 

Let a^ be the acceleration of the level frame along L^. Then 



a Ll 


gu- h 1 

+ a l 

-L ~ 

a L2 

= 


+ a 2 


a L3 



+ a 3 
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PREPROCESSING COMPUTATIONS 


Inputs (Ey) k 2 > (Sy) k i, (E 8 ) k ,and (llnj)for k = 1 , 2 , 3 

The outputs 'B k )dt and J’ t t+ {a-B k )dt (k = 1 , 2 , 3 ) are given by the following 

computations 


Level GC 


if - [(sr) k2 - fcy) kl ] 
p k E 


• L(a:.G k )zit] £ P k (A4>) k - (R) k At 

• [(^A k )At] - P k /(Di) k - (D Q ) k At 

° = C(w-G k ) At J/ At 

• (a- A k ) s [ (a • A k ) At]/ At 

• (a*G k ) = (a- A k ) 

0 0 i"j Ra- 

0 (aRp =10 0 (aTAp 

ioo (|Ta^) 

"o -1 0] RaTAp" 

• (a^) =00-1 (a^) 

0 10 (aTE^) 

.Jt t+At (u.-Gk)dt = :7i-G k )Atj-RB I ) k (i^.) + (B 0 ) k (i^K(B s ) k g^)lAt 

- t( c n ) k ( a.G k ) 2 + ( c ss ) k (a-S k ) 2 - At 

- ^ (G]g) k (a* G k )(a. S k ) + (GQg) k (a* O k )(a- S k ) (GjQ) k («i‘G k )(a*O k )^At 

- [(Q n ) k (^) 2 + (Q IS ) k (iGj;)(^):At 
vf Al R-A k )dt = [(a*A k )At: - (D 2 ) k (i 7 A^) 2 At - (D 3 ) k <£TA^j 3 At 

•J t t+A V*§ k ) dt = 

•J^VSkJdt-SQtiJ^Atfe-AiWt 


-(G r B 2 ) 


Q = -(Gg-B^ 


-(g 2 'B 3 ) 


-%§l) -%§2> 


(Aa'Si) 


-^•Bg) -(A^) 

1 Ra 2 - b 3 ) 


-(A3- B^) -(A3- Bg) 
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Figure 5-1. Level Frame 


where j, «g, a g are the environment-induced accelerations along L-p Lg, Lg. The 
expected value of is zero. To first order, • U = 1 since the rotation from vertical 
is the order of one minute of arc or less. Further, U • Lg and U • Lg can be represented 
as small rotations about Lg and Lg, respectively. Let 6g = -U-Lg and 9g = U*Lg. 
Then, to first order, 


g + Otj 
- g 0 3 + « 2 
g0 2 + a 3 


Further, the angles 9g and 6g have zero expectations. Let T ^ be the orthogonal trans- 
formation from the level frame to the body frame. Then the acceleration in the body frame 
a B is 


g 


0 


v 

0 

+ T 1 

-ge 3 

+ T i 

a 2 

0 


_ ge 2_ 


*3 

— 


The first term is the average gravitational acceleration m the body frame. The second 
term is the variation of the acceleration due to level frame rotation. The last term is 
the environmental acceleration disturbance. The power spectra of a p aig, ag, 9g, Qg 
can be obtained from environmental test measurements. 

Next, consider the angular velocity. Let be the angular velocity of the level frame 
about L^. Then 



“’Ll 


lugu.Li 

+ 

“N — *-l 

+ 

'I " 

— L = 

40 L2 

= 


+ 


+ 

^2 


^L3 


"tk'i's 

+ 

^?-ls 

+ 

^3 


where is the angular velocity of the earth and fiy. is the environment- 

induced angular velocity about L^. To first order, U*L^ = N'Lg = Further assume 
that there is no rotation about U and hence P ^ = 0 and N • Lg = 0. Also, 0g = 9g, P% ~ 0 3 > 
and N*L^ = 9g. Hence > 
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where Wg is the angular velocity on the body frame. These three terms have interpreta- 
tions analogous to the corresponding terms in a^. 

Since the system has been calibrated, the transformation from body axes to sensor input 
axes is known. Let Tg transform a B to the accelerometer input axes; let Tg transform 
aJ B to the gyro input axes. Let 


= T 2-B 


and 


-G = T 3— B 


These vectors represent the sensor inputs. 


There are two alternatives in estimating gravity and angular velocity. First, we can 
estimate the average components in the body frame: 


— — i 


E 

s 


“u 

0 

and T x 

0 

0 


E 
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Second, we can estimate the components at some tune t* in the future, t* > KAt: 


g 


0 

0 

+ T 1 

-gfl 3 (t*) 

0 


g 0 2 (**) 


(5-2) 
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(5-3) 
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T 1 

0 

+ T 1 
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L “N ... 


io®e 2 (t») 


Where the strapdown system is initiated t* - KAt seconds after the last measurement. 
Let t* - KAt = , the prediction interval. 

5.2.2 Observed Sensor Output 

In the development of a processing method, we will assume that the gyros and accelerom- 
eters are linear with unit scale factors and zero bias. Actual values will be used in the 
application to a real alignment problem. Further, we assume that the gyros and accel- 
erometers have relatively large band widths, i.e. , we will neglect the sensor dynamics. 
The preprocessed sensor outputs are integrals of acceleration and angular velocity. 
Namely, the outputs are: 


A t t 

?. W = J a A (r)dr + J n A (r)dT 
t- At t-At 

P G (t) = J c^(T)dT + J* n G (T)dr 
t- At t- At 

where n A (r) and n g (t ) represent noises introduced by the sensors. Further the outputs 
are observed at discrete times At, 2 At, — , KAt. Denote these outputs by P A (i) and 
P (j). It is convenient to transform the outputs P (j) and P y) to the body frame: namely 


?b(J) = 


J At 

s 


"0 

) 

X h 1 

0 

+ T 1 

ge 3 ( T ) 

( dr 

(]-l)At ( 

0 


( 


r ]M \ 



«j(t) 

■f 

n Al (r) 

T 1 

o^(r) 

+ T 2 

n A2^ r ) 


_a 3 (r)_ 


n A3< T > 


dr 


(5-4) 
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?B<J> 


= T. 


^p^]) = J jAt )t x 

(j-l) ' 
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+ T i 


-oj 
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0 9 (r) 


N u 2' 

E 

w u 8 2< r ) 


0 3 (r) 


r jAt 
+ TiJ 

(3-1) At 


L dr 
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) 

0 

♦ 

3 At 1 

n G l( T ) 


e 2 (r) 

dr + | Tj 1 

n G2 (r) 

dr 

a 3 (r) 

(3-1) At 

n G3 ( T ) 
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Since is orthogonal, and since Tg and Tg are nearly orthogonal, T x n(r) + 
and Tg n. G (r) have the same power spectrum as n(r) + n A (r) and n G (r), respectively. 
This simplification will be used in equations (5-4) and (5-5) since only second order 
statistics are used in the following discussion. 

Note that the components of od(r), n^(r ) , andn G (r) have been assumed statistically 
independent and identically distributed. 

In the following sections, estimation techniques are developed based on the above models. 
5. 3 ESTIMATION OF GRAVITY IN LEVEL ALIGNMENT 

Estimations of the components of gravitational acceleration m the body frame are based on 
the observation equation 


p£(5) 


J ]At 

(3-1) At 


i 

g 


0 

K 

0 

+ T 1 

-gs 3 (r) 

( 

0 




+ a(r) + n a (t) > dr 


(5-6) 


with j=l, ,K. Using the observed accelerometer outputs in the body frame, Pg(j) 

(]=1, . . . ,K), our goal is to estimate the average components (5-1) or instantaneous 
components (5-2) in the presence of the disturbances o(t) and n A (t) given in (5-6). The 
former problem is described in the following Section (5.3. 1); the latter problem is 
discussed in (5.3.2). Note that the basic operation is differentiation; we obtain ac- 
celeration from velocity measurements. 
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In the following subsections several estimation methods are developed from mathematical 
statistics. A general discussion of Level Alignment techniques appears m 5. 3. 3; and the 
recommended technique is presented in 5. 6. The characteristics of several techniques are 
described in Section 5 of the trade-off document. These techniques are simple average, 
posterior- mean estimate of average components, posterior- mean estimate of instantaneous 
components, and iterative estimate of instantaneous components. 

5. 3. 1 Estimation of Average Components 


The objective here is to estimate the average gravitational acceleration components in the 
body frame namely: 

Ib = T l 

Let gg(t*) denote the true gravitational components at time t*. Then the rms deviation of 
the estimated average components (say § B ) from the true components can be bounded, 
namely: 



E [(g B ( t!|t )"lB^ T (lB^^“^B^ " E i fefi ^ "Sb^ T _ §B^ ^ + E C ^Ib"^B^ T ^B~^B^ - 


+ 2E[(g B (t*)-f B ) T (g B - g B ) 1 ^ | E l(g B (t*)-g B ) T (g B (t*)-| B )] 


1/2 


E c^ B -i B ) T( lB“iB ) - 


1/2 J 2 


The first term corresponds to the rotational motion of the level frame about the average. 

The second term corresponds to the error mgg as an estimate of gg. The objective is 
to minimize the second term, accepting the first term (the error from the motion about 
its mean). 

In the following subsections three approaches to estimation are considered - simple average, 
least squares, posterior mean. The first approach does not use any a priori information 
about the noise spectra, alignment, gravity or earth rate magnitudes. The second approach 
uses prior measurements of the noise spectra. On the other hand, it does not include the 
prior geophysical measurements of gravity and earth angular velocity. The third approach 
uses a priori information about alignment, gravity magnitude, and earth rate magnitude 
plus measurements of the noise spectra, but the noises are assumed to be gaussian 
processes. This third approach has several advantages: (l) prior geophysical measure- 
ments are included and are weighted with their accuracy; (li) the estimation techniques are 
comparable to those obtained from a least squares approach in complexity; (m) the resulting 
techniques can be used recursively to continuously update the alignment matrix; (iv) the 
posterior- mean estimate is optimum with respect to a large class of loss functions, not 
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just quadratic. From noise simulation, we find that the posterior-mean techniques are 
not sensitive to the gaussian assumption. (See Section 5 of the trade-off document. ) 

5. 3. 1. 1 Simple Average 

This approach is based on the assumption that we do not have any prior information about 
the noise, gravity, earth rate, or alignment. In this case 




KAt ]=1 


Pg(i) 


Note that the same estimate is obtained if K=l, and At is replaced by KAt. 

5.3. 1.2 Least Squares 

Before developing a least-squares technique, it is convenient to define certain notations. 
Let A be the 3K vector whose components are 


x , - p bi«>; 


X j+K P B2^ X j+2K 


- P B3« 


with j = l, . . . , K. Let be the 3Kx3 matrix 


H 


1 


At 0 
♦ « 

• • 

At 6 

0 At 

» « 

9 O 

o a{ 
0 o 

9 • 

* * 

• 9 

0 0 



K-rows 

K-rows 

K-rows 


] At 

Let Vi) = (f.DAt 8r<T)dT > with r ° 2 ’ 3 - 


5-12 



Let $ be a 3K vector with components 


^3 + g ^ T 1^13 ^2^ " g ^ T Pl2 ^3^ 


$ 


]+K + g ^ T l^23 ~ ^ T 1^22 *3^ 


(5-7) 


$ 


]+2K - + g( Tl ) 33 0 2 ( 3 ) - g( Tl ) 32 # 8 (j) 3 = 

Let 0J = (^2(1), • • • ,^2®) an d ®J ~ Then f-Hg 


—2 


L-3J 


where H 2 is the 3Kx2K matrix defined by equation 5-7. Let N be a 3K vector with 
components 


At 


N “ J . . Ca^r) + n .(t)] dr 
1 (]-l)At 1 A1 


jAt 

N j + K - (j I 1)At ^r) + n A2 (.r)-iir 


JAt 


Vi “ J 1)At t“ 3 (r)+n A3 (T)3dT 


with 3 = 1) . . . , K. Then the basic observation equation 5-6 can be rewritten as 


X = Hj g B + £ + N 


(5-8) 


The objective is to estimate gg m the presence of noise ^tN, given the observations X. 
The covariance of this composite noise is the sum of the covariances of $ and N since 
they are independent: 

S $ + N = S< i> + 
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, Further 



where £ ^ is a K x K matrix with elements 

( s 00) i3 = E[0 2 (i) 0 2 (])1 , 1,] = 1, . . . , K. 

Note that ^ and 0g are assumed to be independent and identically distributed. Further, 



where + n is a K x K matrix with elements 


(E a + n )jj = ECN.Nj] , i, j = l,...,K. 

These covariance matrices can be expressed in terms of the correlation functions, 
namely 


E!> 2 (i) «> 2 (j)3 



At 0 

E1N.N1 = J [ At-p] (p + (j-i) At) dp + J [At + p] 0(p + (]-i) At) dp 
1 J 0 -At 

At 0 

+ J [ At-p] 0 (p +(]-i)A t) dp + J [ At - p] 0 (p + (]-i) At) dp 
0 n - At n 
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where 0 Q , 0 a ,3nd <t> n are the^ correlation functions of 0(r), g(t), n(r), e.g. , 0 0 = E(0(O)0 (t)). 
From prior experiments we estimate the noise power spectra and correlation functions. * 
The above integrals may be evaluated numerically or mathematically when the correlation 
function is approximated by a mathematical formula. For example, assume 

= c ie - c 2l r l 


2 Cl e-^lHl At 

= — [ cosh (c 0 At) -1], i ^ j 

c 2 

= ^ [ e ~ c 2 At -1 + Co At 3 

(^ 2 ) 2 

The same methods can be applied to the other covariance matrices. 

Note that H 2 is evaluated by using prior estimates for the value of and g, denoted by T^ 
and g. Precise values are not needed since Hg is used in the noise model. Corrections to 
Hg would be of second order. 

Based on the composite measurement equation (5-8), the objective is to find the unbiased 
linear estimate of g g , say gg(X), which minimizes E|gg - g(X)| 2 as a function of g(X). 

It follows from the Gauss-Markoff theorem that g-JX) is the value of g that minimizes 

9 — tD — rn 

| MX - M Hjg| , where M is the nonsmgular matrix such that = I. ** 


0 e (i) 

Then 


© ) . 




^Spectra data is given by H. Wemstock in "Limitations on Inertial Sensor Testing Produced 
by Test Platform Vibrations", NASA Electronics Research Center, Cambridge, NASA 
TND-3683, 1966. 

**See H. Scheffe’, Analysis of Variance, John Wieley, 1959, p. 14. 


5-15 



In fact g B (X) is a minimum variance estimate for each component of gg. One can show that 


SB® = ( H T Z i 1 + N H ir 1 Hj'^ 1 +N X (5-9) 

Further the ejected value of g-g(X) is gg, even if we have used an incorrect covariance 
matrix + N< The covariance matrix of gg is 

e ««b - Sb> — b - Ib )T] * < H ? e * + n h i> _1 


5. 3. 1.3 Posterior Mean 

In the following discussion we assume that S^t), o(t), n A (r) are gaussian processes. 
Hence the "optimum” estimate of ^ is the posterior mean. This estimate is optimum 
with respect to any loss function L(€) on each component where* 

00 L(0) = 0 

(ii) L(£ 2 ) 2 L(€ j) s o when <= 2 2 € i 2 0 

(m) L(€) = L(-€), 


For example, let gg be an estimate of gg. Then 
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and 


and 
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E 


g Bl" g Bl 

g B2 _g B2 

g B3 _g B3 


*See S. Sherman "Non mean -Square Error Criteria” IRE TRANS. ON INFORMATION 
THEORY, Vol. IT -4, No. 3, p. 125. 
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are all minimized when gg is the posterior mean. As a second example, let £ ^ £ 2 , ^3 
represent the maximum admissible errors in components of gravity. Let 



Then 


and 


and 


E L L i(g B rg B1 )3 = 

E[L 2 <g B2 -g B2 )] = 

E ^ L 3^B3-g B 3) ] = 


p rj|%r%i| > e ij 

P rj|%2"%2| > € 2 J 

P r| | ^B3 - ^B3 | > S j 


are all minimized when gg is the posterior mean. (The expression "Pr {- }" denotes 
"probability that [-} ”.) 


To evaluate the posterior mean, we first determine the conditional distribution of X, given 
gg. In this subsection the notation is the same as that in 5. 3. 1.2. From equation (5-8), 
it follows that X is normally distributed with mean H-^gg and covariance 


From prior observations we have an estimate of orientation of the ISU; and hence we have 
an estimate of Tj, say Tg Also we have an estimate of the magnitude of gg, say g. With 
these estimates, a prior distribution can be defined for gg, namely, gaussian with mean 


/V 

lB 



£ 

0 

0 


and covariance 



"2 2 


0 

0 


~2 2 



where ct o , is the rms error in the estimate of gg|; and a ^ is the rms error in the estimate 
of vertical (expressed in radians). 
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This prior distribution implies that the distribution of gg is isotropic in a horizontal 
plane. 

Applying Bayes formula we find that the posterior distribution of gg, given X, is gaussian 
with mean 


fees (5-11) 

and covariance matrix 

Lg = (H'J’s-^Hj+Bg 1 )- 1 (5-12) 


The estimate gg (X) represents the optimum combination of the measurements X and 
prior data weighted by their respective errors* Note that the covariance is 


If our prior alignment information is poor, the posterior -mean estimate reduces to the 
least-square estimate. Specifically, ascr^-* 00 and erg"*® 5 , then 0 and expression 

(5-11) approaches expression (5-9). Also the covariance (5-12) approaches (5-10). 


The estimate (5-11) and covariance (5-12) are\the basis for an iterative alignment technique. 
Specifically, the initial gg( ^ , si*) are obtained from K measurements based on gp, and 


S~. The second estimate and covariance 2^ are obtained from a second set of K 

measurements based on and etc. This iterative technique is sub-optimal since 
we are summarizing all of the prior measurements in terms of g^ and sQ) . A true 
recursive "least-squares" technique involves significantly more computation since suc- 
cessive measurements are correlated.* Also, all back measurements are used in the 
current computation. The mtermeasurement correlation can be eliminated by augmenting 
the measurement variable. ** This approach also results in a very complex estimation 
procedure. From a practical viewpoint, the sub-optimal technique described above is 
a reasonable compromise. 


si 2 ) 




*P. Gainer, "A Method For Computing the Effect of an Additional Observation on a 
Previous Least-Squares Estimate", NASA Langley Research Center, NASA TN 
D-1599, 1963. 

**See M. Aoki, "Optimization of Stochastic Systems", Academic Press, 1967, p. 38ff. 


5-18 



5.3.2 Estimation, of Instantaneous Components 


In the previous subsection we developed several techniques of estimating the average 
components of acceleration 


T 


1 


g 

0 

0 


In this subsection the objective is to estimate the instantaneous components at time 
t* > KAt; i.e. , 
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T 1 
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+ T! 

~ge 3 (t*) 
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g0 2 (t*) 


In the following discussion we assume that the stochastic inputs 0- (r), a (r), n A (r) are 
gaussian processes. The discussion in Section 5. 3. 1. 3 applies here; the posterior mean 
will be used to estimate instantaneous components. If our prior alignment data is poor, 
the posterior -mean estimate reduces to the least-squares estimate. 

The vectors X, N, 0^ 0 2 , and gg, are defined in Section 5.3. 1. Let S* denote the 
instantaneous components i.e., 
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0 

S* = Tj 

0 

+ T 1 

“g® 3 (t*) 


0 


ge 2 (t*) 


To obtain the conditional distribution of S* given X, we first obtain the joint distribution 
of (S*, X). The components of (S*, X) can be expressed in terms of fundamental random 
variables as follows: 

S * = g B j + < T 1 >j 3S 0 2 - ( T I> 3 2S e 3 * 3 = 1 > 2 > 3 

X i = g B l At + C 1 1) 13**2 W “ (Ti) u fit# 3 (I) + 

X i+K = g B2 At + < T 1>23**2W “ ( T l>22®*3^ + N i+K 

X i+2K = g B3 At + ^ T l ) 33 g0 2^ " ^ T l ) 32 g0 3^ + N i+2K , i = 1, 2,3, •••K 
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or in matrix form as follows: 



T 

where Z is the vector 



Lg B i> gg2’ &B3* N l> 


N. 


0 ; 


* 3K’ 




0 2 (K) } 0*, 




,(K)] 


and where V is the corresponding matrix. Note that V can he evaluated using prior 
estimates g and Tp since they only enter as multipliers of 0g and 0 ^ . A prior distri- 
bution of g B is based on prior alignment data — namely, gaussian with mean 


and eovariance 
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= T. 
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~2 2 
g a 8 
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a similar prior distribution was used in Section 5. 3. 1. 3. The variate Z_ is gaussian 
with mean 

^ i. = ^g]31 5 ®B2’ ®B3’ ***> 

and covariance 
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where 

C Q = E[(9*) 2 ] 

At 

c = E[9|0 2 (])] = J «» e (r - (K+l-j)At - € t )dr 
J 0 

with t* = KAt + € t . 

-bol r I 

If 0g is approximated by b^e , then 

c 0 = b l 

bl b2 C(j - 1) At - t*],. b2 A t _ 

C* “ G L G — 

’ »2 

/ S*\ _ m 

Hence, [ x I is a gaussian variate with mean V Z_ and covariance YS Z V . The condi- 
tional distribution of S* given X is gaussian. + To evaluate the conditional mean and co- 
variance, we must partition the mean and covariance matrix as follows: 

t See T. W. Anderson, "An Introduction to Multivariate Statistical Analysis", John Wiley, 
1958, p. 27ff. 
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" a l" 

and (VS z v T ) _1 = 

A B 
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vz = 

where is 3x1 and A is 3x3. The conditional mean is 

g B = ’■A BX + a ^ + A Bag 

The conditional covariance is A" 1 . Note that 

E C(S* - g B ) (S* - g B ) T ] = A" 1 


(5-13) 


The discussion of iterative techniques in subsection 5. 3„ 1. 3 applies here. The 
estimate (5-13) can be used in such an iterative technique. 

5.3.3 Discussion 

Several level alignment estimation techniques were suggested in this section. A Monte 
Carlo simulation was performed to select the best estimation technique. The simulation 
is discussed in Section 5 of the trade-off report. Three techniques were considered: 
simple average (5. 3. 1. 1), posterior mean (5. 3. 1. 3), and instantaneous estimation 
(5. 3.2. 1). Several values of K, At, and T EB were tried. The effect of nongaussian 
noise was also investigated. The instantaneous estimate is superior to the other estimates, 
in some cases the rms alignment error being one-half the alignment error obtained with 
the simple average. The instantaneous estimate is selected as the recommended technique. 
The simple average is selected as an alternate technique, since it is computationally less 
complex. 

The simulation was also used to investigate the characteristics of the recommended 
estimation techniques. The results of the simulation suggest the following conclusions 
for level alignment: 

8 The instantaneous estimate is probably not sensitive to the noise distribution 
(gaussian or nongaussian). 

• Rotational motion from the environment is most probably the dominant source 
of error for long integration intervals (At > 15 sec). 

• The instantaneous estimate is more accurate than the simple average for 
At 30 sec. 


5-22 



• instantaneous estimation and simple average appear to have comparable accuracy 
for At < 15 sec. 

• If At is held fixed at about 30 seconds and the quantization increased, the 
instantaneous estimate becomes less accurate than the simple average. 

• Low frequency environment noise is not the dominant source of error for short 
integration intervals (At < 15 sec). 

The above points are a summary of the detailed analysis of the simulation results included 
in Section 5 of the trade-off document (Volume 2). 

5„4 ESTIMATION OF GRAVITY AND EARTH RATE IN GYROCOMPASS 

Estimation of the components of gravity and earth rate is based on the observational 
equations (5-4) and (5-5). Using the observed sensor outputs P^(j) and 
] = 1, ••<>, K, we estimate the average components (5-1) or the instantaneous 
components (5-2) and (5-3). The average estimate is investigated in the following 
section, 5„4. 1; the instantaneous estimate is discussed in 5.4.2. The basic estimation 
problem in Gyrocompass Alignment is very similar to estimation in Level Alignment. 

Note that the basic operation is differentiation. We obtain acceleration from velocity 
measurements and angular velocity from angle measurements. 

In the following subsections several estimation methods are developed from a mathematical 
statistics viewpoint. A general discussion of Gyrocompass Alignment techniques appears 
in 5.4. 3, and the recommended technique is presented in 5„ 6. The characteristics of 
two techniques are described in Section 5 of the trade-off document. These techniques 
are simple average and posterior -mean estimate of average components. 

5.4.1 Estimation of Average Components 

The objective here is to estimate the average gravity and earth-rate components in the 
body frame — namely, 



- -i 



E 


g 


“1 * T, 

U 

Ib = T i 

0 

and 

0 


0 



3 


In using an estimate of the average components, we are neglecting the motion about the 

■p 

average. The error bound derived in subsection 5.3, 1 applies to to as well as g. 
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In the following subsections three approaches to estimation are considered: simple average, 
least squares, and posterior mean. The first approach does not use any a prion informa- 
tion about the noise spectra, alignment, gravity, or earth rate. The second approach uses 
prior measurements of the noise spectra. On the other hand, it does not include the prior 
geophysical measurements of gravity and earth angular velocity. The third approach uses 
a priori information about alignment, magnitude of gravity and magnitude of earth angular 
velocity plus prior measurements of the noise spectra. However, the noises must be 
assumed to be gaussian processes. This third approach has several advantages: i) prior 
geophysical measurements are included and are weighted with estimates of their accuracy; 
ii)the estimation techniques are comparable to those obtained from a least squares approach 
m complexity; iii) the resulting techniques can be used recursively to continuously update 
the alignment matrix; iv) the posterior -mean estimate is optimum with respect to a large 
class of loss functions, not ]ust quadratic. From noise simulation, we find that the 
posterior -mean techniques are probably not sensitive to the gaussian assumption (see 
Section 5 of the trade-off document). 

5. 4. 1. 1 Simple Average 


This approach is based on the assumption that we do not have any prior information about 
the noise, magnitude of gravity, magnitude of earth rate, or alignment. In this case. 


A 


1 K 

S 

KAt ]=! 



E 
— B 


1 K 

S 

KAt]=l 




Note that the same estimate is obtained if K = 1 and At is replaced by KAt. 


5. 4. 1. 2 Least Squares 

Before developing a least-squares technique, it is convenient to define certain notation. 
Let X be the 6K vector whose components are 


with 3 = 1, 2, 


Xj = p^o), 

X j+K “ P B2^)> 

X ]+2K ' P B3^ 

j+3K " P 31 l]l -- 

X ]+4K = 

X j+5K = P B3( ] ) 

K. Let Hg be the 6Kx6 matrix 
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j At . 

*£(]) = I s r (T)dT 


with r = 2, 3. 

Let £ be a 6K vector with components 

f = *2> —> *6K> 
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where 


* 3 = g( T l>13 0 2^ - g( T l)l2V ]) 

# ]+K = g( T l)23 0 2^ " g( T 1^22 0 3^ 

$ j+2K = ^^1^33 0 2 ^ " g ^ T 1^32 0 3^^ 

^3+3K = ^-(Ti)ii 4 + t T l)l3 w U^ 0 2^^ “ ( T 1^12 w U 0 3^^ + ( T l)l2 0 2^ 

+ ^ T 1^13 0 3^ > (5-14) 

$ j+4K = *-'^ T 1^21 w N + ( T l)23 a> U'' 0 2^ " ( T 1^22 W U 0 3 ^ + ^ T 1^22 0 2^^ 

+ ^ T 1^2 3 0 3^) 

# j + 5K = 1 '^ T l)31 t0 N + ^ T 1^33 a) U- i 0 2^ ' ^ T l)32 w U 0 3^) + ( T l)32 9 2^ 

+ (T^gg^JU) y 


with j 

= 12** 

• 5 K. 



Let 

*i- 

^ 2 (1)> 

# 2 ( 2 )» 

..., 0 2 (K)3 

and 

(tp T - 

C^(l), 


• • • , 0 2 (K)3 

and 

4 ?- 

^ 3 ( 1 ); 

0 g (2) , 

0 3 (K)3 

and 


: C0 + (1), 

#J(2), 

•••, ^(K)] 


Then equation (5-14) can be rewritten as follows. 



where is the 6Kx4K matrix defined by equation (5-14) and is introduced for mathe- 
matical convenience. Let N be a 6K vector with components 
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N. 


N. 


]+K 


N. 


]+2K 


N j+3K 


N. 


J+4K 


N 3+5K 


jAt 

J la At) + n A1 

(j-l)At 

] At 

I 1* 2 (t) + nA2 

(3-1) At 

] At 

J , Ca 3 < T > + n A 3 

(J-l)At 

3 At 

I n G l( T ) dT 

(i-i)At 

3 At 

J n G 2 (r)dr 

(3-l)At 

] At 

I n G3^ dT 

(3-1) At 


(r)] dr 
(r)3dr 
(t) ] d r 


with 3 = 1, 2, • •♦,]£. The vector N represents the environment and sensor noises. The 
basic observation equations (5-4) and (5-5) can be rewritten as: 

+ £ + N (5-15) 

- - E 

The objective is to estimate gg and uig in the presence of noise + N. 


X = H, 


Ib 


-E 
— B 


The covariance matrix of + N is the sum of the covariance matrices, since the noises 
are independent; i.e. , 

Si „ = S, £. T 

®+N # N 


Further, 


% = ^ 


111 

S 00_]VJ_O_' o_ 

|^b 0 + I ^0+0 +j^ o I o 
JJ^ + 

_° | 0 j s 00+| E 0+0 -*J 


h 
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where 


(*W)i3 = EC0 2 (i)0 2 (])] 
(S 0+0+ ). 3 , E[0 + (i)0+(j)] 

Pvt+hl = E[0 2 (i)«>J(])] 


with i, ] = 1, 2, • ••, K. Note that e 2 and 0 3 are assumed to be independent and 
identically distributed. Further, 



where 


Po+n)!] = ELN^j], i, 3 = 1, K 

PqJjj = ECN^] , i, j = 3K+1, •♦*, 4K 


These covariance matrices can be expressed in terms of the correlation functions (see 
subsection 5.4. 1.2). The following identities* are useful in simplifying £ 0+0+ and 

S 0P+ : 

d 2 

S [0 *(0)e '(T)] = Et0(O) 0(T)] 

dr 2 

d 

E[0(O) 9' (t)] = Et0(O) 0(r)] 

dr 


where we have assumed that 6 2 (t) and 0g(r) are stationary processes. 


*E. Parzen, "Stochastic Processes", Holden-Day, 1952, p.83. 
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Note that is evaluated using prior estimates for the value of and g, denoted by 
T-^ and g. Precise values are not needed since is used in the noise model. Corrections 
to would be of second order. 


Based on the composite measurement equation (5-15), the objective is to find the unbiased 
linear estimate of (gg, Wg), say (j|p(X) , Wg(X)), which minimizes e|§ b -g(X)j 2 as a function 
of g(X) and minimizes e|w|j - w(X)P as a function of w(X). It follows from the Gauss-Markoff 
theorem that (gg(X) , uig(X)) is the value of (g, a)) that mmimizes 


MX - MH, 


w 


rn 

where M is the nonsingular matrix such that N M = I. In fact gg(X) and 


a T? T? 

w g(X) are minimum variance estimates for each component of gg and cog. One can 
show that 


^(X) 


- (hTs ; 1 


3 4>+N n 3 


Ho) 


-1 


-1 


X 


(5-16) 


Further the expected value of gg(X) is gg, and expected value of fig (X) is w®, even if 
we have used an incorrect covariance matrix The covariance matrix of 

(Ib> © 1S 


( H 3 S $+N H 3 ) 


(5-17) 


5. 4. 1. 3 Posterior Mean 


In the following discussion we assume that the stochastic inputs are gaussian processes. 



The "optimum” estimate of gg and cog is then the posterior mean as shown in Section 

5.3. 1.3. 


To evaluate the posterior mean, we first determine the conditional distribution of X 

T? — 

given (gg, Wg). In this subsection the notation is the same as that in 5.4. 1.2. From 
equation (5-15) it follows that X is normally distributed with mean 


£b 


H 3 



and covariance 2^-^ 
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From prior observations we have an estimate of the orientation of the ISU, and hence we 
have an estimate of T, , say T 1 . Also we have an estimate of the magnitude of g„ 

(say g), the magnitude of Wg, (say w ) and latitude (say A.) from geophysical and 
astronomical measurements. With these estimates a prior distribution can be defined 
for gg — namely, gaussian with mean 



and covariance 



where is the rms error in the estimate of j gg j and is the rms error in the 

estimate of vertical (expressed m radians). Similarly, we can define a prior distribution 
- W 

for Wg — namely, gaussian with mean 



and covariance 



■p I 

where cr^ is the rms error in the estimate of t5g | and cr^ is the rms error in the 
estimated direction of Wg (expressed in radians). Let 
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Applying Bayes’ formula, we find that the posterior distribution of ( gg , dg) given X 
is gaussian with mean 


f B (X) 

*1® 

and covariance 

5 - (HjEilNHs+Si 1 )' 1 


* H 3 S $+N H 3 + S T/ f H 3 S 4>+N— + ^ 


-B I 


0) 


:E 

B 


(5-18) 


(5-19) 


"R 1 

The estmiates g(X) and Mg (X) represent the optimum combination of the measurements 
X and prior data weighted by their respective errors. 


The alignment procedure described by (5-18) and (5-19) reduces to a least squares 
procedure when the prior measurements are very inaccurate. Also, this procedure 
can be used recursively to update the alignment matrix (see Section 5. 3. 1. 3). 


5.4.2 Estimation of Instantaneous Components 


The earth's angular velocity is small compared with gyro quantization, in contrast to 
gravity and accelerometer quantization. Hence, it is reasonable to estimate the average 
angular velocity and instantaneous gravity. In Section 5. 3. 3 we concluded that the 
posterior -mean estimate of the instantaneous gravity components is best, based on a 
Monte Carlo simulation. On the other hand, based on the same simulation, there is no 
advantage in using a posterior -mean estimate of earth rate as opposed to a simple 
average (see Section 5.4.3). 


5.4.3 Discussion 


Several alignment estimation techniques are suggested. A Monte Carlo simulation was 
performed to select the best estimation technique. The simulation is discussed in 
Section 5 of the trade-off document. Two techniques were considered - simple average 
(5.4. 1. 1) and posterior mean (5.4. 1.3). Several values of K, At, and T were tried. 
The effect of nongaussian noise was also investigated. The simple average was superior 
to the posterior mean. An alternate technique is to use an instantaneous estimate of 
g (5.3.2) and an average estimate of to . The accuracy will be improved but at the price 
of a significant increase in the computation requirements. Therefore, the recommended 
technique is simple average of both accelerometer and gyro data. 
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The simulation is also used to investigate the characteristics of the recommended 
techniques. The accuracy of gyrocompass alignment is strongly dependent on quantization 
errors of the gyro. The alignment error is of the order of 100 seconds of arc. 

5. 5 CALCULATION OF ALIGNMENT MATRICES FROM ESTIMATES OF GRAVITY, 
EARTH RATE AND OPTICAL ANGLES 

The final operation in alignment is the calculation of the alignment matrix (see Chart 
5-1). The basic equations are developed in Section 2.3.3 and are repeated here for 
completeness. The Mirror -Alignment matrix is presented in Chart 5-3, the Level- 
Alignment matrix in Chart 5-4, and the Gyrocompass matrix in Chart 5-5. 

5.6 RECOMMENDED ALIGNMENT TECHNIQUES 

Referring back to Chart 5-1, we find that there are four basic types of equations: 
alignment matrix, preprocessing, estimation, and estimation matrix equations. The 
alignment matrix equations are presented in Charts 5-3, 5-4, and 5-5 for Mirror 
Alignment, Level Alignment, and Gyrocompass. The alignment matrix computations 
are the only computations needed for Mirror Alignment. The preprocessing equations 
for level alignment and gyrocompass are presented in Chart 5-2 of Section 5. 1. Note 
that the dots on the left indicate which equations are used for Level Alignment and 
Gyrocompass. 

The estimation equations for Level Alignment are presented in Chart 5-6. The estima- 
tion matrix equations are presented in Charts 5-7 and 5-8. The estimation equations 
and matrix equations for Gyrocompass Alignment are presented m Chart 5-9. 

The procedures required to implement the preceding alignment techniques are presented 
in the Procedures Manual, Part 3. The estimation equations were programmed for the 
Monte Carlo simulation, which is described in Section 5 of the trade-off document. 
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CHART 5-3 


MIRROR ALIGNMENT MATRIX 


Inputs 9-^ Qg and 


From these quantities the alignment matrix is given by. 


(y-M^ 

(Mj X U)*(Mi X Mg) 

(E x N)- (M x x Mg) 

iMjL xM 2 l 

iMiXMgl 

(E-Mi) 

(Mi X E)-(Mi X Mg) 

(Nx U)-(Mi XMg) 

|Mi xMgl 

| M 1 x Mg I 

(N-Mj) 

(Mi x N)- (MjX M>) 

(U x E)*(Mi x M 2 ) 

|Mi xMg l 

|Mi x Mg | 


where 

IJ^x M 2 | = Cl - (M x - Mg) 2 ] 1 / 2 

(Mj_ • M 2 ) = (M^ U)(M 2 . U) + (M^ E)(Mg. E) + (M x . N)(Mg- N) 


(u-Mi) 


COS ©i 


(u-m 2 ) 


cos 0g 

(E-Mi) 

= 

cos o^i sin0i 


(e-m 2 ) 

= 

cos^g sin0g 

(N-M^ 


sm i sin 0 i 

7 

(N-Mg) 


sin ®g sin 0g 


An optional technique might utilize the value of lM^ x Mg I from a previous alignment 
thus eliminating the aforementioned dot product and square root operations. 






CHART 5-4 


LEVEL ALIGNMENT MATRIX 
Inputs (g- B^, (g ° B 2 ), (g- B 3 ) and 

From these quantities the alignment matrix is given by: 




1 

o 

o 

1 


0 

T 

- 

0 

sin u ^ cos a ^ 


0 



0 

-cosa^ sina^ 


1 




iMjXUi 


1 0 


I 

O 

o 

rH 

1 

1 



0 Im^uI 


(U-Bj) (U-Bg) (U*B 3 ) 

(Mi-u) - 1 “ 



iMjxul 0 


0 -(U-Bg) (U-Bj) 


where 

• (Mj- U) = (U- b x ) 

• IMjxUl = Cl - (M-j^-U 2 ] 1 ^ 2 

c (U- Bjj) - fe* B k )/g 

• g = C(g- B^ 2 + (g. B 2 ) 2 + (g • B 3 ) 2 ] l/2 

An optional technique might utilize any of the following additional inputs: 

• The zenith angle (0^) of mirror one might be utilized to find (M^ • U) from 

(Mi* U) = cos 0 X 

• The magmtude of gravity (g) might be supplied from a local survey. This piece 
of information can be utilized to reduce the number of required accelerometers 
to two. 
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CHART 5- 5 



GYROCOMPASS MATRIX 


Inputs (g° (g. R 2 ) j (g • B 2 ) ; (« E - B^) (oj E . B^ and (^ E ° Bg) 

From these quantities the alignment matrix is given by: 


- - 


0 

1 

0 


(W-Bj) 

(W-Bg) 

(W-Bg) 

T 

= 

0 

0 

1 

|wxu[ 


(U*Bi) 

(U'B 2 ) 

(U-Bg) 



1 

(W-U) 








|WxUl 

|WxU| 

0 


(WxUj-^xB 

3 ) (WxUMBgXB^ (WxUMBjXBg) 


where 

(Wo u) = (w.BjXu.B^ + (w- b 2 )(u- b 2 ) + (w° b 3 )(u. b 3 ) 

|WxU| = Cl - (W-U) 2 ] 1 ^ 

(W.B^ = (w E . B^/oF 
(U-y = (g-B^/g 

« E = C (w E ° B x ) 2 + (w E - B 2 ) 2 + (co E . Bg) 2 ] 1 / 2 
g = [(g. B 1 ) 2 + (g- B 2 ) 2 + (g. B 3 ) 2 ] 1//2 

An optional technique might utilize any of the following additional inputs: 

• The local latitude (X) might be utilized to find (W • U) from 

(W . U) = cos X 

• The magnitude |of gravity (g) might be supplied from a local survey. 

• The magnitude of earth rate j(w E ) | might be supplied from a local survey. 

A use of all additional inputs could reduce the number of necessary instruments to three 
(either two accelerometers and one gyro, or one accelerometer and two gyros). 
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CHART 5-6 


ESTIMATION ROUTINE COMPUTATIONS - LEVEL 


Inputs: Preprocessed accelerometer measurements, X, estimation matrix, M, 
and vector, b 

/\ 

Output Estimate of acceleration components in body frame, g° B , 1 = 1, 2, 3 
at time t* 


MX + b 


KAt At 

••• J a‘B.,dt,j' a-Bgdt,-**, 
(K-l)At 0 

a • Bg dt 


• Posterior Mean Technique (Instantaneous): Computations of b and M from 
the Estimation Matrix Computation Chart 

• Simple Average Technique' 


The basic estimation computation is 


r g% (t*) 

A, (t*) 

g-B 3 (t^) 


where 


X T 


'At 2 At 

| a-B dt, | a-B-dt, 
0 At 


KAt JKAt 

J a • B 2 dt , • • • j 

(K-l)At (K-l) At 


At = Intersample time 
K = Number of samples 



K 


K 


K 
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ESTIMATION MATRIX COMPUTATIONS - LEVEL 

Inputs: intersample time, At (sec) 

number of samples, K 

~ 9 

estimate of gravity, g (ft/sec ) 

9 

rms error in gravity estimate a (ft/sec ) 

o 

estimate of T^ B , T-^ 

rms angular error in prior estimate of vertical, (radians) 
noise covariance functions (tabular) 

• accelerometer noise $ n (t) (ft /sec ) 

2 4 

• translational acceleration noise 0 a (t) (ft /sec ) 

2 

® rotational noise 0g (t) (radian ) 
prediction time (sec) 

Outputs alignment parameters M and b 

The intermediate quantities E n , E^, Cq, c^. and Eg are computed from the 
inputs. 

• is Kx K matrix with components 

At 

(S a )y = J* [At- u]0 a (u + (j-i)At)da 
„0 

+ J [At + u] 0Ju + (]-i)At)du 
-At 

• S n is K x K matrix with components 

At 

(E ) =| [At - u]0 (u + (]-i)At)du 

11 0 

0 

+ J [At + u] 0_(u + (j-i)At)du 
-At n 
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CONTINUATION OF CHART 5-7 


• isaKxK matrix with components 


(S )■ = J [At - u]0 Q (u+ (j-i)At)du 


+ J [At + u] 0 e { u+ (j-i) At) du 
-At 


• c g = 0Q (°) 


c- = f ( ui-5- (j-l)At - KAt - € f )da 

] o 


1 = 1, 2, •••, K 


where the integrals are evaluated by a convenient integration technique such as 
trapezoidal rule or Simpson's rule. 

• S~ is a 3 x 3 matrix 


m “2 2 mT 

T 1 g a 0 T 1 

n ~2 2 

L° S 


From these intermediate quantities, £ z , V, A, B, , and ag are computed. 


• Ejj is a 3K x 3K matrix 


£ +2 
n a 


£ + S 
n a 


£ + £ 
n a 
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CONTINUATION OF CHART 5-7 


• £ z is a (5K + 5) x (5K + 5) matrix 


%l i 0 

; 


s z “ 


"N i ^ 

Tfj 


i c oQ”*_ C |^_ 


"i c oFr*- c K _ 

T c .i| 

| : ! 

|C K 


3 3K 1 K 


Numbers at edges of matrices denote dimension of submatrices. 


• Matrices A (3 x 3) and B (3 x 3K) are submatrices 


A B ’mi 

T = (V^V 7 )' 1 

b a d 6 


where matrix V is the (3K + 3) x (5K -f 5) matrix given on the following chart. 


i(Ti>i 

g(T n )« 


^tg(Tj) n ^ 
A tg(T 1 )j 1 \ 

AtgCT^i \ 
Atg(T x ) 31 | 

At g 1)31 ^ 
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CONTINUATION OF CHART 5-7 


Then, the outputs are given by 

• M = -A -1 B 

_i 

Note that A is the covariance matrix of the estimate. 

• b = + A~ 1 Ba 2 
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CHART 5-8 



CHART 5-9 


ESTIMATION ROUTINE COMPUTATIONS - GYROCOMPASS 


Inputs: Preprocessed accelerometer measurements, X, estimation matrix, M, 
and vector, b 


Outputs. Estimates of gravity and earth rate components m body frame, g • IT and 
i = 1, 2, 3 


The basic estimation computation is 


g-Bi 


E R 

co • 


E ^ 
co • B 9 

“ -ft 


= MX + b 


where 


At KAt At KAt 

J* a- B^ dt, • • • . J a* B-^dt, J* a* B 2 dt, • • • , J* a* Bgdt, 

0 (K-l)At 0 (K-l)At 


At KAt KAt 

J* co • dt , * * • , J co • B^ dt , • • • , J co • B^ dt 

0 (K-l) At (K-l)At 


At 

K 


Intersample time 
Number of samples 
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CONTINUATION OF CHART 5-9 


• Simple Average Technique; 


b = 0 



K K K K K K 

Diagonal 6 x 6K 


• Hybrid Technique; Use posterior -mean estimate of gravity as given for Level Alignment and simple 

average for earth rate. 








APPENDIX A 

THE MATHEMATICAL MODEL OF THE 
VIBRATING STRING ACCELEROMETER 


A -1.0 INTRODUCTION 

This appendix describes the operation of the Vibrating String Accelerometer (VSA) that 
has been selected for the ERC Strapdown Inertial Guidance System and develops a 
mathematical model to be used to relate the output of this type instrument to an estimate 
of applied acceleration. 

A-2.0 DESCRIPTION OF THE ACCELEROMETER'S OPERATION 

A functional block diagram of the Vibrating String Accelerometer (VSA) is shown in 
Figure A-l. The accelerometer consists of a seismic mass (mass 1 and mass 2 separated 
by a spring) which is supported by. 1) two taut strings that function as oscillator "tank" 
circuits, and 2) ligaments as shown in Figure A-l and normal to the plane of Figure A-l. 

When the VSA is at rest or moving with constant velocity, the sum of forces acting on its 
seismic mass is zero. When the VSA is accelerated, the resultant force acting on the 
seismic mass changes so that it accelerates with the case. The displacement of the 
seismic mass, relative to the case, that is produced by this resultant force is negligible 
except along the sensitive axis, A, as shown in Figure A-l. The tension in the strings as 
a result will not be affected by any motion other than that along the sensitive axis. This 
change in tension (from the at rest tension) of each string is, therefore, a function of the 
acceleration acting along the A axis of the instrument. 

Since the natural frequency of a vibrating -string is a function of its tension, the vibrating 
frequencies of the strings in the accelerometer are directly related to the applied accel- 
eration along A. 

Each of the strings of the VSA passes through a magnetic field supplied by the two permanent 
magnets of Figure A-l. When set to vibrating in its field, an electric signal is generated 
by the string. This signal is amplified and fed back to the string in such a manner that a 
sustained vibration occurs. The electric signals so generated are nominally sinusoidal 
with frequency equal to the resonant frequency of the individual string. The vibrating string 
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Model: ARMA D4E Vibrating String Accelerometer 

Axis: A is a unit vector directed along strings S-j and S 9 

(the sensitive axis) 


Figure A-l. A Schematic Diagram of the Accelerometer 







acts as a high Q tank cricuit for the oscillator within the associated feedback amplifier 
electronics. The vibrating frequency of each string is read by using Schmidt triggers to 
generate pulses corresponding to zero crossings of the respective sine waves. The 
frequency of the zero crossing pulse train is proportional to the frequency of the vibrating 
string. 

A-3.0 KINETICS OF VSA 
A-3. 1 COORDINATE AXES 

The accelerometer coordinate axes (A, O, P) used in the derivation of the fundamental 
mathematical model are illustrated in Figure A-2. The unit vector A is along the nominal 
position of the string, while O and P are unit vectors arbitrarily defined to make A, O, 
and P a right handed, orthogonal system. 

A-3. 2 THE TENSION IN THE STRING 

The forces acting on and Mg along A are shown in Figure A-3. The lateral supporting 
forces along P (normal to the page) are not shown in this figure. 

The equations of motion for the two masses can be written 

Ff - pF u = Mja ; F 2 = E F 2j = M 2 a (A-l) 

As we are interested in the tension of the strings, only the A. component of equation A- 1 
will be considered. 


(pyA) = Mjfe • A) 
(F 2 - A) = M 2 la • A ) 


(A-2) 


Because the supportmg force from the ligaments acts orthogonally to the string, we have 


(FfA) = Tg-T 1 
£ 2 ‘A) = T 2 -T 3 


(A-3) 
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Therefore from equations A-2 and A-3, we have 


T 3 -T 1 = M^a-A) 
T 2 -T 3 — M 2 (^* * A ) 


which can be combined as 


T 2 -T x = (T 2 -T 3 ) n-CTg-Tj) = (Mj + M 2 ) (a • A) (A -4) 

When the accelerometer is stationary a = 0 and equation A-4 gives T 2 = T x . Let this tension 
in the strings be defined as T . When the VSA experiences an acceleration a, T 1 and T 2 
will be changed to cause the seismic masses to accelerate with the case. If (a • A_) is 
positive (T 2 -T 1 ) will also be positive. The tension in each of the strings can be written as 


T 2 = T o + AT 2 


T 


1 


T o + AT l 


(A— 5) 


and 


AT 1 + AT 2 = (Mj + M 2 ) (a • A) 


The amount of change in the tension of both strings will be the same if the strings are 
identical. In practice, however, the strings cannot be made to be identical. 

For the range of accelerations that is within the proportional limit of the strings, we can 
write 


AT X = Kj(a * A) 
AT 2 = K^a* A) 


(A -6) 


where 


K^ + Kg — ^2 



Beyond that range of applied acceleration, the strings creep and the mathematical model 
derived m the following pages will not apply. 


A -4.0 MODEL DEVELOPMENT 

A -4. 1 THE RESONANT FREQUENCY OF A VIBRATING STRING 

The resonant frequency of a uniform string under tension is directly proportional to the 
square root of the tension. 

For the vibrating strings S^ and Sg, their pulse tram frequencies f ^ and fg (which are pro- 
portional to the respective resonant frequencies) can be written as 

f i = c iV^~ 

(A-7) 

f 2 = C 2V^ 

C j and C 2 are proportionality constants determined by the dimensions, density and other 
physical properties of the strings. Combining equations A-5, A-6 and A-7, we have 

f l = c l 1 /T 0 -K I (a.A) 

( 2 = C 2 V T o +K 2(i-^ 


By Taylor's expansion, this becomes 
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The above expansions will converge rapidly since the instrument is constructed so that T Q 
is much greater than either K^(a* A) and KgCa • A). The frequency difference of the 
two strings is 


f 2 " f l - ( C 2 ~ C l)^o + 2 + C l K l)te*A) “ 


8T 


o^: 


(C 2 k| - CjK^Xa'A)* 


1 b 

+ (C 2 k| + CjKjXa'A) 3 (C 2 k| - C^Xa -A) 4 • • • (A-8) 

iftft 128 T o#o 


The series given in equation A-8 converges rapidly because T Q is made large. However, 
Cj, C 2 , Kj and are constants determined by the dimension and material of the vibrating 
strings S x and S 2 . The accelerometer is manufactured so that (C^-Cg) and (K^-Kg) are 
kept as small as possible (a highly symmetric instrument). For this reason, the even 
order terms are very small and the linear term is the most significant. The approximated 
frequency difference obtained by truncating equation A-8 after the third degree terms may 
be written: 


f 2" f l =D 1 D 0 +D l(i'^:) +D 1 D 2 (fL* A) 2 +D x D 3 (a. A) 3 


(A-9) 


where, by definition, 


D 1 D 0= ( C 2- C l>l/^ 
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A -4. 2 THE ACCELEROMETER READOUT 


Let the signal generated by the string S be 


t 

e^t) = sin[©j + J t irf^dt] 

where f is the frequency of the pulse train from string S (in pulses per second) and 0 
is a constant. The number of zero crossings in the interval (t , t b ) is 

t K - A t 9 .tb 

N > f 6 1 fdt = Eq] + 5 f,dt 

‘a +At l t a 

where (t & + At^) and (t b - At 2 ) are the times of the first and last zero crossings in the 
interval (t a , t b ) and Eq ^ is the quantization error given by 


t + At-. t 

Eli - J a fjdt - J l.dt 

‘a VS* 


The time increments At^ and At 2 are defined by illustration in Figure A-4. In the same 
way, for string S 2 , we have 




Figure A- 4. At^ and Atg 


A- 9 


The difference in the number of zero crossings of the two strings can be written, using 
equation A -9, as 


t. 

N 2 - N x = (Eq 1 - Eq 2 ) + J °(f 2 - f j) dt 




= Eq + J D 1 D Q dt+ J* D 1 (a-A)dt 


+ J tb D 1 D 2 (a-A) 2 dt+ J tb D 1 D 3 (a«A) 3 dt 


A-4. 3 THE VSA FUNDAMENTAL MATHEMATICAL MODEL 

In summary, the readout of the VSA is two pulse trams corresponding to the zero crossings 
of the sinusoidal signals from the two vibrating strings. The input to the VSA is the 
acceleration of the case along its sensitive axis. The accelerometer readout is related 
to its input by the mathematical model given in Chart A-l. 
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CHART A- 1 


THE FUNDAMENTAL ACCELEROMETER MODEL 


THE ACCELEROMETER MODEL IS: 


t b *b 

J f 2 dt - J f l dt = ( N 2' N l) + Eq = D x J (a-A)dt 


f> , « I 

+ D x U CD 0 +D 2 (a-A)' 1 +D 3 (a<A) ,:i ]dt > 
< *a ^ 


a is the acceleration applied to the accelerometer 

t fit £ is the time interval over which a is measured 

A is a unit vector directed along the input axis of the accelerometer 

and N 2 are the number of zero crossings detected in t^tst^ 
from both strings of the accelerometer 

Eq is the instrument quantization error due to the fact that t and t^ 
do not correspond to zero crossings 

D^ is the accelerometer scale factor 

Dq is the accelerometer bias 

D 2 is the second order coefficient 

Dg is the third order coefficient 

f 2 and f^ are string frequencies in pulses/second 
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APPENDIX B 

THE MATHEMATICAL MODEL OF THE GYROSCOPE 


B-1.0 INTRODUCTION 

The purpose of this appendix is to find the mathematical expression that relates the outputs 
of the Honeywell GG 334A gyroscope to the environmental input to which this gyro is 
subjected. 

Section B-2.0 is devoted to the general description of the 334A gyro and its principle of 
operation. The mathematical model is then developed in Section B-3. 0. 

B-2.0 DESCRIPTION OF THE GYRO OPERATION 

The Honeywell GG 334A gyro contains a gimballed rotor spinning at a very high angular 
rate (see Figure B-l). A hydrodynamic gas bearing is used to support the rotor. The 
gimbal is restricted by the gimbal bearing to rotate only about the output axis relative to 
the case. The signal generator of Figure B-l consists of a moving coil attached to the 
gimbal and a stationary wound stator attached to the gyro case. It generates an a-c voltage 
with an amplitude that is directly proportional to the angular displacement of the moving 
coil from its null position. In this way the gimbal deflection relative to the case is 
measured. At each sampling cycle (3 0 6 KHz rate), the gimbal deflection is detected, 
sampled and compared to two thresholds (positive and negative of equal level) to determine 
if a positive, zero or negative rebalance torque is to be generated. A current switch and 
associated electronics provide the torque generator with correct torquing current pulses 
of constant strength. The timing information (3. 6 KHz) used to derive the cycle periods 
is furnished. 

Any angular motion of the gyro case about the input axis, G, will generate a gyroscopic 
torque that tends to rotate the gimbal about the output axis, O. The signal generator 
senses the resulting gimbal deflection and produces the signal to the gyro electronics 
necessary to generate the correct torquing current pulses to the torque generator. In 
this way, the gyroscopic torque developed initially about the gimbal axis is rebalanced by 
the pulsive torque produced by the torque generator. The average rebalance torque is 
proportional to the average gyroscopic torque which is in turn proportional to the gyro 
angular rate about G. A readout of the pulse train of the rebalance current is used as the 
instrument's output. 
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Model, Honeywell GG 334A single -degree-of -freedom, pulse 
rebalance gyroscope. 

Axes: S is a unit vector along the spin axis of the rotor. 

O is a unit vector along the output axis as defined by 
the gnnbal. 

G = O x S is the sensitive axis of the gyro. 


Figure B-l. A Schematic Diagram of the Gyro 










B-3. 0 ]VIATH MODEL DEVELOPMENT 


B-3. 1 COORDINATE AXES 


The coordinate system used in the following derivations is illustrated in Figure B-2. The 
gimbal axes (G_, O, S) are defined as fixed to the gyro gimbal with O (the output axis) 
directed along the gimbal rotary axis, S (the spin axis) directed along the gyro rotor spin 
axis and G (the input axis) directed along the direction ofOxS. 

The set (G, O, S) is right handed, orthogonal and is assumed to be coincident with the 
gimbal principal axes. 

B-3. 2 THE GIMBAL DYNAMICS 

The gimbal angular momentum can be expressed in the gimbal coordinate axes as 

H = Cl GG ( w g • G) JG + [I Q0 ( w g • O) JO + Cl gs (w g • S) + H r ] S (B-l) 


where I GG is the moment of inertia of the gimbal and the rotor about G, Iqq is the 
moment of inertia of the gimbal and the rotor about O, Igg is the moment of inertia of 
the gimbal about S , and H r is the constant rotor spinning angular momentum. (p g -G), 
(w g * O), and (cu g *S) are the components of the gimbal angular velocity about G , O, and 
S, respectively. 


Since G, O, S are assumed to be the prmcipal axes of 


inertias I 


GO ; 


: OS 5 ^S 5 


etc. are assumed to be zero. 


states that the torque applied to the gimbal is equal to 
momentum. 


the gimbal, all the products of 
The second law of rotational motion 
the derivative of the gimbal angular 


dll 

dt 


Using equation B-l and writing equation B-2 in component form, we have 


(T-G) = I GG (^ S *G) + (w g -p)(« g .S)(I ss -I 00 )4-H r (a- g .O) 


(3-2) 


(B-3) 
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(T-0) = I 00 (^ S '0) + (« g -S) (cU g .G) (I GG -I ss )~H r (w g .G) 
( 1 • S) = Igg(w § * S) + (u> g - G) (w g - O) (I OQ - I GG ) 


(B-3) 


(T*G) and (T-_S) are reaction torques from the gimbal bearing. Since we are concerned, 
only with the torque about the output axis, O, we may write 

T q = (T-O) -I 00 (^ S .0) -H r (oi g .G) + (to g S)(tu g -G)(I GG -I gs ) 

For a single -degree-of -freedom gyro, the gimbal can only move relative to the case about 
the gimbal axis O. Thus the gimbal angular velocity can be expressed in terms of the case 
angular velocity and the relative angular velocity between the gimbal and the case. 

u> g = + (u) g ~co) =w + 80 

where oj is the gyro case angular velocity and 9 is the gimbal deflection with respect to 
the case (see Figure B-3), Therefore, 

T o = (T-0) = I 00 9 +I 00 fe*0) -H r (w-G) +-( i o-S)(w-G)(3 GG -I ss ) (B-4) 

B-3. 3 THE GIMBAL TORQUE 

The gimbal torque, T q , is the sum of all torques applied to the gimbal about the output 
axis, O. T q includes a dampening torque, T^; a rebalance torque, T , provided by the 
torque generator; and error torques. 

* 

The dampening torque is proportional to the rate of change of gimbal deflection angle, 9 . 

T d »-ce 


The rebalance torque is 


T 


r 


= -L5 k 


where 6 k is the logic value of the pulse at the instant t k . 
negative or zero pulses, respectively. 


& k = +T, -1, or 0 for positive. 
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L 


L 0 LU(t -t k ) - U(t -t k - h)J 


where L q is the amplitude of the pulse 
h is the pulse width 
U{t) xs the unit step function 
ty, is the k 01 sampling period. 


L o h = 


A * 1 


Ldt is the strength of the pulse and would be constant for a linear rebalance 


loop. However, the torque rebalance loop will not be linear in reality, and L Q h will be a 
function of u>. To take into account the effect of nonlinearity, let us assume 


L o h = L 0 h(l+a(to.G t )) 


where L Q h is constant and a very small. 

2 

The error torques include a constant torque, an a-sensitive torque, an a -sensitive torque, 
and other torques considered as noise (for example, reaction torque from the signal 
generator) . 

The constant torque is denoted by R'. 

The a-sensitive torques are mainly due to the fact that the center of support of the gimbal 
is not coincident with its center of mass. If the gyro acceleration is a, the a-sensitive 
error torque is 


T 1 = B| (a- G) + B^(a- O) + B'(a.S) 


where Bp Bq and Bg are gyro unbalance coefficients, 
o 

The a -sensitive error torque is due to the fact that the gyro gimbal is not a rigid body. 

To make the gimbal follow the motion of the gyro case, there are forces acting on the 

gimbal through the gnnbal bearing. The gimbal deforms when subjected to these forces. 

Because of this deformation, the center of mass of the gimbal will be displaced from the 

q 

center of support and therefore produce an a -sensitive torque about the output axis. It 
is assumed that the deformations also occur in the lateral direction as well as along the 
direction of the acceleration. 
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If the acceleration is a, the a -sensitive error torque is 

T 2 = C fj(a-G) 2 + Cgg(a-S) 2 + C^(a-G) (a-S) + C^ s (a-Q)(a.S) + C| Q (a. G)(a- O) 
where C^, C^g, C^g , C G g , and Cjq are so-called compliance coefficients. 

The total gimbal torque is then 

T q = -C0-L6 k +R’ + Bj (a- G)+B^(a- q)+B£ (a- S)+C^.(a- G) 2 +Cg S (a- S) 2 + C^ g (a- G) (a- S) 


+ C 0s(-“ 0)(a- S) +C > Q (a- G)(a- O) + T n 


(B-5) 


where T n is the torque due to other effects and is considered as a noise component. 


B-3.4 CONCLUSION 

Combining equation B-4 and equation B-5 we have 

-ce - L 6 k +R '+Bj (a.’ G) +B^(a-0)+B^(a.S)+C^(a.G) 2 +Cg S (a.S) 2 +C| g (a.G)(a. S) 

+ Cq S ( a. 0)(a* S)+C ' Q (a- G)(a- 0)+T a = I 00 (w-0) - H r (w-G)+(ui.S)(co .G)(I GG -I gs ) 

where & is the gyro case angular velocity, and the component of the reaction torque, 
*00^’ has been neglected since it is small compared with the damping torque, CQ . 

A rearrangement of equation B-6 gives 

— 6 , ( 0 ,.G) + -hrR'tBJ(a.G)tBJ ) (a.O)*B^a.S) + C : f I (a.G) 2 

H r H r 


+ c ss + + 


(I GG ‘ *<?) - W“*2>] 


T n Cf5 


H r H r 


(B-7) 
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Let us now integrate equation B-7 over the IN sampling periods starting at time t and 
ending at time t n . We have 


L o h 


li 


( E 6 - J (tu • G) dt + J R + Bj (a • G) + Bq (a *^) + Bg (a • S) 
k=l + L ” 


+ Cjj(a*G)^ + C gs (a • S) 2 + Cjg (a » G)(a « S) + CQg (a » Q)(a • S) 

+ C lO^-‘— )(-'—) +Qi S («J‘G)(w.S) + J(co .oj|dt + A n + Eq (B-8) 
where co is the gyro case angular velocity 


j-t- 

= 1,0, -1 is the k m rebalance pulse 
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pj pi p*t pi p t 

u ss u is u os °io 




, , , , and , respectively 





H r 3 H r H r H r 


“ (I GG " W 


J = 


*00 


EL 


f *n Tn 

An = j dt is the effect of noise torques 

to H r 


„t„ C 


p n * 

Eq = J — 6 dt is a quantization error, 
to H r 
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However, 


L Q h = L Q h[ 1 + a(cu*G)J 

With ol very small, we can write 

1 

L h = L h ** ( i - (co • G)) L h 

i + «(fe-G) ° 

Multiplying both sides of equation B-8 by (1 - -G)) produces the gyi’o model as given 

on. Chart B-l. 

L o h V 

In the model, = is the instrument scale factor. ©„ = -a is the coefficient 

H r 11 

of the term of the scale factor nonlinearity, a is very small, and the higher order effects 
of it have been ignored. 
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CHART B-l 



WHERE 


• u- is the angular velocity applied to the gyro 

• a is the acceleration applied to the gyro 

• tg £ t £ t N is the time interval over which a and a- are measured 

• - tg = Nr, where N is an integer, andT is the gyro sampling 
period 

• S is a unit vector along the spin axis of the rotor 

• O is a unit vector directed along the output axis as defined by the 
gimbal 

• G is a unit vector along OxS (that is, the sensitive axis of the gyro) 

• 6^ is the kth gyro pulse, equal to +1, -1, or 0 for positive, negative, 
or no pulse 

• A<t> is the gyro scale factor 

• R is the gyro bias 

• Bj Bq and Bg are the gyro unbalance coefficients 

• ^II Sis ^IS ^OS anc ^ ^IO are gy ro compliance coefficients 

• Qjg and Qjj are dynamic coupling coefficients due to gimbal deflection 
and scale factor nonlinearity, respectively 

• J is the angular rate coefficient 

• An is the effect of gyro noise over the [tg, tj^] interval 

• Eq is the gyro quantization error 






APPENDIX C 

ALTERNATE FORM OF Q MATRICES 


In Section 2.2.4 the Q G and Q A matrices were expressed in terms of (Q A ) ” ^ and 

(Q G )~\ which are the matrices calibrated in the ERC laboratory. The calibrated elements 

were seen in Section 2.2.4 to have either the form • B ^ or G^ * B? . Because the body 

axes (B^) are defined, the elements do not equate directly to physical ISU angles — angles 

like the angle between, say, two gyro axes. It is possible, however, to express the Q 

matrices as a function of physical angles only. Those expressions are found in Chart C-l 

and C-2. In Chart C-I we see the general expression, and in Chart C-2 we see the first 

G A 

order approximation of the matrices. (Recall that the nominal Q and Q matrices are 
identity matrices.) 

The form of the two matrices (Q A and Q G ) in Chart C-l and C-2 are, naturally, the same. 
In Chart C-l the Q matrices have been separated into sums and products of submatrices, 
where each submatrix is a function of only one type of ISU angle. For example, the first 
submatrix is a function of only the angle between the mirrors, the second submatrix is a 
function of only the angles between the mirrors and instruments, the third submatrix is a 
function of only the angles between the accelerometers or gyros, and so forth. 

The calibrated Q A and Q G elements can be equated to the elements found-m Chart C-2, 
allowing for the solution of the physical ISU angles. Such solutions could be useful for 
the determination of the satisfaction of design requirements. 
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CHART C-l 
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CHART C-2 


Q MATRICES (OPERATIONAL) 


a N _ N _ N 

h ~ Si 

a N _ ,,N _N 

—2 -2 “ —2 

N N 

M ^ be perpendicular to M 2 


where the superscript N denotes 
the nominal vector 


IGNORING SECOND ORDER TERMS, WE HAVE- 


r 1 (Mj-Ag-A^Ag) (Mj-Ag-Ag.A^ 

Q A = (Mg-Aj - M r M 2 - Aj-Ag) 1 (M 2 -A 3 - A 2 -A 3 ) 

-(M^Ag) -(“2*^) 1 

Or, letting A^) = M^M^A^) 

and A^l^-Ag) = 
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APPENDIX D 

COMPUTER SYSTEM DESCRIPTION 


This appendix contains a description of the laboratory computer and its associated equip- 
ment. Section 1 describes the laboratory computer. Section 2 describes the Interface 
Electronics Unit (IEU), the device that interfaces the computer to the ISU. Section 3 
describes briefly the devices used for the computer manual interface. 

D-1.0 COMPUTER 

The Honeywell DDP-124 computer is a small scale scientific/control digital computer 
with a 1.75 microsecond memory access time. The memory is an 8, 192 word, 24-bit/ 
word core memory. Arithmetic is performed on 24-bit sign- magnitude (not complement) 
data with the left-most bit of the data word containing the sign and the other 23 bits con- 
taming a binary representation of the magnitude. The basic arithmetic register is a 24- 
bit A Register which is extended by a 24-bit B Register for multiplication, division and 
shifting. 

The instruction repertoire contains 47 instructions allowing fairly flexible fixed point 
processing. Unique instructions include a step multiple precision, store address portion 
of A, output and input to A (may be ANDs) as well as input/output to memory, direct con- 
trol pulse outputs and sense line skips. Because of the sign magnitude number representa- 
tion, the computer has both arithmetic shifts (sign bit(s) do not shift) and logical shifts. 

Indirect addressing may be performed by use of one bit in the instruction Three index 
registers are available. 

A Fortran IV Compiler is available and is considered preferable by NASA for calibration 

» > 

programming. The 124 is not equipped with floating point hardware so use of the Fortran 
Compiler will necessitate use of time consuming floating point software routines Because 
of real-time considerations, Fortran shall not be used for alignment 

The computer interfaces with the Interface Electronics Unit, the displays and magnetic tape 
unit via a direct memory access (DMA) subunit. This allows direct transfer of data from 
and to memory under buffer control in one of two modes. These modes are the time sharing 
mode and the hog mode. In the hog mode, the input/output will hold the memory until the 
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entire transfer is completed. In the time sharmg mode the input/output and the processor share 
memory with either locked out for one memory cycle while the other completes one transfer. 

Execution times of instructions are as given in the DPP- 124 Programmers Reference Manual, 
and execution times and memory sizing for standard arithmetic subroutines are as given m 
the DDP-124 Users Guide. One magnetic tape handler is available for program storage 
and/or other uses. 

D-2.0 INTERFACE ELECTRONIC UNIT (IEU) 

The IEU provides the computer an interface to the system equipment. A block diagram of 
the IEU is shown in Figure D-l . 

The IEU counts information in its counters from the gyros, the accelerometer strings and 
the timer. Each counter is compared to a manually selected interrupt condition. This 
condition Is selected as any number for the time counters or any power of two for gyro 
and accelerometer inputs. When an interrupt condition is met, a signal is sent to inter- 
rupt logic 3 if time counter 2 has satisfied it condition, interrupt logic 2 if time counter 
1 has satisfied its condition or interrupt logic 1 if a gyro or accelerometer register has 
satisfied its condition. The interrupt logic generates an interrupt to the computer on its 
own interrupt channel and sends a reset signal to the counters Interrupt logic 1 and 2 
send reset signals to all of the counters other than time counter 2 and interrupt logic 3 
sends a reset signal only to time counter 2. When a counter receives a reset signal, it 
will hold the contents of the main register, clear an auxiliary register and begin to ac- 
cumulate data in the auxiliary register. 

When the computer has received an interrupt, it will initiate a direct memory access 
(DMA) controlled input from the IEU of the counters and ISU status registers. The main 
registers of the counters are read. After the reading process has been completed, a 
resume signal will be sent to all the counters from the DMA input control. This signal 
will cause any counter that is counting in an auxiliary register to clear the main register, 
add the auxiliary register to the main register, and continue accumulating in the main 
register 

While the IEU has the capability of using any of the inputs to determine sampling rate as 
described above, it is not expected that any criterion other than time counter 1 is needed 
for the main calibration and alignment routines. The IEU interface program should verify 
that the time criterion has been met (Interrupt i). 
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Figure D-l. IEU Block Diagram 
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Since the IEU sampling determines the time to sample ISU data on the basis of a time 
criterion, use of the IEU results in the maximum worst-case quantization error. 

The IEU will transfer up to four 24-bit parallel data words from (and to) die ISU to (and 
from) tiie computer. Output from the computer is via DMA transfer. 

D-3. 0 COMPUTER MANUAL INTERFACE DEVICES 

The operator interfaces with the computer via the display panel, a keyboard and typewriter 
and a paper tape reader and punch. 

The display panel can display nine numbers. Each number has a signed one decimal digit 
mantissa and a signed five decimal digit characteristic. This capability will be used to 
display results or intermediate results or request and to display normalized data output 
from the ISU during real-time data collection by the computer 

The display panel has three rows of eight buttons each to be used to select parameters to 
be displayed and 24 buttons to select program options. 

The keyboard and typewriter may be used to enter data into the computer in small amounts 
and to furnish the operator with information such as desired settings of test table axes. 
Maximum transfer rate is 15. 5 characters per second. 

The paper tape reader and punch will be used to enter large amounts of data into the com- 
puter and for output of the results of the procedure Maximum transfer rates are 110 
6-bit characters per second for the reader and 300 6-bit characters per second for the 
punch. 


ENb 
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